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CD;0D and CD;OH have been prepared and their infra-red spectrum in the range from 2.54 
to 184 has been measured. The Raman spectrum of CD;OD has been observed. CH;OH and 
CH;0D have been reinvestigated and an assignment for all four methyl alcohols has been given. 


ETHYL alcohol and methyl alcohol-d have 
been subjected to a thorough investigation 
in the Raman! as well as in the infra-red spec- 
trum.? Methyl alcohol can be treated theoreti- 
cally since it is a comparatively simple molecule. 
One degree of freedom involves the torsion about 
a single bond. From an accurate assignment of 
all the other fundamental frequencies it should 
be possible to evaluate the value of the torsional 
vibration with the use of heat capacity data. 
The hindrance potential to free rotation around 
this single bond could thus be determined. 
Methyl alcohol is an asymmetric top molecule 
of only small asymmetry, thus differing only 
slightly from a molecule of the point group C3». 
Molecules of this group have a class of doubly 
degenerate frequencies. The degeneracy would 
naturally be lifted, if the asymmetry of the 
molecule had any effect on the methyl group of 
the molecule. Then the highest symmetry of the 
methyl alcohols would be C,. The splitting of 
these frequencies therefore gives a measure of 
the asymmet lecule. 


Chem. Ph 
(1938) Borden and E. F. Barker, J. Chem. Phys. 6, 553 


To obtain more data the isotopic molecules 
CD;0D and CD;OH were investigated in the 
infra-red. CD;0D was also investigated in the 
Raman spectrum. 


EXPERIMENTAL PART 


The infra-red spectrometer used has been 
described before.* The polished rocksalt windows 
of the absorption cells (30 cm long, 2.5 cm in 
diameter) were sealed on with Apiezon wax. 
Blanks were recorded over the whole range. 
The alcohol was kept in a side arm and was 
sealed off to avoid as far as possible any contact 
with air and stopcock grease which tend to 
exchange with the deuterium of the hydroxyl 
group. Later on, when different fractions of the 
alcohol were measured to obtain a change in 
relative intensities of CD;0D and impurity 
bands, the system was closed by a stopcock. 
This caused some exchange in the OD group. 
The vapor pressure was regulated by immersing 
the side arm into cooling mixtures. The Raman 
spectrum was taken with a cell containing 
~0.4 cc CD;0D. On some plates the 4047A 


3H. Gershinowitz and E. B. Wilson, Jr., J. Chem. Phys. 
6, 194 (1938). 


693 


en 
nt 
he 
ti- | 
to 
Id- 
als 
ow | 
1ey 
not 
on- 
ent 
1in, 
ne- 
sly. 
ant, 
| 
the 
‘hey 
lace 
| 


HERMANN D. NOETHER 


% TRANSMISSION 
80 


1700 1800 1900 2000 2100 2200 2300 2400 2600 2800 30003100 31003100 4000 
FREQUENCY IN CM” 


1350 1400 


Fic. 1. Absorption bands in CD;OD. Curves (a) to (g): These curves show mainly the 
bands of the impurity (bands at ~580, 836, 933, 1069, 1160, 2056, 2196, 2925, and 3000 
cm“). (a) First fraction (0°C); (6) same, CD,0D had been stored in a °300-cc flask and 
only the gas phase was condensed in the trap of the absorption cell. Thus mostly the im- 
purity was collected. (c) (~—50°C); (d) high ressure run of the sample, collected as 
described under (5); (e) (—50°C); (f) (—32°C); (g) sample as under (6) high pressure run. 
Note the slight absorption in the OD stretch range; there is no absorption in the range 
from 3100 to 4100 cm. Curves (hk) to (p): middle and third fraction containing the im- 
purity with bands at 1260, 1370, and 1762 cm™. (hk) (O°C); (z) (—8°C); (Rk) (5-10°C ?); 
(2) later record with some CD;0H present; (m) (0°C) early record with nearly no CD;0OH in 
sample; (m) sat. pressure; (0) (15°C ?); (p) sat. pressure ?; © (0°C). (7) to (v): third fraction 
and nearly pure CD;0D with some CD;OH. (r) (—1°C); (s) low pressure nearly pure 
CD;0D; (¢) 11°C); (u) (12°C); (v) (~25°C) some CD OH present. The temperatures 


given are those of the ‘cooling mixtures in which the cell trap was immersed. 


excitation was eliminated by a Wratten 2A filter, 
while a Rhodamine 5GDN Extra (DuPont)! 
solution cut down the green mercury lines. 

CD;,0D was prepared from CD;Br by hy- 
drolysis with a D.,O solution of NaOD on the 
steam bath in a sealed tube. The D.O and 
CD;OD were separated by fractional distillation. 
CD;OH was prepared from CD;0D by exchange 
with H:O followed by fractional vacuum distil- 
lation of the CD;OH. CH;0D was prepared 
by reacting D.,O (99.6 percent) with dry 
Mg(OCHs)2.° The CH;OH used in this prepara- 
tion and for some infra-red measurements was 
chemically pure (Merck) and had been distilled 
from Mg(CH;0). after refluxing for several 
hours. 


_4J. T. Edsall and E. B. Wilson, Jr., J. Chem. Phys. 6, 
12 (1938). 
5 See reference 1. 


Frequencies of CD,;0D and CD,0H 


The frequencies of the absorption maxima 
measured for CD;OD in the range from 2.5y to 
18u are given in Table I. However, many of 


TABLE I. Infra-red and Raman frequencies in the CD;0D 
spectrum (frequencies of impurities in parentheses). 


Raman 


(576) m (6) 
~821 v.w. (4) 

989 s (9) 
1072 v.w. (7) 
1131 v.w. (2-3) 
(1168) v.w. (2-3) 


2082 s (10) 
2162 m (3) 
2245 br (8) 

In brackets num- 
ber of times ob- 
served. / and k 
excitation. 


Infra-red 
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Fic. 2. Absorption bands of CD;0H. Curves: 750-1400 cm™ range: 
(a) (~—35°C) from exchange of first fraction of CD;0D; (6) same fraction 
at 25°C; (c) from exchange of third fraction of CD;O0D (—14°C) (strong 
change in the intensity of the 1060 and 1165 cm= bands); (d) third fraction 
(4°C); (e) third fraction (25°C); (f) third fraction (—10°C); (g) third 
fraction (0°C). There is no absorption from 1400 to 1700 cm™. Range 


1700-4100 cm: (a) second fraction (25°C); (6) third fraction (—3°C); 
note the change in relative intensities in the 2100 cm™ bands; (c) second 


fraction (~15°C); (d) first fraction (25°C). The first two fractions contain 
a considerable amount of the impurity with bands at ~580, 833, 933, 
1060, 1165, 2062, 2186, and 2888 cm~. The third fraction contains only a 
small amount of it. 


these apparently belong to two impurities CD,;OD and Fig. 2 for CD;OH). Therefore it is 


which were difficult to eliminate completely 
because of the small quantity of material avail- 
able. Certain bands were strong in the first 
distillation fraction (either plain vacuum distil- 
lation or vacuum distillation from charcoal),® 
weak in the middle fraction, and practically 
absent in the last fraction. (See Fig. 1 for 


TABLE IT. Frequencies in the infra-red spectrum of CD;0H 
(frequencies of impurities in parentheses). 


1270 s 
1289 s 
1308 s 
(1355) v.w. 
(1767) w 
1975 w 
(2062) s 
2094 s 
(2186) s 
2232 s 
(2888) m 
(2940) m 


(1165) s ~3680 m 


°H. D. Noether, J. Chem. Phys. 10, 664 (1942). 


fairly certain that these bands belong to a more 
volatile impurity. This is substantiated by low 
temperature runs in which only the bands of 
this impurity could be observed. They are at 
~580, 836, 933, 1069, 1160, 2056, 2196, 2925, 
and 3000 cm~. This assignment is confirmed by 
the fact that when these bands are strong the 
987 and 2725 cm~ bands which quite definitely 
belong to CD;0D are weak or missing. The 
bands at 1290 and 3700 cm also occur in the 
spectrum of CD;OH (see below) and are only 
found in later records of CD;0D where some 
exchange in the OD group had taken place. 
The bands of the other impurity which is present 
in the last fraction of the charcoal distillation of 
CD;0D and CD;OH had already been deter- 
mined and separated in the measurement of the 
spectrum of CD;Cl and CD;Br.* These bands 
occurred originally quite strongly at 1263, 1362, 
and 1760 cm~, but became much weaker after 
the charcoal treatment. 
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FREQUENCY IN CM’ 
Fic. 3. Infra-red bands of CHsOH and CH;OD. Curves: (a) CH;OD (20°C); (6) CH;OD 


(- 16°C); (¢c) CHs3OH (—22°C) the curve is shifted upwards by 10 


percent to avoid super- 


position with curve (6); (d) CHsOH (19°C); (e) CH;OD (20°C); CH;0H (4°C); (g) 
CH;OH (12°C). Temperatures i in brackets are those of the cooling mixtures in which the 


cell trap was immersed. 


For the CD;OH measurements (Fig. 2, Table 
II), three different fractions of CD;0D were ex- 
changed with H.O. The first fraction gave strong 
bands at 580, 833, 1060, 1165, 2062, 2186, and 
2900 cm=. These are exactly the bands which 
were found to belong to the impurity in CD,OD. 
The CO band at 990 and the OH band at 3680 
are medium strong, proving the presence of 
CD;OH. In the other exchanged fractions the 
same shifts in intensities of the bands occur as 
already observed in CD;OD, thus supporting 
the conclusions drawn for the bands of the 
CD;0D spectrum. The second impurity is 
present only in quite small concentration as 
shown by the weakness of its 1760 band. 


Frequencies of CH;0H and CH;0D 


The infra-red spectra of CH;0H and CH;OD 
were measured to find the exact shape and 
intensities of the fundamental bands of these 
compounds in our spectrograph (see Tables III 
and IV, Fig. 3). The bands agree in general very 
well with those measured by Borden and Barker? 
and Barker and Bosschieter.’? There are a few 
very weak bands in the 1230 region which have 
not been reported by these authors. They are 
at 1209 and 1260 cm for CH;OH and at 1207 
and 1232 for CH;OD. They were taken at 
saturation vapor pressure at 20°C. Comparing 
the 1450 range of CH;0H and CH;OD, it could 
be assumed that there is another weak band 
around 1430 cm which is partly covered by the 
very strong 1340 absorption. This band is quite 


7E. F. Barker and G. Bosschieter, J. Chem. Phys. 6, 
563 (1938). 


clear in the CH;0D molecule where the two 
bands of Barker and Bosschieter at 1459 and 
1480 cm! seem to form one strong band at 
1468 cm due to our insufficient resolution. 
The intensities of the bands in CH;OH and 
CH;0D are as follows: The 1033 and 1040 bands 
have exactly the same shape, the R branch 
being stronger than P and Q. The 863 band is 
very strong, and the 1340 band is much stronger 
than those in the 1460 range. The CH stretches 
are strong and have about the same intensity. 


Assignments 

We can expect 12 frequencies for the methyl 
alcohol molecule of which 3 pairs should be 
either doubly degenerate, if the molecule behaves 
like a symmetrical top molecule, or should show 
a slight splitting, if the OH or the OD group 
makes the molecule asymmetric. The degenerate 
or nearly degenerate frequencies represent the 
carbon-hydrogen stretches, bendings, and rock- 
ings of the methyl group, all connected with a 
change of electric moment perpendicular to the 
C axis. Then there are three frequencies for 
which the electric moment changes parallel to 


TABLE III. Infra-red frequencies of CH;OH. 


Own 
measurements Borden & Barker 


1455 
1477 


2054 


2845 
2978 
3683 


Own 
measurements Borden & Barker 


s 1034 
v.s 1058 
1209 w 
1260 w 
P 1322 
Q 1340 }s 
R 


1365 
? 1425-1430 


P 1014 1010 
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the C axis, representing the symmetrical stretch- 
ing and bending of the carbon hydrogen bonds 
and the carbon-oxygen stretch. They should 
show the structure of a parallel band with P, Q, 
and R branches. There are finally three vibra- 
tions of the hydroxyl group: the oxygen- 
hydrogen stretch, the oxygen-hydrogen bending, 
which because of the nearly rectangular angle of 
the oxygen hydrogen bond with the C axis of the 
molecule should give rise to a parallel type band, 
and the torsion of the hydroxyl against the 
methyl group, where the hydroxyl hydrogen 
vibrates perpendicular to the OH(OD) bond and 
the C axis. At higher energies this vibration 
changes to a rotation of the hydroxyl group 
around the CO axis. This frequency is very low 
and probably undetectable with our instrument. 


(a) C—D and C—H Stretches 


The assignment of the observed bands to the 
fundamental frequencies is given in Table V. 
v(r) and »(c) are easily assigned by comparing 
with the corresponding frequencies of the ordi- 
nary and the deutero-methyl! halides. In these 
compounds two parallel type bands are found in 
this region due to resonance of 276(¢) with v(x). 
In CH;OH and CH;0D only one parallel type 
band is observed in this range. Because of the 
impurity band at 2190 in the CDs alcohols, the 
presence of the resonance doublet cannot be 
ascertained. In the Raman spectrum of CH;O0H 
and CH,OD, Halford, Anderson, and Kissin! 
find the resonance doublet 24(c), v(x) at 2837, 
2942 and 2882, 2948, respectively. In the Raman 
spectrum of CD;0D this doublet occurs at 2082 
and 2162 


(b) Methyl Group Bendings 


From a comparison with the methyl halides 
and especially with CH;F, Barker assigns the 


TABLE IV. Infra-red frequencies of CH;OD. 


Barker& 


measurements Bosschieter 


1459 
1480 


2065 


2720 
2850 
2965 
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parallel type bending 6(7) for CH;O0H and 
CH;0D as given in Table V. The assignment 
of the 1120 frequency in CD;0D and CD;OH 
as 6() is based on the following considerations: 
5(3) increases from 1252 for CHs3I, to 1305 for 
CH;Br, 1355 for CH;Cl, 1455 for CH;OH, and 
1459 for CH;0D. CD;Br and CD;Cl have at 
987 and 1029 cm~. From CH;Cl to CH;OH the 
frequency increases by ca. 100 cm~'. Expecting a 
similar change for the deutero series on the basis 


TABLE V. Frequency assignment of the methyl 
alcohols (in 


Characteristic of band 


(I) 
C—H stretch ( 1 ) 
bending (||) 


H—C—H bendings ( 


CH;OH CH;0D CD;0H CD;,0D 


C—O stretch (||) 
CHs-rockings (_L ) 
¥Q—H(o) O-—H stretch (L) 
Torsion 


O-—H bending 


* The difference in shape of 59 44(+) for CHsOH from 69_y4(x) for the other 
methylalcohols is perhaps due to overlapping by 7’’(o). 


of the nearly equal bromide-chloride change in 
the two series, we obtain the assigned value for 
in CD;0D and CD;OH. For the perpendicu- 
lar type bands in CH;OH and CH;0D the bands 
at 1477 and 1480 can be assigned, but this 
would assume the molecules to behave like 
symmetrical top molecules and this band to be 
doubly degenerate. In CH;SH, which is probably 
less asymmetric than CH;OH, Thompson and 
Skerrett® find a splitting of about 45 cm 
(1430, 1475 cm) for these frequencies. In 
CH;0D a band at 1427 and approximately 
1480 could account for the splitting. In CH;OH 
it is hard to find an exact position for the lower 
part of this doublet since the 1340 band overlaps 
very strongly. From comparison with CH,OD, 
it could be assumed to be somewhere around 
1430, the other part being at 1477 cm~. In the 
methyl halides 6(¢) remains nearly constant from 
to CH;F (1445 to 1476 cm-).® In CD3Br 


5H. W. Thompson and N. P. Skerrett, Trans. Faraday 
Soc. 36, 812 (1940). 
(1928). H. Bennett and C. F. Meyer, Phys. Rev. 32, 888 


|| 
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wo 2978 2965 2232 2297 
d 1455 1459 «1123-1121 
in ~1430 1427 71055 71045 
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at P~1010 ~1012 967 969 
1040 © 990 S987 
1058 ~1060 1009 1005 
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and CD;Cl it is found and calculated from the 
product rule” at 1053 and ~1058 cm-, respec- 
tively, and it can be expected at ca. 1070 cm~! 
for the alcohols. Unfortunately, the impurity in 
CD;0D and CD;OH has a strong band at 1060 
cm, and this makes a definite assignment 
impossible. The flat part of the absorption curve 
in this range for nearly pure CD;OD suggests 
absorption bands in this region, and also a band 
is found at 1072 cm in the Raman spectrum of 
CD;0D which could account for 6(¢). There is 
some evidence that these vibrations may have 
the frequencies indicated in Table V. By the 


-ratio rule" 6’(¢) and for CD;0D and 


CD;0OH are calculated at ~ 1074 and 1036 cm“. 
These values are probably not very far from the 
correct frequencies. The lower part of the 
doublet may form the shoulder of the strong 
990 band of CD;0D and CD;0OH or come 
directly on top of the extremely strong R branch 
at 1007 and 1009 cm-, respectively. 

The assignment 6(7)>46(c) for the CD; 
alcohols is supported by calculations of these 
frequencies for the CD3X series with different 
masses X, using Linnett’s! set of force constants, 
which yields the result that for small X 6(2) >(c). 
Intensity considerations in the methyl halide 
series where is always stronger than 
polarization measurements on CD;NO:,"* where 
the higher frequency is polarized, thus proving 
the assignment and Stitt’s'* assignment of the 
bendings in C2D., where the parallel bendings 
are both higher than the perpendicular bendings, 
support the deutero-methyl alcohol assignment 
of >6(c). 


(c) The Methyl Rocking Frequencies 


In the methyl halides the rocking frequencies 
change from 885 for CH;I to 957, 1021, and 
1200 cm! for the bromide, chloride, and fluoride. 
The alcohols are very similar to the fluoride, 
and therefore these fundamentals could be 
expected somewhat higher than 1200 cm. 
Borden and Barker assign the very strong 1340 
band to the methyl rocking in CH;OH. An 
objection to this assignment is the fact that the 

10 E. B. Wilson, Jr., J. Chem. Phys. 9, 76 (1941). 

11H. D. Noether, J. Chem. Phys. (to be published). 

12]. W. Linnett, J. Chem. Phys. 8, 91 (1940). 


83 T, P, Wilson, J. Chem. Phys. (to be published). 
4 F, Stitt, J. Chem. Phys. 7, 297 (1939). 
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rocking frequencies are usually not very strong, 
being at least weaker than the bending fre- 
quencies. Furthermore, the 1340 band cannot 
be found in CH;OD in spite of its great intensity 
in CH;OH. We suggest that the weak bands at 
1209, 1260 in CH;0H and 1207 and 1230 in 
CH;0D be assigned to rocking frequencies. 
There is considerable uncertainty about the 
exact position of these bands, at least in the 
case of CH;OH, because they are probably 
strongly overlapped by the 1340 and to a lesser 
extent by the 1034 bands. It is possible that in 
CH;0OH the higher frequency of the doublet 
comes on top of the 1340 band. 

A comparison with the rocking frequencies of 
the halides and deutero-halides places these 
motions of the CD; alcohols in the neighborhood 
of 900 cm~ and thus leads to the assignment of 
and for CD;OH as given in Table V. 
The rocking frequencies of CD;OD are probably 
in the same range, but they cannot be found 
because of the strong absorption at 775 and 
987 cm and the impurity band at 833 cm. 


(d) The C—O Stretch 


This band is the strongest band in the alcohol. 
It is a parallel type band showing P, Q, and R 
branches. In all four bands, the R: branch is the 
strongest. The first harmonic is found at 2054, 
2065, 1975, and 1966 cm~, respectively. In 
the Raman spectrum yc-o(r) appears at 1029 
(CH;OH), 1032 (CH;OD), and 989 (CD;OD). 


(e) The Hydroxyl Frequencies 


The hydroxyl stretches are easily assigned on 
comparison with other molecules containing the 
OH or OD group. 

The hydroxyl bendings should show bands of 
parallel type since the vibration is nearly parallel 
to the C axis. The assignment of frequencies for 
hydroxyl bendings is in general quite difficult due 
to association of molecules with hydroxyl groups 
and due to the complexity of the spectrum of 
such molecules in the bending range. For 
isotopic molecules where only the hydrogen of 
the hydroxyl group can be exchanged by 
deuterium, a criterion for a correct hydroxyl 
bending assignment is the disappearance of this 
band on exchanging the hydrogen and the 
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TABLE VI. Doublet spacings for the methyl 
alcohols (in cm~), 


CH;0H CH;0D CD;0H 


48 43 
46 38 
45.9 44.6 


appearance of a very similarly shaped band in 
a lower frequency range. With this rule, it was 
possible to obtain the hydroxyl bending fre- 
quency for the isotopic methyl alcohols. The 
spacings agree quite well with the spacings of 
vc-o(r) (Table VI). They are calculated, using 
Dennison and Gerhard’s formula,!®> which is 
strictly correct only for symmetrical top 
molecules. 


Av=(S/n)(kT/A)! 
logio S=0.72/(8+4)"" 
B=[(A/C)—-1] 
C-—H=1.09A, C—O=1.42A, O-—H=0.96A, 


where 


C—O-H angle=105°. 


The above discussion leads to an assignment 
on a purely experimental basis. A theoretical 
check could be obtained by the application of 
the product rule, treating the molecules as of 
symmetry C, (C—O bond is assumed to be 
rigid, thus splitting out the torsional motion in 
the molecule). It is found that the shifts in 
frequencies observed cannot be accounted for by 
this rule. Applying the ratio rule, the frequencies 
responsible for the failure are found to be the 
hydroxyl bending 59—1(7) and the methyl rocking 


(1933) L. Gerhard and D. M. Dennison, Phys. Rev. 43, 197 


RAMAN SPECTRA 699 
frequencies r(c). It should be noted that the 
hydroxyl bending frequency is reduced by more 
than v2 on substitution of deuterium. Because 
of the failure of the product rule a theoretical 
calculation on the basis of a quadratic potential 
function becomes impossible, at least in the 
approximation in which the torsional motion is 
split out. 

It is interesting to note that only frequencies 
which represent motions influenced by the 
hydroxyl group show disagreement with the rule. 
A very strong shift in the hydroxyl stretching 
motion is usually observed in measurements on 
liquids. This is due to association. Here, however, 
in the measurements of the methyl alcohols in 
the gas phase, the O—H stretch occurs in its 
normal place, thus excluding association as a 
basis for the discrepancy. It is possible that the 
lack of agreement with the product rule can be 
explained by a difference in the amount of 
coupling between the hydroxyl bending and the 
torsional motion for the isotopic species. This 
would simultaneously explain the failure of the 
rule in the case of the methyl rocking, because 
this motion is certainly influenced by the 
torsional motion. A solution for these questions 
can be given only when more data on the 
isotopic hydroxyl groups become available. 
Unfortunately, heat capacity data for the iso- 
topic methyl alcohols cannot be used for the 
solution of the problem, because the torsional 
motion has to be known accurately first, to 
permit any conclusions concerning the hydroxyl 
bendings. 

In conclusion I should like to thank Professor 
E. Bright Wilson, Jr. for suggesting this work 
and for his continued interest and for many 
discussions during the course of this study. 
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IV. The Absorption and Fluorescence Spectra of ms-Tetra(3’,4’-methylenedioxyphenyl)porphine of 
HCl Number Four and of Its AgC;H;N, Zn, and Ni Complex Salts 


V. M. Avpers, H. V. Knorr, AND Davin L. Fry*t 


JOURNAL OF CHEMICAL PHYSICS 


Spectroscopic Studies of the Simpler Porphyrins 


VOLUME 10 


The C. F. Kettering Foundation for the Study of Chlorophyl! and Photosynthesis, Antioch College, Yellow Springs, Ohio 


(Received September 4, 1942) 


The molecular absorption coefficients of ms-tetra(3’,4’-methylenedioxyphenyl) porphine and 


the AgC;H;N, Zn, and Ni complex salts have been measured over the visible region of the 
spectrum. The curves representing the absorption coefficients as a function of the wave-length 
are shown. The spectra of this porphyrin and its metal complex salts are quite similar to those 
of ms-tetraphenylporphine and its corresponding metal complex salts. The fluorescence spectra 
of these substances have also been studied. All four substances fluoresce in the red region of 
the visible spectrum. Curves showing relative intensity of fluorescence as a function of wave- 
length are shown for the parent substance and the Zn complex salt. The fluorescence of the 


AgC;H;N and Ni complex salts is so weak that intensity measurements could not be made. 


INTRODUCTION 


HE preceding publication! of this series de- 
scribed the absorption spectra of ms-tetra- 
phenylporphine** and a series of its metal com- 
plex salts. The successful synthesis of ms-tetra- 
(3’,4’-methylenedioxyphenyl)porphine, Fig. 1, 
and some of its metal complex salts by Rothe- 
mund and Londergan,? at this laboratory, made 
it possible to study the effect, on the absorption 
spectra, of substitution in the phenyl groups 
which occupy the meso positions of the porphine 
ring. The fact that the zinc and nickel com- 
plex salts of both the ms-tetraphenylporphine 
and the ms-tetra(3’,4’-methylenedioxypheny])- 
porphine were available made it possible to com- 
pare not only the spectra of the parent substances 
but also the spectra of two corresponding metal 
complex salts of these parent substances. The 
silver salt of the ms-tetraphenylporphine con- 
tained two atoms of silver while the silver salt 


*From the dissertation submitted by David L. Fry to 
the faculty of the graduate school of The Ohio State Uni- 
versity in partial fulfillment of the requirements for the 
degree of Master of Science. 

Now at the General Motors Research Laboratory, 
Detroit, Michigan. 

1V. M. Albers and H. V. Knorr, J. Chem. Phys. 9, 497- 
502 (1941). 

** The system of nomenclature, for these substances, 
suggested by Rothemund will be followed in this paper, i.e., 
ms-tetraphenylporphine is the porphine with phenyl groups 
substituted in the a, 8, y, and 6 positions in the porphine 


ring. 
(1940) Londergan, Ph.D. Thesis, The Ohio State University 
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of the ms-tetra(3’,4’-methylenedioxypheny])por- 
phine contained one atom of silver and one mole- 
cule of pyridine. The nature of the linkage of the 
pyridine molecule with the porphine ring is not 
known. 

The formula for ms-tetra(3’,4’-methylenedi- 
oxyphenyl) porphine is shown in Fig. 1. The metal 
complex salts are formed by substituting the 
metal for the hydrogens which are attached to 
the nitrogen atoms in the pyrrole rings. Only one 
of the isomers of the ms-tetra(3’,4’-methylene- 
dioxyphenyl)porphine, that of HCI number four, 
was present in the preparations. 


H 
g 
cH 
HC. CH 
¢ 
Cc 
pac 
C-O 
WAN 
H Cc H 
Ho “CH 
HC, 
He 


Fic. 1. Structural formula of ms-tetra(3’,4’- 
methylenedioxypheny]) porphine. 
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Fic. 2. Absorption spectrum of ms-tetra (3’,4’-methylenedioxypheny])- 
porphine in benzene solution. 


EXPERIMENTAL 


The absorption spectra were measured with a 
recording spectroradiometer which has been pre- 
viously described." The spectra were examined 
in all cases from 900 my to 370 mu. Since none of 
the substances showed characteristic absorption 
above 700 mu, only the region between 700 mz 
and 370 my was actually measured. The spectra 
were all measured for dilute benzene solutions 
and the absorption coefficients were calculated 
from Beer’s law in the form of Eq. (1), 


(1) 


where I7/Io is the ratio of transmission of the 
absorption cell filled with solution to its trans- 


*V. M. Albers and H. V. Knorr, J. Opt. Soc. Am. 28, 121- 
123 (1938). 


mission when filled with the pure solvent, c is the 
concentration of the solution in moles per liter, 
x is the thickness of the absorption cell in cm, 
and a is the molecular absorption coefficient of 
the substance. 

The value of a@ was determined at 2.5-myz 
intervals between 700 mu and 370 mu. A con- 
centration of from 150 to 300 mg per liter was 
sufficient to determine @ in the regions of low 
absorption, and a series of dilutions was then 
made in order to obtain concentrations suitable 
for measurement of a in the regions of greater 
absorption. Values of a for each region of the 
spectrum were determined for concentrations of 
solutions which yielded values of I7/Io between 
0.15 and 0.50. Whenever the values of I7/Io for 
two concentrations of solution fell within this 
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Fic. 3. Absorption spectrum of the ms-tetra (3’, 
_ porphine Zn complex salt in benzene solution. 


range at any one wave-length, both were used in 
determining the value of a and a mean of the two 
values was used. 

Figures 2—5 show the absorption spectra of the 
substances studied. In all cases the absorption in 
the violet end of the spectrum is much greater 
than in the rest of the region so the curves are 
broken at the point where this heavy absorption 
sets in and the strong band is plotted on a one- 
twentieth fold scale. The ordinates for the part 
of the curves to the left of the break are indicated 
at the left edge of the graphs and those for the 


part of the curves to the right of the break are 


indicated at the right edge of the graphs. 

The fluorescence spectra were measured in 
benzene solution made up to a concentration of 
10 mg per liter in anhydrous benzene. The solu- 


‘-methylenedioxypheny])- 


tions were placed in a sample tube under an 
atmosphere of nitrogen and irradiated by a spe- 
cial illuminator containing four Pyrex mercury 
arcs arranged parallel to the sample tube. The 
fluorescence radiation emitted from the flat end 
of the sample tube was photographed with a 
Hilger E3 spectrograph. The experimental ar- 
rangement is the same as that used previously 
by the authors* * in the study of the fluorescence 
spectra of the chlorophylls and some of their 
derivatives. The fluorescence spectra were photo- 
graphed on Eastman I-S spectroscopic plates. 
These plates have their maximum sensitivity at 


4V. M. Albers and H. V. Knorr, Cold Spring Harbor 
Symposia on Quantitative Biology 3, 87-97 (1935). 
5H. V. Knorr and V. M. Albers, Cold Spring Harbor 
Symposia on Quantitative Biology 3, 98-107 (1935). 
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Fic. 4. Absorption spectrum of the ms-tetra (3’,4’-methylenedioxyphenyl)- 
porphine Ni complex salt in benzene solution. 


about 690 mu. In order to determine relative in- 
tensity of fluorescence as a function of wave- 
length, the plates were calibrated and standard- 
ized by the method described by Harrison.* The 
lamp used for the standardization was a standard 
lamp operated with its filament at a color tem- 
perature of 2393°K and the calibration was per- 
formed by the use of a series of neutral screens 
varying in transmission from 74.0 percent to 
10.4 percent. The distribution of energy in the 
spectrum of the lamp filament was calculated by 
the use of Planck’s radiation law. The exposure 
times for the calibrating exposures were five 
minutes. On each plate there were thirteen cali- 
bration exposures, a five-minute and a ten- 


°G. R. Harrison, J. Opt. Soc. Am. 19, 267-316 (1929). 


minute exposure to the fluorescence of the parent 
compound, a ten-minute and a twenty-minute 
exposure to the fluorescence of the zinc complex 
salt, a twenty-minute and a forty-minute ex- 
posure to the fluorescence of the silver and also 
of the nickel complex salt. The five-minute ex- 
posure to the fluorescence of the parent com- 
pound and the ten-minute exposure to the fluor- 
escence of the zinc complex salt yielded sufficient 
darkening to obtain microphotometer traces from 
which intensity measurements could be made. 
Experience has shown that, for this type of photo- 
graphic plate, one is justified in assuming that 
the reciprocity law holds for the range of ex- 
posure times used. The relative intensities for 
the spectra of the two substances were reduced 
to the same value for the regions of highest in- 
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Fic. 5. Absorption spectrum of the ms-tetra (3’,4'-methylenedioxyphenyl)- 
porphine AgC;H;N complex salt in benzene solution. 


tensity. Actually, the intensity of fluorescence at 
the position of maximum intensity for the parent 
substance is about four times as great as the in- 
tensity at the maximum of the strongest band of 
the zinc complex salt. The fluorescence of the 
nickel and silver complex salts is so weak that it 
would be impossible to obtain exposures with 
sufficient darkening on the plate without making 
the exposure times so long that there would be 
danger of considerable photo-decomposition in 
the sample. The positions of the darkening on 
the plate produced by the longer exposures to the 
fluorescence of the silver and nickel complex salts 
indicate that their fluorescence spectra consist of 
one band only in about the same region as that 
of the parent compound. The graphs of relative 
intensity of fluorescence as a function of wave- 
length for the fluorescence spectra of the parent 


substance and the zinc complex salt are shown in 
Figs. 6 and 7. 
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Fic. 6. Fluorescence spectrum of ms-tetra(3’,4’-methylene- 
dioxyphenyl)porphine in benzene solution. 
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SPECTROSCOPIC STUDIES OF PORPHYRINS 


DISCUSSION 


In the study of the absorption spectra of ms- 
tetraphenylporphine and its metal complex salts,! 
it was demonstrated that the variations in the 
spectra of the different metal salts of this par- 
ticular ms-tetrasubstituted porphine are much 
greater than the variations which had been pre- 
viously observed’ in the absorption spectra of 
different ms-tetrasubstituted porphine parent 
substances. On the basis of these observations one 
would expect relatively little difference between 


TABLE I. 


Wave-length 
mu 


Name 


ms-Tetra(3’,4’-methylenedioxy- 
phenyl) porphine 


ms-Tetra(3’,4’-methylenedioxy- 
phenyl)porphine Zn complex salt 


ms-Tetra(3’,4’-methylenedioxy- 
phenyl) porphine Ni complex salt 


ou 


ms-Tetra(3’,4’-methylenedioxy- 
phenyl) porphine AgC;H;N com- 
‘plex salt 


resS 
nwo’ 


~ 
~ 


*ms-Tetraphenylporphine 


*ms-Tetraphenylporphine Zn com- 
plex salt 


co 
SERSx 


sos 
Dw 


*ms-Tetraphenylporphine Ni com- 
plex salt 


Hwan’ 


*ms-Tetraphenylporphine Ag2com- 
plex salt 


SE 
on. 


* These data are reproduced from the table in a previous publication 
(reference 1). 


7V. M. Albers and H. V. Knorr, J. Chem. Phys, 4, 422- 
425 (1936). 
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Fic. 7. Fluorescence spectrum of the ms-tetra(3’,4’- 
methylenedioxyphenyl)porphine Zn complex salt in ben- 
zene solution. 


the spectrum of ms-tetraphenylporphine and ms- 
tetra(3’,4’-methylenedioxyphenyl)porphine, and 
one would also expect corresponding metal salts 
of these two compounds to have similar spectra. 
These expectations are supported by the data 
obtained. In Table I the wave-lengths and ab- 
sorption coefficients of all of the bands in the 
spectra of ms-tetra(3’,4’-methylenedioxypheny])- 
porphine and its zinc, nickel, and silver complex 
salts between 700 my and 370 muy are tabulated. 
In addition, the corresponding data for ms-tetra- 
phenylporphine and its zinc, nickel, and silver 
complex salts are reproduced from the table in 
the previous paper.' One of the striking similari- 
ties is that between the spectra of the silver salts 
of the two porphines. One would expect that, 
since the silver salt of the ms-tetra(3’,4’-methyl- 
enedioxyphenyl)porphine contains one molecule 
of pyridine per molecule, this might have an 
effect on the absorption spectrum. This retention 
of pyridine by the porphyrin when pyridine is 
used as a solvent in one stage of the preparation 
is a common occurrence. The fact that the pyri- 
dine is retained in a definite molecular proportion 
even after recrystallization from a chloroform- 
methanol mixture indicates some definite linkage 
of the molecules but this linkage must be of such 
nature that it does not appreciably affect the 
absorption spectrum. 
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The differences in the vapor pressures of natural nitrogen 
and samples containing 34.6 percent nitrogen-15 were 
measured by a differential method and were found to be 
given by log P:/P2=(0.2474/T) —0.001994 where 1 refers 
to the natural nitrogen and 2 refers to the heavy sample. 
The difference in the heats of vaporization of the two 
samples was calculated to be 1.14 calories, that of the heavy 
sample being higher. The triple point pressures were 
measured and found to be 9.386 cm of Hg and 9.378 cm 
‘of Hg, respectively. From this the difference in the triple 
points was calculated to be 0.020°K. Assuming that 
Raoult’s law may be used, the ratio of the vapor pressures 
and the difference in the heats of vaporization of the two 
pure isotopes, nitrogen-14 and nitrogen-15, were found 
to be 


log (Pis/P 15) = (0.7230/T) — 0.005822 
AH yap 15— AH yap 14 = 3.33 calories. 


The triple point of the pure nitrogen-15 was found to be 
0.058° higher than that of nitrogen-14 and the boiling 


point 0.052° higher. The differences in the vapor pressures 
of solid natural ammonia and of solid trideuteroammonia 
containing 98 percent deuterium were measured and found 
to be given by 


log (PNH;/PNDs) = (49.69/T) —0.1305. 


The vapor pressures of the liquid ammonias were found to 
follow the equation 


log (PNH3/PNDs) = (46.25/T) —0.14003. 


These two equations agree at the triple point. The differ- 
ence in the heats of sublimation was calculated to be 227 
calories, while the difference in the heats of vaporization is 
212 calories, the values for the trideuteroammonia being 
higher than those for ammonia. The triple point pressures 
were found to be 4.557 cm of Hg for ammonia and 4.822 
cm of Hg for the trideuteroammonia. This corresponds to 
triple point temperatures of 195.68°K and 198.79°K, re- 
spectively. The boiling point of the trideuteroammonia 
was found to be 2.37° higher than that of the ammonia. 


INTRODUCTION 


HE separation and concentration of iso- 
topes in recent years has made possible the 
study of the physical properties of isotopic com- 
pounds. This study is of great importance and 
interest for both theoretical and practical pur- 
poses. It has been found, for example, that there 
exist small differences in the vapor pressures of 
H,0'8 and H,O'*,!~ in the vapor pressures of 
and in the vapor pressures of 
and N“H;,° and in many other isotopes and their 
compounds. Much larger differences have been 


t Dissertation submitted by Isidor Kirshenbaum in 
partial fulfillment of the requirements for the degree of 
Doctor of Philosophy in the Faculty of Pure Science, 
Columbia University. Publication assisted by the Ernest 
Kempton Adams Fund for Physical Research of Columbia 
University. 

-1G. N. Lewis and R. E. Cornish, J. Am. Chem. Soc. 55, 
2616 (1933). 

2M. H. Wahl and H. C. Urey, J. Chem. Phys. 3, 411 
(1935). 

3D. F. Stedman, Can. J. Research 13, 114 (1935). 

*E. H. Riesenfeld and T. L. Tschang, Zeits. f. physik. 
Chemie B33, 127. (1936). 

5H. G. Thode, J. Am. Chem. Soc. 62, 581 (1940). 


observed in the deuterium compounds.®*’ These 
differences cannot be explained classically and 
must be due to second-order quantum mechan- 
ical effects. 

One might, on first thought, expect that such 
a simple principle as Trouton’s rule would be 
obeyed by two isotopic substances with greater 
exactness than by most compounds. Examination 
of the literature,* nevertheless, shows this not to 
be the case. There exist considerable deviations 
from the rule when isotopic compounds are com- 
pared. Upon a more detailed consideration, how- 
ever, these deviations are understandable in that 
Trouton’s rule (or its Hildebrand modification’) 
can be derived from classical statistical me- 
chanics'® whereas the differences in the physical 


6 E. H. Riesenfeld and T. L. Tschang, Zeits. f. physik. 
Chemie B33, 122 (1936). 

7Champetier, Bonnet, and Magat, ‘Deuterium et 
Composés de Deuterium,’’ Tables Annuelles de Constantes et 
Données Numerique (1937), p. 57. 

8J. R. Bates, J. O. Halford, and L. G. Anderson, J. 
Chem. Phys. 3, 531 (1935). 

9 J. H. Hildebrand, J. Am. Chem. Soc. 37, 970 (1915). 

1K. S. Pitzer, J. Chem. Phys. 7, 583 (1939); see also 
oo in J. H. Hildebrand, J. Chem. Phys. 7, 233 
1939). 
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properties of isotopes must be quantum me- 
chanical in origin." It will be shown later that a 
consideration of the quantum mechanical effects 
in isotopes leads to a partial understanding of 
the differences in the vapor pressures of the ni- 
trogen isotopes. It is, therefore, possible that a 
study and correlation of physical propefties of 
isotopes may lead to a better understanding of 
the liquid state. 

It was with this in mind that the differences in 
the vapor pressures and other thermodynamic 
properties of the liquid nitrogen isotopes and of 
solid and liquid ammonia and trideuteroammonia 
were measured. Nitrogen may be considered a 
pseudo-monatomic non-polar liquid, and conse- 
quently lends itself to a theoretical treatment 
more readily than do many other available iso- 
topic substances. High concentrations of nitro- 
gen-15 in the form of NH; have, moreover, 
recently been made available.’* With the results 
of this research, it is interesting to note, there is 
now available a complete set of data on the iso- 
topic effect on vapor pressures and other thermo- 
dynamic properties; namely, the data on deu- 
terium,'* ammonia-N”,> ammonia-d; and nitro- 
gen (15). The differences in the vapor pressures of 
liquid ammonia and ammonia-d3; have been re- 
ported previously by Jungers and Taylor." Their 
measurements, at only five temperatures, gave 
each of the vapor pressures to the nearest milli- 
meter. It was, therefore, thought advisable to 
repeat the measurements on ammonia-d;, at 
many more temperatures, measuring the differ- 
ences directly, in order to obtain them more 
accurately. 

Isotopes are, moreover, of great value in chem- 
ical, physical, and biological research. The work 
of Schoenheimer, Bonhoeffer, Urey, Rittenberg, 
H. S. Taylor, and others has amply demonstrated 
the value of and the need for isotopes in such 
researches. The separation and concentration of 
isotopes has, therefore, been a major problem for 
many years.’ Probably the most successful 


1 K.F. Herzfeld and E. Teller, Phys. Rev. 54, 912 (1938). 
(1939); G. Thode and H. C. Urey, J. Chem. Phys. 7, 34 
RB. Scott, F. G. Brickwedde, H. C. Urey, and M. H. 
Wahl, J. Chem. Phys. 2, 454 (1934). 
(19355 Jungers and H. S. Taylor, J. Chem. Phys. 2, 373 
16 See F. W. Aston, Mass Spectra and Isotopes (Long- 
mans, Green, New York, 1933). 
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method as yet developed for the large scale con- 
centration of light isotopes other than hydrogen 
is the one depending upon chemical exchange 
reactions.'* This process, which takes advantage 
of the slight differences in che chemical properties 
of isotopes,'’ has been used to produce relatively 
high concentrations of the rare isotopes of car- 
bon,'® nitrogen,” sulfur,'® etc.2° This method of 
separation produces these isotopes at a lower 
cost than do other methods, but is, nevertheless, 
expensive. A more economical process would be 
welcome. 

Much work has been done in recent years on 
close fractionation columns and many highly 
efficient columns have been developed. If, there- 
fore, two isotopes or their compounds have an 
appreciable difference in vapor pressure, they 
may be separated efficiently and economically 
by fractional distillation. As far back as 1931, 
Keesom and co-workers partially separated neon 
isotopes by fractional distillation of liquid neon. 
Deuterium was concentrated by this method the 
following year.” Huffman and Urey” have suc- 
ceeded in increasing the concentration of H,O"* 
by a factor of 4.5 using such a process, although 
the fractionation factor is only 1.003-1.006.*4 
The fractionation of the oxygen isotopes by dis- 
tillation of the liquid was observed by Klar and 
Kraus and by Smith.*%* Stedman* has concen- 
trated oxygen (18) using this process. He also suc- 
ceeded in partially separating the chlorine iso- 
topes in the form of carbon-tetrachloride. It has 
also been observed that the nitrogen isotopes 
are concentrated to a slight extent by distillation 


16 H. C. Urey, J. Wash. Acad. Sci. 30, 277 (1940) ; J. App. 
Phys. 12, 270 (1941). 

(1985). C. Urey and L. J. Greiff, J. Am. Chem. Soc. 57, 321 

18 C, A. Hutchinson, D. W. Stewart, and H. C. Urey, J. 
Chem. Phys. 8, 532 (1940). 

( 240). W. Stewart and K. Cohen, J. Chem. Phys. 8, 904 
1 

2 G. N. Lewis and R. T. MacDonald, J. Am. Chem. Soc. 
58, 2519 (1936); T. I. Taylor and H. C. Urey, J. Chem. 
Phys. 6, 429 (1938). 

21W. H. Keesom and H. van Dijk, Proc. Acad. Sci. 
Amst. 34, 42 (1931); W. H. Keesom, H. van Dijk, and J. 
Haantjes, Proc. Acad. Sci. Amst. 37, 615 (4934) ; Physica 1, 
1109 (1934). 

2H. C. Urey, F. G. Brickwedde, and G. M. Murphy, 
Phys. Rev. 40, 1 (1932). 

23]. R. Huffman and H. C. Urey, Ind. Eng. Chem. 29, 
331 (1937). 

24 R, Klar and A. Kraus, Naturwiss. 22, 119 (1934); E. R. 
Smith, J. Chem. Phys. 2, 298 (1934). 
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of liquid ammonia.” The fractionation factor is 
in this case only 1.002.5 If, therefore, the nitrogen 
isotopes have a large enough difference in vapor 
pressures, this method may be used. It is one 
object of this paper to supply the fractionation 
factors for the nitrogen and ammonia isotopes, 
i.e., the ratios of the vapor pressures. These 
factors are absolutely necessary for determining 
the feasibility of separating these isotopes by 
distillation. 


APPARATUS 


In order to measure with precision the small 
differences between the vapor pressures of the 
nitrogen isotopes, a differential method must be 
used. The apparatus used in this investigation is 
similar to the one developed by Thode’® in his 
work on ammonia-N™, but with several essential 
modifications which were necessitated by the 
much lower temperatures of liquid nitrogen. 

The apparatus is shown diagrammatically in 
Fig. 1. A small glass cell of Jena 16 III A was 
divided into two vacuum tight compartments by 
a platinum partition. The compartments con- 
tained small platinum beads. The platinum 


26 M. H. Wahl, J. R. Huffman, and J. Hipple, J. Chem. 
Phys. 3, 434 (1935); H. C. Urey and A. H. W. Aten, Jr., 
Phys. Rev. 50, 575 (1936). 


partition and beads were necessary for the rapid 
equalization of temperature between the two 
isotopic liquids whose vapor pressures were being 
compared. For a 1 percent or better precision in 
the differences between the vapor pressures of 
the nitrogen isotopes, a temperature difference 
of +0.0002°C or less had to be maintained 
between the two compartments. 

The cell was suspended by 6-mm Pyrex leads 
L, L’ in a glass jacket which could be evacuated 
through B to 10~* mm or better. The leads L, 
L’, and B were then passed through a silvered 
Dewar flask S and thence to the measuring ap- 
paratus. Liquid nitrogen was used in the Dewar 
during the nitrogen runs to precool the leads. 
Dry ice and chloroform were used in the am- 
monia runs. The leads Z and L’ were joined by 
glass rod at several points. This resulted in a 
sturdier apparatus and also aided in the tem- 
perature control. The whole apparatus from O 
down was immersed in a cryostat consisting of a 
large Dewar flask F and a steel case C. Both the 
Dewar and case had windows through which the 
cell, and samples could be observed. The windows 
were covered with aluminum foil during measure- 
ments in order to minimize radiation effects. The 
temperature of the cryostat was controlled by 
regulating the pressure under which the cryostat 
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liquids, nitrogen or ammonia, boiled by means 
of a barostat and Hy-vac pumps at V. The 
nitrogen bath was stirred by nitrogen gas using 
a lift pump device P. It was found advisable to 
modify the apparatus for the ammonia runs. The 
cryostat was then connected to the barostat by 
means of a large modified Dewar dry ice trap. 
This served as a condenser for the ammonia, 
which was stirred by a small nichrome heater 
using about 16-20 watts. 

The leads L, L’ were connected to the gas 
reservoir bulbs and to the arms of a differential 
manometer D which contained Apiezon B oil for 
the nitrogen runs and mercury for the ammonia 
runs. The Apiezon oil was well outgassed before 
use and was found to have a density of 0.8582 at 
27.55°C. The manometer had a vacuum-seated, 
mercury cup stopcock separating the arms. This 
type of stopcock was used because it was found 
that the Apiezon oil would dissolve the grease 
on the plug, and water would leak through an 
ordinary, well-ground stopcock in a relatively 
short time. The manometer was immersed in a 
water bath maintained at 27.55°+0.02C. 

The side containing the natural samples was 
also connected to a mercury manometer M, 
which was used to measure the absolute pres- 
sures. The absolute temperatures of the samples 
were then calculated from the observed vapor 
pressure curves of ordinary nitrogen®*® and am- 
monia.?” The differential manometer was ob- 
served though a small cathetometer precise to 
+0.02 mm. The measurements of the absolute 
pressure were made with a Gaertner M-900 
Cathetometer precise to +0.05 mm with a scale 
accuracy of 0.003 mm. The small cathetometer 
was calibrated using the scale of the larger one 
as standard. 

The leads from the cryostat to the manometers 
were lagged with asbestos cord because it was 
found that a sudden change of temperature at 
any point in one of the leads would upset the 
liquid-gas equilibrium in the corresponding cell 
compartment. A relatively long period of time, 
as compared to the time of reading, was then 
needed to re-establish thermal equilibrium. All 
stopcocks in the apparatus were precision ground 


*® F, Henning and J. Otto, Physik. Zeits. 37, 634 (1936). 
*7 Roy Overstreet and W. F. Giauque, J. Am. Chem. Soc. 
59, 254 (1937). 
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and vacuum seated and were greased with 
Apiezon M. 


PREPARATION OF GASES 


The nitrogen samples were prepared by circu- 
lating ammonia over hot copper oxide. The 
ammonia was obtained from C.P. ammonium 
chloride using C.P. sodium hydroxide. Concen- 
trated solutions of both were boiled under 
vacuum and then mixed. The ammonia was 
dried by passing it through a column containing 
sodium hydroxide pellets; it was collected in a 
liquid nitrogen trap, evacuated, and then led into 
the oxidizing system which had previously been 
outgassed to better than 10-*-mm pressure. 

After the ammonia was oxidized, the water and 
any oxides of nitrogen that may have formed 
were frozen out in a liquid nitrogen trap, and 
the nitrogen gas was pumped into a storage bulb 
which was then sealed on to the measuring 
apparatus. Two samples of ordinary and two 
samples of heavy nitrogen were made on different 
occasions. 

The samples of nitrogen were analyzed for 
purity and isotopic content on a Bleakney type 
mass spectrometer. The isotopic content was 
checked on a Nier mass spectrometer. The two 
analyses agreed within experimental error. No 
precautions were taken to remove the traces of 
methylamine in the samples of heavy ammonia 
as reported by Thode’ since the methylamine was 
completely oxidized by the hot copper oxide. No 
impurities were found in the samples. The purity 
of the samples was further checked by blank 
determinations (v.i.) and measurement of the 
triple point pressure. 

The ammonia for the ammonia measurements 
was prepared by reacting water with magnesium 
nitride. The nitride was made in a steel com- 
bustion tube”* from reagent pure magnesium and 
tank nitrogen. The nitrogen, which was 99.8 per- 
cent pure, was further purified by passing it 
first over hot copper wire and then through a 
dry ice-trichlorethylene trap. The steel tube con- 
taining the magnesium was well evacuated and 
then heated to 700-800°C. The nitrogen was 


28 The Mg;N2 had been first prepared in a Corning 702 
combustion tube. The tube blackened considerably due to 
the formation of Si and Mg.Si. The product resulting from 
the reaction with water contained a considerable amount of 
a foul smelling gas believed to be silane or its homologs. 
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TABLE I. A typical set of readings. 


Time after 


evacuating jacket P:i—P:cm of oil Picm of Hg 


1.439 
1.434 
1.438 
1.441 
1.443 
1.438 
1.440 
1.443 
1.444 
1.443 


40 minutes 
41 25.025 


25.055 
25.090 
25.120 
25.150 
25.175 
25.200 
25.245 
25.265 


passed in at about one atmosphere. Precautions 
had to be taken to prevent the loss of too much 
magnesium by distillation to cold parts of the 
system. 

The nitride was transferred to another system 
and was outgassed at 400°C for several hours. 
The system was evacuated to 10-* mm and water 
was added to the nitride in vacuum. The am- 
monia formed was fractionally distilled several 
times in vacuum from dry ice temperatures to 
liquid nitrogen temperatures. The samples were 
then evacuated at —95°C to remove any silanes 
that might have formed during the outgassing 
process. Each fraction was again fractionated 
and led into the storage bulb. Only about one- 
half of the ammonia formed was used for 
measurements. The water used to make the 
deuteroammonia contained 98 percent D,O. Two 
samples of ammonia and one of ammonia-d3 were 
prepared. The purity of the samples was checked 
by blank determinations and measurement of 
triple point pressures.”® 


ISOTOPIC CONCENTRATION AND EQUILIBRIUM 


The equilibrium constant for N2“+N,%— 
2N“N”® at 300°A was calculated using the 
method of Urey and Greiff'’ to be 3.993. The 
spectroscopic data was obtained from Herzberg.*° 


2° The Mg(OH)2 remaining from the reaction of MgsNz 
= was analyzed qualitatively for silicon. None was 
ound. 

% G. Herzberg, Molecular Spectra and Molecular Structure 
(Prentice-Hall, Inc., 1939), vol. 1, p. 490. 
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The samples of heavy nitrogen were then ana- 
lyzed on a Nier mass spectrometer for isotopic 
equilibrium. The observed equilibrium constant 
was 3.78. This non-equilibrium condition is ex- 
plicable in light of the observation by Thode’ 
that the heavy nitrogen ammonia contained 
methylamine, which had a natural nitrogen con- 
centration and the observation by Joris and 
Taylor* that hot copper will not establish iso- 
topic equilibrium between nitrogen isotopes. 
A sample was then heated to 1650-—+1800°C for 
a long time by a tungsten filament. The analysis 
then showed the equilibrium constant to be 
3.98+0.03. 

Measurements of the vapor pressure were 
made on the original samples. If one assumes 
that the vapor pressure of N“\N® is the arithmetic 
mean of the vapor pressures of N,“ and N,", 
then using a non-equilibrium sample introduces 
no error. Since the ratio of the vapor pressures 
of the pure isotopes is about 1.01, the geometric 
mean and arithmetic mean differ by only 0.00001. 
Assuming, therefore, that the vapor pressure of 
N*“N® is the geometric mean of the other vapor 
pressures, an error of less than 0.1 percent is 
introduced by using the non-equilibrium samples 
for measurements. Moreover, since the geo- 
metric mean and the arithmetic mean of the two 
vapor pressures are the same well within experi- 
mental error, Raoult’s law may be applied to the 
experimental results without concern as to 
whether the law applies merely to mixtures of 
N," and N,", or to equilibrium samples as well. 

Analysis on a mass spectrometer showed that 
the heavy nitrogen samples contained 34.6 per- 
cent+0.2N”. The ordinary nitrogen, having been 
prepared from ordinary ammonium chloride, had 
an isotopic concentration of 0.38 percent N".' 


TABLE II. Same sample of nitrogen (14) in 
each compartment. 


AP cm oil 


—0.010 
+0.008 
—0.013 
— 0.004 
+0.002 
— 0.008 
—0.014 


P cm Hg 


11.0 
13.7 
16.1 
22.8 
23.5 
35.0 
37.6 


AP cm Hg 


— 0.00063 
+0.00050 
— 0.00082 
— 0.00025 
+0.00013 
— 0.00050 
— 0.00080 


31 G. G. Joris and H. S. Taylor, J. Chem. Phys. 7, 893 
(1939). 
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Samples of the ammonia-ds were decomposed by 
a hot tungsten filament, and the resulting gases 
were analyzed on a Bleakney type mass spec- 
trometer for their deuterium-hydrogen ratios. 
These indicated that the ammonia-d; contained 
97.9 percent deuterium. The ordinary ammonia, 
having been prepared from distilled water, con- 
tained 1 part deuterium per 6500 parts of 
hydrogen.® 


THE VAPOR PRESSURE DIFFERENCES OF THE 
LIQUID NITROGEN ISOTOPES 


A. Method of Measurement 


About half a gram of nitrogen gas was con- 
densed from the reservoir bulbs into each of the 
compartments of the cell A by means of solid 
nitrogen. The solid nitrogen was formed by 
pumping on liquid nitrogen in the cryostat. 
Precautions were taken to prevent condensation 
of nitrogen in the leads. The levels of the liquefied 
nitrogen samples in each compartment, after 
condensation was complete, were a trifle above 
the top of the platinum beads. The jacket J was 
filled with dry hydrogen gas and the Dewar S 
was filled with liquid nitrogen during condensa- 
tion. The Dewar was kept filled during the 
measurements. 

The compartments were then connected to the 
oil manometer and the jacket was evacuated. 
The cryostat was filled with liquid nitrogen and 
regulated to the desired temperature. Hydrogen 
was again added to the jacket and the samples 
in the cell were permitted to come to the tem- 
perature of the cryostat. When the correct tem- 
perature was attained, the jacket was evacuated 
to 10-* mm or better and after forty minutes 
readings were made with the cathetometers. 


TABLE III. Different samples of nitrogen (14) in 
each compartment. 


AP cm oil 


+0.006 
—0.013 
— 0.008 
— 0.004 
—0.012 
—0.008 
—0.023 
—0.027 


cm Hg 


+0.00038 
— 0.00082 
— 0.00050 
—0.00025 
—0.00076 
— 0.00050 
—0.00145 
—0.00170 


* N.F. Halland T. O. Jones, J. Am. Chem. Soc. 58, 1915 
(1936). 


TABLE IV. 


Pi—P: cm oil 


0.744 
0.766 
0.745 


1.202 


* Observed in same run. 


Readings were taken for twenty or thirty min- 
utes, with a reading of the difference in vapor 
pressures being taken every second minute. The 
absolute vapor pressure of the natural nitrogen 
was measured before and after each reading of 
the vapor pressure difference. Thus ten to fifteen 
observations of the difference were made. During 
a set of measurements the temperature of the 
cell never changed by more than 0.08°C, usually 
by 0.04°C. The average deviation from the re- 
ported temperature is, of course, much smaller— 
of the order of +0.01°C. The two compartments, 
however, never differed by more than +0.0002°C 
in temperature, as may be calculated from the 
blank determinations. The cryostat was then 
regulated at another temperature and the above 
procedure was repeated at the new temperature. 
Liquid nitrogen was added to the cryostat when- 
ever one quarter of the bath had evaporated. 
A typical set of readings is shown in Table I. 
The temperature during the set of readings 
varied 0.06°C. The average deviation from the 
average temperature, used as a final value, is 
much less than +0.01°C. The following symbols 
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will be used in discussing the results of the 
nitrogen measurements. P, is the pressure of the 
sample of natural nitrogen in cm of mercury; 
Pz, is the pressure in cm of mercury of the heavy 
sample containing 34.6 percent N", Py, and Pis 
refer to the vapor pressures of the pure isotopes, 
and the symbol N, refers to the gas. 

The accuracy of the apparatus was tested by 
running a blank with one sample of natural 
nitrogen. Both the precision of the apparatus 
and the purity of the samples were ascertained by 
measuring the differences in the vapor pressures 
from 64°K to 75°K of two different samples of 
natural nitrogen. The measurements on the 
heavy nitrogen were also made on two different 
samples to assure purity and reproducibility. 


B. Results 


Tables II and III show the results of the 
blanks. Although the absolute error in the meas- 
urement of AP increased with temperature, the 
percent error did not. It is also evident that in 
no case did the differences in the vapor pressures 
of the two natural samples of nitrogen (14) exceed 
the experimental error using one sample. 

Table IV gives the results obtained with the 
heavy samples. It includes the vapor pressure of 
the natural nitrogen P;, as determined on the 
mercury manometer, the vapor pressure differ- 
ence P,;—P, in cm of oil as measured on the 
differential oil manometer, where P2 is the vapor 
pressure of the nitrogen containing 34.6 percent 
N". The table also indicates on which sample of 
heavy nitrogen the measurement was made I, II, 
or III, where III is a mixture of I and II. Every 
P,—/2 is the average of at least ten readings. 

By a series of approximations, introducing 
errors negligible compared to the experimental 
errors, the values for P,;—P: in set A were re- 
duced to give values for P;—P» at a pressure of 


TABLE V. Vapor pressures of the nitrogens. 


(Pi1—P2) 
obs.— 
Picm Pi-P2 A.D. (Pi-P2) Piu/Pis 
Hg TK cmHg % Pi/P2 calc. calculated 
11.250 64.11 0.0480 0.8 1.00428 —0.0002 1.01264 


11.900 67.50 0.0771 0.9 1.00389 +0.0004 1.01132 
25.100 69.00 . 0.0921 0.5 1.00368 +0.0003 1.01078 
31.110 70.46 0.1081 0.5 1.00349 —0.0003 1.01027 
42.022 72.61 0.1368 0.6 1.00327 +0.0008 1.00957 
54.500 74.61 0.1646 0.6 1.00303 —0.0005 1.00894 
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Fic. 2. Vapor pressure of the nitrogens. log P:/P2 vs. 1/T. 


11.250 cm, those in set B were calculated to 
give values for P;—P: at a pressure of 19.900, 
etc. These values were then converted to centi- 
meters of mercury, averaged, and are reported 
in Table V. The temperatures were calculated 
from the observed absolute vapor pressure of 
nitrogen (14) using?® 


log = 7.781745 
—0.006249T — (341.619/T) (1) 


and are, as has been indicated, good to +0.01°. 
The table includes the A.D. in P,—P2. It also 
includes the ratio Pi/P, and the differences 
between the observed (P:—P2)’s and those calcu- 
lated by Eq. (3). 

Assuming the differences in the vapor pressures 
of the two nitrogens follow an equation of the 
type 


logio Pi/P2=(A/T)+B (2) 


and applying the method of least squares to the 
data in Table V, it is found that the vapor 
pressure data for the liquids can be repre- 
sented by 


log P:/P2=(0.2474/T) —0.001994. (3) 


The P;/Pz2 is obtained from the observed differ- 
ence P,—P»2. The accuracy with which this equa- 
tion fits the data can be seen from Fig. 2 and 
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from the results in next to the last column of 
Table V. The mean difference between calcu- 
lated and observed values of the vapor pressure 
differences (P;—P2) is +0.4 percent. The experi- 
mental error was estimated to be about 1.1 
percent. 

It has been shown® that H.—De mixtures 
follow Raoult’s law. It is, therefore, reasonably 
safe to assume that nitrogen isotopes do the 
same within our experimental error of 1 percent. 
Making this assumption that the vapor pressure 
of this isotopic mixture is a linear function of the 
mole fraction, it is found that 


log Pis/Pis= (0.7230/T) —0.005822 (4) 


expresses the relationship between the vapor 
pressures of pure liquid nitrogen (14) and pure 


TABLE VI. 


_ Sample A—ordinary Nitrogen—14 
Melting curve reezing curve 


9.384 cm Hg 9.384 cm Hg 
9.388 9.384 

9.388 9.390 
Average: 9.386+0.002 cm Hg 


Sample B—Nitrogen—34.6% N'5 
Melting curve 


Cooling curve 
9.386 cm Hg 9.377 cm Hg 9.372 cm Hg 
9.374 9.377 9.382 
9.379 9.377 


9.374 9.377 
Average: 9.378+0.003 cm Hg 


liquid nitrogen (15). The values of Pi4/P15 calcu- 
lated at several temperatures by means of this 


_ equation are given in the last column of Table V. 


C. Differences in the Heats of Vaporization 
and Boiling Points 


Using Eqs. (1) and (4) and extrapolating to 
P1;=76 cm, it is found that nitrogen (15) boils 
0.052° (77.404°K) higher than does nitrogen (14) 
(77.352°K).26 It is interesting to note that the 
difference in the normal boiling points (A7;) of 
ammonia-N“ and ammonia-N"™® is 0.055°, that 
of ammonia-N" being higher. If the vapor 
pressures of the isotopes follow an equation of 
the form 


logio (AH, —AH,/2.303RT)+C, (5) 


where the subscripts have the same significance 


8 G. N. Lewis and W. T. Hanson, J. Am. Chem. Soc. 56, 
1000 (1939). 


as before, then the difference in the heats of 
vaporization of the two isotopic samples of 
nitrogen can be calculated from Eqs. (3) and (5) 
to be 

AH,—AH,= 1.14 cal., 


where AH, was found by Giauque and Clayton* 
to be 1332.9 at the boiling point. Again assuming 
an ideal solution, a simple extrapolation shows 
the heat of vaporization of pure nitrogen (15) to 
be 3.33 calories higher than the heat of vaporiza- 
tion of pure nitrogen (14). 


DETERMINATION OF THE TRIPLE 
POINT OF NITROGEN 


The vapor pressures at the triple points were 
obtained for each of the two isotopic samples 
from both melting and freezing curves, obtained 
by plotting vapor pressure at and near the triple 
points against time. A sharp break in these 
curves was observed at the triple points. This 
break could be reproduced with great precision. 
The results obtained are given in Table VI. The 
values are corrected to 0°C. The most recent 
values for the triple point pressure of ordinary 
nitrogen are 9.401, 9.401,** 9.360,3* and 9.39157 
cm of mercury, omitting the 9.46 reported by 
Henning and Otto.*® The average of these values 


TABLE VII. The vapor pressures of the solid ammonias. 


(PNH; 
—Pyp, obs.- 
PNH; No. of (PNH; 
cm of observa- Pyy,/PNpD; —PND; 
Hg tions calc. calc. cm tte 


2.57 189.7 
2.80 190.5 
2.98 191.1 
191.4 
191.9 
192.1 
192.4 
192.8 
» 193.2 
193.5 
194.2 
194.7 
195.0 
195.1 
195.3 


+0.00 
+0.02 
+0.01 
+0.00 
+0.01 
+0.01 
—0.01 
—0.02 
+0.02 
+0.00 
—0.02 
+0.00 
+0.01 
+0.00 
+0.03 


— 

NWN WR Whee 


34 W. F. Giauque and J. O. Clayton, J. Am. Chem. Soc. 
55, 4875 (1933). 
36 W. H. Keesom and A. Bijl, Physica 4, 305 (1937). 


( ped » T. H. Verschoyle, Proc. Roy. Soc. London A230, 189 
1931). 


87 E. Justi, Ann. d. Physik [5] 10, 983 (1931). 
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Fic. 3. Vapor pressures of the liquid ammonias. 
log PNH3/PNDs vs. 1/T. 


is 9.388 cm of mercury. The value given in this 
paper, 9.386 cm of Hg, is in excellent agreement 
with this average. 

The difference in the triple points of the two 
samples can be calculated from the above data. 
With the two samples at the triple point of 
ordinary nitrogen (63.145°K), the vapor pres- 
sure of the ordinary sample would be 0.0415 cm 
higher, as calculated from Eq. (3). It was found, 
however, that at their respective triple points the 
light sample has a vapor pressure that is only 
0.008 cm higher. The heavy sample must, there- 
fore, melt at a higher temperature, and this in- 
crease in temperature corresponds to the increase 
in pressure of (0.0415—0.008) cm of Hg. Since 
from Eq. (1), dP/dT=1.7 cm/deg. at the triple 
point, the difference in the triple points of ordi- 
nary and nitrogen-34.6 percent N" is 


(0.0415 — 0.008) /1.7 =0.020°K. 


Assuming that the differences in the triple points 
for isotopic mixtures of nitrogen are proportional 
to the mole fractions, nitrogen (15) is found to 
melt 0.058° higher than does nitrogen (14). Since 
the triple point pressures were reproducible to 
+0.002 and +0.003 cm, respectively, an error 
of 0.0036 cm may exist in the pressure difference. 
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This corresponds to a 10 percent error in the 
difference between the triple point temperatures. 


THE VAPOR PRESSURE DIFFERENCES OF THE 
AMMONIA ISOTOPES 


A. Solid 


About half a gram of ammonia gas was con- 
densed into each of the compartments by means 
of liquid nitrogen. Precautions similar to those 
in the nitrogen runs had to be taken. The Dewar 
S was filled with dry ice and chloroform. The 
jacket, which contained dry hydrogen gas, was 
then evacuated and the cryostat containing 
liquid ammonia was screwed into place. After 
evacuating the jacket for twenty minutes, read- 
ings were taken. Since the temperature changed 
rather rapidly, no more than one observation 
could be made, at any one temperature at a 
given time. Many of the points, however, were 
observed more than once. Blank determinations 
indicated that the differences may be too low by 
at most 0.04 cm of mercury. This was due to 
the very low thermal conductivity of the non- 
metallic solid. Table VII summarizes the results. 
All pressures are in centimeters of mercury and 
have been corrected to 0°C. The temperatures 


TABLE VIII. The vapor pressures of the liquid ammonias. 


(PNH;—PND, 
PN Hs 
ND; 


1.41 
1.56 
1.93 
2.15 


PNH; 
cm Hg 


obs. 
(PNHi~ 
average calc. 


ND;) 


1.227 
1.225 
1.210 
1.206 
1.203 


1.197 
1.195 


1.189 


1.185 
1.175 


1.172 
1.168 


—0.01 
+0.01 
—0.01 
—0.01 
+0.00 


—0.02 
+0.00 


—0.01 


—0.02 
—0.02 


+0.00 
—0.01 


7.63 
8.51 
11.11 
12.59 
13.48 


15.16 
16.20 


18.34 


20.25 
25.40 


27.50 


30.20 
41.50 —0.01 


+0.02 
+0.06 


52.93 
62.16 
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were calculated from the observed vapor pres- 
sures of solid ammonia using?’ 


log PNH3= — (1630.700/T) +9.00593. (6) 


Assuming that Eq. (2) may be used, and applying 
the theory of least squares to the data in Table 
VII, it is found that 


log (Pxus/PND3) = (49.69/T) —0.1305. (7) 


The accuracy with which this equation describes 
the experimental data is also indicated in the 
table. The equation for the vapor pressure of 
solid trideuteroammonia may be obtained by 
subtracting (7) from (6). 


log PND3= — (1680.4/7) +9.1364. (8) 


This equation agrees exactly with that found for 
the liquid ND; (Eq. (11)) at the triple point. 
Since altogether different data were used in calcu- 
lating each of these equations, this agreement is 
a check on the consistency of the measurements 
and justifies extrapolating the equations to the 
triple point. Equation (7) will give PNu;— PND; 
accurately to about +0.05 cm and may be used 
to calculate T to 0.12°K or better. Assuming 
that the vapor pressures of the solids follow 
Eq. (5), the difference in the heats of sublimation 
may be calculated to be 6AH,.,»=227 calories, 
that of ammonia-d; being higher, 6 indicating 
differences between isotopes. 


B. Liquid 


After determining the difference in the vapor 
pressures of the solids, the measurements were 
continued into the liquid. Here thermal equi- 
librium was much more easily attained and the 
error in the measurements is consequently much 
lower. The method of measurement was similar 
to that used in the nitrogen experiments and 
similar precautions were taken. The accuracy of 
the apparatus was tested by running blank de- 


TABLE IX. 


Triple point pressure 
of ammon' 


Melting 
curve 


4.562 cm of Hg 
4.557 


Triple point pressure 
of ammonia-d3 


Freezing 
curve 


Melting 
curve 


Freezing 
curve 


4.564 cm of Hg 
4.550 
4.554 


4.823 cm of Hg 4.816 
4.824 


4.822 


4.823 
Average: 4.557 cm of Hg Average: 4.822 cm of Hg 
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Fic. 4. Vapor pressures of liquid ammonias. 
T log PNH3/PND; vs. 1/T. 


terminations with one sample of ammonia at 
fifteen points between 196°K and 236°K. Each 
point was observed ten times. The average error 
was +0.015 cm of Hg with practically all of the 
errors being negative, i.e., PNH;— PND; observa- 
tions are too low by about 0.02 cm of Hg. The 
maximum error observed was —0.028 cm of Hg. 
This was observed only once and under the most 
adverse experimental conditions and is, there- 
fore, the greatest possible error. Precautions were 
taken to avoid such conditions during the meas- 
urements. The precision of the apparatus and the 
purity of the samples were checked by measuring 
the differences in the vapor pressures from 201°K 
to 236°K of the two different samples of ordinary 
ammonia. The error was found to be negative, 
i.e., —0.007 cm average. 

Table VIII summarizes the results obtained 
with liquid ammonia-d;. Each PNu3— PND; is the 
average of at least ten readings. The table also 
includes the differences between the observed 
vapor pressure differences and those calculated 
from Eq. (10). All pressures are again given in 
centimeters of mercury and have been corrected 
to 0°C. The temperatures were calculated from 
the observed absolute vapor pressures of am- 
monia using?’ 


log Pnus= — (1612.500/T) —0.0123117 
+0.0000125217?+10.83997 (9) 
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Fic. 5. Triple point pressures of ammonia. 


and are correct to about +0.02°-+0.03°K. 
Again assuming that Eq. (2) may be used and 
applying the method of least squares to the data 
in Table VIII, it is found that the vapor pressure 
data for the liquids can be represented by 


log PNu3/PND3= (46.25/T) —0.14003. 
Subtracting Eq. (10) from (9) gives 


log Pnp3= — (1658.75/T) —0.012311T 
+0.00001252172+10.98000. (11) 


The accuracy with which Eq. (10) fits the data 
can be seen from Fig. 3 and from the results in 
the last column of Table VIII. An examination 
of Fig. 4 and the trend in deviations in Table 
VIII shows that log Pnu3/PND3; may not be a 
linear function of 1/T but rather of 1/7*. Since 
the 1/T curve seems to fit better, it was not 
thought advisable to find any correction terms. 
Using Eq. (11) and extrapolating to P=76 cm, 
it is found that the ammonia-d; boils at 242.05°K 
compared to a boiling point of 239.68°K re- 
ported?’ for ammonia. The difference in the boil- 
ing points is thus 2.37° and should be compared 
to a difference of 2.6° reported by Jungers and 
Taylor for ammonia-99 percent d;. If the vapor 
pressures of the liquid ammonias obey Eq. (5), 
then the difference in the heats of vaporization 
of the two ammonias is AHNp3;—AHNnu3=212 


(10) 


calories, where AHNuz is 5581 calories*’ at the 
normal boiling point. Jungers and Taylor re- 
ported a difference in heats of vaporization of 
193 calories. 

The vapor pressures at the triple points were 
obtained in a manner similar to that used in the 
nitrogen experiments. A greater tendency for 
supercooling was observed with the ammonias. 
The results obtained are given in Table IX. 
(See Fig. 5.) The triple point pressure of ammonia 
is in excellent agreement with the 4.558 cm re- 
ported by Giauque.?’ Substituting the observed 
triple point pressure in Giauque’s equation gives 
195.68°K as the triple point temperature of 
ammonia. The triple point pressure for ammonia- 
d; is smaller than the 4.86 cm reported by de 
Bruyne and Smyth.* Using Eq. (11) and the ob- 
served triple point pressure for ammonia-d;, the 


TABLE X. 


“NUNS 


NHs 
“NH3—d 


77 1.070 
110 1.006 1.065 
310 1.005 1.050 
550 1.004 1.042 
760 


H:O 
~H2O—d 


1.062 


1.025 


3% J. M. A. de Bruyne and C. P. Smyth, J. Am. Chem. 
Soc. 57, 1203 (1935). 
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triple point temperature of ammonia-98 per- 
cent d; is calculated to be 198.79°K. Using Eq. 
(8) for the solid ammonia-d;, the triple point is 
also calculated to be 198.79°K. Jungers and 
Taylor" reported the triple point temperature to 
be 199°K. De Bruyne and Smyth reported it as 
199.56°K. This, however, must be in error for if 
their observed triple point pressure is substituted 
into Eq. (11), one obtains 198.9°K as the triple 
point temperature. 


SEPARATION OF ISOTOPES 


In calculating the feasibility of separating 
isotopes by fractional distillation, the fractiona- 
tion factor a must be known. This is the ratio 
of the vapor pressures and can be calculated for 
the nitrogen and ammonia isotopes from the 
data in this paper. The assumptions will be 
made? that and that 


PNH3-d=(PNH3 The fractionation 
factors are listed in Table X, together with 
values for HOO—H,O-d and H.O — Com- 
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parison of the fractionation factors for the nitro- 
gen isotopes with those for HAO—H,O"'8? and 
consideration of the results obtained by Huffman 
and Urey** in separating the latter by distillation 
indicate that, on the same scale, the chemical 
exchange reaction™'® is a much more satis- 
factory method of concentrating the nitrogen 
isotopes. The chemical reaction 


NH,(g) +N“H,+ (sol) 
(sol) 


has a fractionation factor equal to 1.023. 

Once the fractionation factors are known, the 
decision as to whether a distillation method can 
be used depends upon the engineering problems 
involved and these in turn depend upon the scale 
of operations. The great difficulty in carrying out 
distillations at liquid nitrogen or liquid ammonia 
temperatures makes this method definitely un- 
satisfactory for any but large scale operations. 
And whether this method can then be used de- 
pends upon the solution of such problems as 
efficient heat transfer. 
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The differences in the thermodynamic properties of liquid nitrogen (14) and liquid nitrogen (15) 
are discussed on the basis of several Debye-like models. The quantum mechanical effects are 
found to be greater than one could expect from an harmonic oscillator model. Other factors 
affecting these properties are considered. 


ISOTOPE EFFECT IN THE GAS 


HE thermodynamic properties of the nitro- 

gen isotopes in the gas phase will be needed 

in the discussion that follows. These have been 
calculated assuming that nitrogen is a perfect 
gas. This assumption introduces smaller order 


t Dissertation submitted by Isidor Kirshenbaum in 
partial fulfillment of the requirements for the degree of 
Doctor of Philosophy in the Faculty of Pure Science, 
Columbia University. Publication assisted by the Ernest 
Kempton Adams Fund for Physical Research of Columbia 
University. 


effects than those to be discussed. The results of 
these calculations are shown in Table I. The 
constants for the heavier isotope were calculated 
from those observed for nitrogen (14). The nitro- 
gens are in their lowest vibrational state at the 
temperatures concerned and the rotational de- 
grees of freedom are practically classical. The 
equations used? are 


1G. Herzberg, Molecular Spectra and Molecular Structure 
(Prentice-Hall, Inc., 1939), Vol. 1, p. 490. 

2J. E. Mayer and M. G. Mayer, Statistical Mechanics 
(John Wiley and Sons, Ltd., 1940). 
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TABLE I. Isotope effect in the gas. 


Entropy of gas: 
N.!4 Sis-14 
Rotational S, 6.749 e.u. 6.886 e.u. 0.137 e.u. 
Translational S; 32.358 32.587 0.229 
Total S, 39.107 39.473 0.366 


Heat content of gas: 
Zero point energy 3362.11 3248.39  —113.72 cal. 
Rotational 123.612 123.739 0.127 cal. 


Strans = R(3/2 In M+5/2 In T 
—In Pytmos— 1.157), (1) 
Srot = R(1—In ¢—2/90— ++), (2) 


where 
o=Bhe/kT =h?/8x?IkT 


Herans=5/2RT, (3) 
Trot = RT (1 — 0/3 —0?/45). (4) 


The differences for the isotopes may be obtained 
directly by differentiation with respect to the 
mass. The following symbols will be used in this 
discussion : 6 indicates a difference between iso- 
topes; A indicates the difference between phases, 
e.g., AHyap is the heat of vaporization; Pi, and 
Ps; are the vapor pressures of the two isotopic 
compounds; and the subscript g will be used for 
gas, L for liquid, and D for Debye; e.g., Sp is the 
Debye entropy. For convenience, the experi- 
mental data to be explained (and that which is to 
be used in calculations) are outlined in Table II. 
The 6AS,.p was calculated from 6AHy,,/T. All 
calculations are at 63.203°K (triple point of 
nitrogen (15)), unless indicated otherwise. Sz_: is 
the entropy of the liquid due to the three transla- 
tional degrees of freedom. 


DISCUSSION 


The question as to why there is a difference in 
the vapor pressures of isotopic compounds di- 
vides itself into two parts. One is the problem 
of the liquid state and the other is the more 
general problem of the isotope effect. The latter 
raises no theoretical difficulties even in the most 
complicated molecules.? The methods of calcu- 


3 Teller’s product theorem ewe y by W. R. Angus 
et. al., J. Chem. Soc. p. 977 (1936); E | Bright Wilson, J. 
Chem. Phys. 9, 81 (1941). 


lating shifts in vibrational frequencies, moments 
of inertia, etc. are well known. As to the former, 
however, it is still beyond the power of statistical 
mechanics to deduce quantitatively and rigor- 
ously the properties of a liquid from a given 
molecular model. It is, therefore, necessary to use 
various approximate models when attempting to 
explain properties of liquids. The type of model 
used depends, of course, upon the properties that 
have to be explained. When dealing with liquids 
near the triple point, a very useful, although 
crude, approximation is one assuming a quasi- 
crystalline state. X-ray analysis on symmetrical 
molecules leads to the conclusion that a certain 
amount of regularity still exists in liquids near 
the triple point. Kratkey* has shown that the 
x-ray diagram of simple liquids can be explained 
by assuming any atom to be the center of a 
crystal lattice with a somewhat diffused struc- 
ture. The atoms thus vibrate in a cell of neigh- 
boring molecules in a manner similar to the 
vibrations in the solid. One can then, in principle, 
set up the problem in a convenient set of coordi- 
nates and then transfer to a set of normal 
coordinates. If © is the number of degrees of 
freedom per molecule, then 


ON. (5) 


Assuming an harmonic oscillator, then 
S= N(v)Sp(v)dv, (6) 
0 


where Sp(v) is the Planck expression for the 


TABLE II. Experimental data. 


log Pis/P1s = (0.7230/T) —0.005822 

6A Hvap = 3.33 calories (independent of temperature) 

6A Svap=0.053 e.u. 

= 1438 calories (value calculated from Giauque’s 
ata 

AHvapi;= 1441.33 calories 

ASvap14= 22.752 e.u. 

ASvapis = 22.805 e.u. 

Difference in entropy of liquids =0.313 e.u. 

Entropy of liquid 14 Pen at 63.145° = 16.346 e.u. 

Stu— 


sas. F. Giauque and J. O. Clayton, J. Am. Chem. Soc. 55, 4875 


40. Kratkey, Physik. Zeits. 34, 482 (1933); see also 
B. E. Warren et. al., J. Chem. Phys. 4, 236 (1936) ; zbid., 6, 
666 (1938) ; ibid., 7, 141 (1939). 
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TABLE III. 


A B 


OSL iquid 15 —14 


644 yap 15-14 


log Pus/P1s = 


log Pus/P 1s = [Constant 6AH and 64S] 


0.313 e.u. 
68.7 0.329 


6: =68.7 66.4 
ur =0.503 0.486 0.326 
6: =68.7 4 
8 0.325 
6 


0.337 


Or =44.3 4 
57.6 


3.33 cal. 
1.48 


(0.7230/T) —0.005822 


10.2540/T? (0.3245/T) —0.002567 


2.13 


(13.094/7?) +0.00044 (0.44206/T) —0.003278 


2.03 


(12.29/T?) —0.00044 (0.41654/T) —0.003077 


entropy of a linear oscillator. The chief difficulty 
lies in evaluating N(v). Most authors who have 
used this method have assumed some form of the 
Debye function. Eucken,® for example, intro- 
duced into the Debye function a limiting fre- 
quency which increased with temperature. Other® 
assumptions have also been made. 

In discussing differences between isotopic com- 
pounds and various other quantum mechanical 
effects,’ the ability of a molecule to move from 
cell to cell is probably of secondary importance. 
It can be arbitrarily compensated for by the 
introduction of the communal entropy.* The 
assumption will, therefore, be made that N(yv) is 
proportional to v*; i.e., a Debye frequency spec- 
trum® will be used. Since the acoustical fre- 
quencies in the liquid are lower than the corre- 
sponding frequencies in the solid, this assumption 
per se should not prove too bad. Although one 
cannot expect to predict the exact quantum 
mechanical effects using such a method, one can 
obtain some idea as to the magnitude of these 
effects. This assumption leads directly to the 
Debye expressions 


Sp=®R{4/3 D(6/T)—In (1—e-*'7)}, 


where 
6 3 x3dx 
= 
(0/T)? 
5 A, Eucken, Ber. Akad. Wiss. Berlin 22, 682 (1914). 
a ag Symposium, Trans. Faraday Soc. (September, 
7K. S. Pitzer, J. Chem. Phys. 7, 583 (1939). 
§ H. Eyring and J. Hirschfelder, J. Phys. Chem. 41, 249 
(1937); J. Hirschfelder, D. Stevenson, and H. Eyring, J. 


Chem. Phys. 5, 896 (1937). 
® P. Debye, Ann. d. Physik 39, 789 (1912). 


(7) 
(8) 


The introduction of the Planck expression in 
Eq. (6) assumes an harmonic oscillator. The 
amplitudes of the molecules vibrating in the 
liquid cells are so large, however, that this 
assumption is a very poor approximation. The 
effects of the anharmonicity must be taken into 
account. This, unfortunately, requires a theory 
of the liquid state, the lack of which has forced 
the use of very approximate models of this sort. 
These effects must consequently be discussed 
separately. This will be done briefly and qualita- 
tively later. 

Pauling’® has shown that there most probably 
is free rotation in solid nitrogen near the triple 
point. It is then safe to assume that the molecules 
are rotating freely in the liquid and that there- 
fore the rotational degrees of freedom in the 
liquid are the same as in the gas. The internal 
vibrational degree of freedom is of course en- 
tirely unexcited. @ is consequently 3. The appro- 
priate 6,4 for liquid nitrogen (14) may then be 
found by equating 


Sp=Si+—-R, 


where Sp is given by Eq. (8), Sz_: was defined 
above and is given in Table II, and R is the 
communal entropy. The value of 6:5 for liquid 
nitrogen (15) is obtained from 614 by assuming 
615= 014, where p in this case is the ratio of the 
square roots of the masses. The data in Table III, 
Model I, were obtained using these 6’s and the 
aboveequations. 5Sriquia = ThedAH ya, 
was calculated assuming that the potential ener- 
gies of the liquids are the same for the two iso- 
topes and that Exi,= Hig. Beattie’s tables" and 


LL. Pauling, Phys. Rev. 36, 430 (1930). 
uJ. A. Beattie, J. Math. Phys. (M.I.T.) 6, 1 (1926-7). 
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Fic. 1. 


the approximation formulae given by Mayer? were 
used in the calculations. Log Pu/Pis in column 


_ A was calculated from SAF °=8AH yap — T5ASvap 


= —RT In Ps. The log Pis/Pis the last column 


B wascalculated assuming that /AHy,, and 5ASvap 
are independent of temperature. This was found 
to be so experimentally and can be seen from 
Fig. 3. Figures 1 and 2 show the agreement of 
these equations with experiment. Figure 3 shows 
the deviation of the experimental curve from line- 
arity when T log P14/P1s is plotted against 1/T. 
It is to be noted that this model and the others 
to be discussed do not per se give a vapor pressure 
curve of the form found experimentally, i.e., 


In Pis/Pi1s= —(64H/RT)-+constant. 


The additional assumptions of the constancy of 
54H and 6AS° are necessary in order to obtain 
such an equation. This strong dependence on 
temperature does not, however, invalidate the 
use of these models as a means of determining 
the quantum mechanical effects. A @ is calculated 
at one temperature and from this 6 the AH yap 
and 6ASyap are calculated for that one tempera- 
ture. The temperature variation of the calculated 
64H.» and 6ASyap are of concern here only inso- 
far as they show the limitations of these models. 
It is unreasonable to expect simple models like 
those being discussed to give a complete under- 
standing of the dependence of the vapor pressures 
on the temperature. The configuration of the 
liquid changes with temperature from a quasi- 
crystalline state at the triple point to an almost 


gaseous one near the critical point. It is, there- 
fore, doubtful that anything short of a theory 
of liquids can predict satisfactorily this tempera- 
ture dependence. 

If the liquids were classical, then the observed 
5Szis-14 would have been the same as that of 
the gas; namely, 0.366 e.u. with dASyap of course 
equal to zero. Experiment has shown this not to 
be the case and a model correctly introducing all 
of the quantum mechanical effects would give the 
experimental values of 0.313 e.u. for 6Sris—14 and 
0.053 e.u. for 5ASyap. The values of 0.329 e.u. 
and 0.037 e.u., respectively, obtained for Model I 
are thus a measure of the quantum mechanical 
effects when the liquids near the triple points are 
treated as harmonic oscillators. The agreement 
with experiment is surprisingly good considering 
the crudeness of the assumptions involved. 

Pitzer’ has shown that there may be reason to 
believe that liquid nitrogen has two degrees of 
restricted rotation. If this be the case, the rota- 
tional degrees of freedom may be represented by 
either Einstein or Debye terms. When dealing 
with molecular crystals, the usual assumption” 
made is that the translational degrees of freedom 
may be treated as in Model I, whereas the re- 
stricted rotational degrees of freedom are vibra- 
tional and representable by Einstein terms. 
Making this assumption and the additional one 
that the two Einstein terms are the same, the 


results shown in Table III, Model II, and Figs. 1 


2a. R. H. Fowler and E. A. Guggenheim, Statistical 
Mechanics (Cambridge Press, 1939), p. 150. b. M. Born, 
Dynamik der Kristallgitter (Teubner, Peibaig, 1915). 
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and 2 were obtained. Lord,” however, has 
assumed that restricted rotation in the molecular 
crystal may be represented by oscillations which 
are described by a Debye function with the 
appropriate degrees of freedom. This assumption 
was here applied to the liquid nitrogen with 
@=2. The results obtained are labelled Model 
III. The Debye 6 for the rotational degrees of 
freedom of nitrogen (14) in this model (or the 
corresponding Einstein frequency in the previous 
model) was obtained by equating the Debye 
function (or Einstein function) for two degrees 
of freedom to S, and solving for @ (or u). The 
agreement between these two models is under- 
standable if one realizes that if T>6, i.e., the 
rotation, if hindered, is almost classical, the Ein- 
stein and Debye equations give the same result. 
Calculations were also made assuming that all 
degrees of freedom have the same Debye @. 
This is given as Model IV. . 

The assumption upon which these calculations 
have been made is that the molecules vibrate in 
a field of the type shown in Fig. 4a. It has been 
shown," however, that the potential is actually 


“as 


Fic. 2. 


wR, C d, J. E. Ahlberg, and D. H. Andrews, J. 
Chem. Phys. 5, 649 (1937). 
nnard -Jones and Devonshire, Proc. Roy. Soc. 
aed, 53 (1937); A165, 1 (1938); A169, 317 (1939) ; A170, 
464 (1939), and others. 


more like that in Fig. 4b. Lennard-Jones and co- 
workers" have calculated approximately the po- 
tential energy of a molecule, due to its immediate 
neighbors, as a function of position and also of 
size of available volume. They have been able 
to calculate surprisingly well the thermodynamic 
properties of liquid nitrogen. Their treatment is 
completely classical in character. Their model 
can, however, be used if quantum mechanical 
effects are introduced.'® This has not been done 
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Fic. 3. Nitrogen-experimental curve. 


here. It must, moreover, be remembered that 
these models are not pictures of the liquid state 
but are merely devices used in calculating the 
quantum effects. These models are in perfect 
agreement with the generalization'® that the 
quantum corrections always diminish the parti- 
tion function and that its absolute value in- 
creases with the mean square values of the 
forces and diminishing mass of the particles 
concerned. 

The tacit assumption has been made that the 
vibrational frequencies in the gas and liquid are 
the same. There is, however, generally a small 
shift in the internal vibrational frequencies when 
a vapor condenses and a corresponding change in 


16 J. G. Kirkwood, Phys. Rev. 44, 31 (1933); E. Wigner, 
(1982) Rev. 40, 749 (1932) ; F. Block, Zeits. f. Physik 74, 295 
(1938 K F. Herzfeld and E. Teller, Phys. Rev. 54, 912 
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the internal zero point energy of the molecule!” 
occurs. This shift is nearly always in the sense 
that the frequencies are lower in the liquid 
phase. Such a shift has been observed in the 
case of nitrogen. The Raman spectrum of liquid 
nitrogen has been studied by Sutherland!* and 
McLennan.!® They report shifts of 2326 cm 
and 2330 cm~, respectively, as compared to the 
shift of 2330.7 cm! observed in the gas.?° This 
corresponds to a higher heat of vaporization for 
the lighter isotope of about 0.28 cal., using 
Sutherland’s value. It is assumed that each 
isotope has the same fractional shift. It must be 
remembered, however, that the Raman fre- 
quencies may differ considerably from the mean 
value of the frequencies in the liquid and that it 
is the latter which determine the effective zero 
point energies. The actual effect may, moreover, 
be in the other direction.'® Herzfeld and Teller 


174, E. Bauer and M. Magat, J. de phys. et rad. 9, 319 
(1938) ; b. W. Buchheim, Physik. Zeits. 36, 694 (1935);c. R. 
Mecke and O. Vierling, Zeits. f. Physik 96, 559 (1935). 

18 G. B. B. M. Sutherland, Proc. Roy. Soc. London 141, 
535 (1933). 

19 J, C. McLennan and J. M. McLeod, ‘eo 123, 160 
1 Trans. Faraday Soc. 25, 797 (1929) 
. Rasetti, Phys. Rev. 34, 367 (1929) ; b. S. Bhaga- 
weihen ‘Ind. J. Phys. 6, 319 (1931). 


have also shown that a difference in the polariza- 
bilities and molecular volumes of two isotopes 
will affect 5AH,,, and consequently the vapor 
pressures. Tle sign of the change depends upon 
the properties of the molecules and unfortunately 
there is no such data available for the nitrogen 
isotopes. 

Furthermore, a gas below its boiling point is 
not ideal. Association can be expected in the gas 
phase. This association gives rise to a zero point 
energy ; and quantum effects, similar to those in 
the liquid, will enter in the gas. This will decrease 
the difference between the two isotopes in the 
gas and tend to make the vapor pressure of the 
lighter isotope lower. This effect, however, should 
be of second order unless the association is 
accompanied by high intermolecular vibrational 
frequencies. This is unlikely in the nitrogens. 

Another factor which tends to decrease the 
vapor pressure of the lighter isotope is the 
anharmonicity of the lattice vibrations in the 
liquid. These anharmonicities will reduce the vi- 
brational frequency in the liquid. Since the reduc- 
tion is greater for nitrogen (14) than for nitrogen 
(15) the difference in the entropies of the two 
liquids is increased and the 6ASyap is decreased. 
The effect on the heat of vaporization is in the 
same direction. One should, nevertheless, not gen- 
eralize about the effect of anharmonicity in the 
lattice vibration on 6AS,.p without knowing more 
about the molecular fields in the liquid. 


ACKNOWLEDGMENT 


I should like to take this opportunity to thank 
Professors Joseph E. Mayer, Edward Teller, and 
Harold C. Urey for a number of most helpful and 
inspiring personal discussions. 


} 722 
Q 
| 
1 
q 1 
’ 
( 
4 
| —— 
( 
a | 
q 


DECEMBER, 1942 


JOURNAL OF CHEMICAL PHYSICS 


VOLUME 10 


Potential of Silver Concentration Cells with Liquids of Low Dielectric Constant 


ANDREW GEMANT 
The Detroit Edison Company, Detroit, Michigan 
(Received September 17, 1942) 


The purpose of this work was to obtain, by means of measurements of electromotive force, 
information on the behavior and nature of ions present in liquids of low dielectric constant. 
These measurements were made on silver and silver-silver chloride concentration cells con- 
taining solutions whose dielectric constants were as low as 2.6 and whose specific conductances 
were as low as 10-" mho/cm. The results of these tests were supplemented by conductivity 
measurements on all the solutions involved. From the data obtained, approximate values for 
the dissociation constants and electrochemical transference numbers of the electrolytes for 
various electrolyte-solvent combinations were deduced. 


1. PURPOSE OF THE WORK 


N spite of considerable recent progress in the 
field of electrical insulating oils,! one of the 
fundamental questions is still unanswered. This 
question is: What is the chemical nature and the 
physical nature of the particles that conduct 
the current? It is very likely that in liquids of 
standard purity the negative particles are anions 
of complex organic acids, and the positive ones 
are hydrogen ions, but there is no direct experi- 
mental proof for these suppositions. The residual 
conduction in highly purified material is even less 
well understood. Also, it is not known whether 
ions are of molecular dimensions (say below 10-7 
cm), or larger colloidal aggregates. There is only 
indirect evidence as to the size of the charge 
carriers.” 

In view of these uncertainties it is concluded 
that further substantial progress will be achieved 
only by the establishment of new methods. It 
appears that the methods available today have 
been utilized pretty nearly to capacity and will 
not yield much more regarding fundamentals 
beyond what is already known. 

Comparing, by way of analogy, the develop- 
ments in the field of aqueous and other highly 
conducting solutions, it is striking that one of the 
most powerful tools in that field, namely, the 
measurement of electromotive forces in cells 
built of metals and liquids, had not been investi- 
gated in the realm of liquids of low conductivity. 
The reason for this is the large number of tech- 
nical difficulties involved in the measurements. 

* National Research Council, Contributions of the Chemist 


to Insulation Research (1942). 
* Andrew Gemant, Phys. Rev. 58, 904 (1940). 


Still, the author believes that it is worth opening 
up this new field, as some day the method might 
prove to be just as valuable with insulating 
liquids as with highly conducting solutions. The 
present paper is intended as a contribution along 
these general lines. 

Successful measurements with liquids having 
dielectric constants down to 17 have been carried 
out by Hartley, Harned, and co-workers.’ The 
conductance of those systems, however, is above, 
say, 10-° mho/cm, whereas the conductance of 
the solutions investigated here extended from 
10-*° to 10-" mho/cm. Thus the present investi- 
gation covered the range of so-called semicon- 
ducting liquids, it having been found impracti- 
cable to reach the range of insulating liquids, 
say below 10-" mho/cm. 

Plumley‘ reported a few measurements on cells 
with purified octane. The present author® carried 
out a research on cells composed of copper elec- 
trodes and liquids such as kerosene, xylene, etc., 
with and without copper oleate. The construction 
of the cells and boundaries and the composition 
of the liquid phases were such that the research 
could be considered only as preliminary. The 
chief conclusions were of a qualitative nature, 
viz., (1) such cells possess a definite e.m.f., 
measurable by means of a potentiometer, (2) the 
magnitude of the e.m.f. depends on the concen- 
tration of the copper salts in the solutions; thus 
it can be interpreted as electrochemical potential. 


3D. A. MaclInnes, The Principles of Electrochemistry 
(Reinhold, New York, 1939). 

4H. J. Plumley, Phys. Rev. 42, 140 (1937). 
a _ Gemant, Trans. Electrochem. Soc. 78, 49 
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The present paper carries this work one step 
further. The cells, the liquid boundary, and the 
composition of the solutions have all been better 
defined, so the attempts to evaluate the test data 
in a quantitative manner have been justified. 
Still, it should be kept in mind that the work is 
of an exploratory type, that the field is quite 
new as yet, and that in consequence the experi- 
mental data should not be considered as final. 


2. TECHNIQUE OF MEASUREMENTS 


For all e.m.f. measurements in this work, 
silver was used as the electrode, and the cells 


were of the following type: concentration cells 


with transference, i.e., a diffusion potential was 
allowed to build up at the boundary of the two 
liquid phases. The electrodes were partly bright 
silver, thus of the ‘‘first-order’”’ type in which the 
electrode reaction is: Ag=Ag*++e-; and partly 
silver coated with silver chloride, thus of the 
“second-order” type in which the reaction is: 
Ag+Cl-=AgCl+e-. The latter is in reality a 
reversible chlorine electrode. 

Strips of silver foil were used as electrodes. 
Analysis of the material revealed 100 percent of 
silver within the accuracy (about 0.1 percent) of 
the analysis. Coating with AgCl was done electro- 
lytically in 0.1” HCI solution according to a 
method described by A. S. Brown.® 

In some of the tests, chiefly in the first stages 
of the work, open glass beakers were used as half- 
cells, the bridge between being established by 
means of a semicircular quartz capillary of about 
1-mm diameter. In view of the volatility of the 
solvents, the need for a closed cell was soon 
recognized, and the one shown in Fig. 1 was used. 
It consists of a U-tube with a stopcock (A) of 
wide bore in the middle, which is closed only 


A 
Fic. 1. Diagram of glass cell for e.m.f. measurements. 


6 A. S. Brown, J. Am. Chem. Soc. 56, 646 (1934). 
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during certain phases of filling of the tube. The 
liquid of higher concentration fills the left branch, 
and the central right branch up to the opening 
(B), the rest of the right branch being filled with 
the liquid of lower concentration. The diffusion 
boundary is formed at the opening (B) of the 
narrow tube within the wider one. From platinum 
wires (C) are hung silver strips (D). Ground glass 
stoppers (Z) close the vessel. During measure- 
ments the stopcock (A) is, as already mentioned, 
open and the cell is placed upon a block of paraffin 
wax within a grounded metal box. Connections 
from the platinum wires lead to an electrometer 
tube device of standard construction employing 
an FP54 tube. It was estimated that cells of an 
internal resistance up to 10% ohms could be 
measured by this equipment. 

The solvents and electrolytes used in the work 
were the best commercial preparations, mostly 
chemically pure, but no further laboratory purifi- 
cations were carried out, neither were they pur- 
posely dried before use. The reasons are as 
follows: 

1. The more the solvents are purified and 
dried, the higher are the specific resistances 
attained. This would have increased the diffi- 
culties of measurement considerably. Thus it was 
thought better to compromise. It is important to 
keep in mind that the materials were of medium 
purity and in equilibrium with average labora- 
tory humidity. 

2. Since the work is of an exploratory nature, 
it was not deemed essential to undertake ex- 
cessive purification. That would be more justified 
in cases of confirmatory or precision measure- 
ments in a field of already established technique. 

From among the electrolytes used, only two 
were prepared in the laboratory. Tetrabuty]- 
ammonium chloride was prepared from the corre- 
sponding iodide via the hydroxide that was ob- 
tained by adding silver oxide. Silver picrate was 
produced electrolytically, silver being used as 
anode in a 0.1 normal picric acid solution in ethy] 
alcohol. Yellow crystals soon covered the anode, 
acting as a barrier to the further passage of cur- 
rent, and lowering the latter’s initial value in a 
ratio 1 : 20. 

In order to supplement the e.m.f. data and 
help in the evaluation of them, parallel con- 
ductivity measurements were carried out on all 


| 
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the solutions that were subjected to e.m.f. meas- 
urements. An all-glass-and-platinum cell of closed 
construction, devised by J. C. Balsbaugh at 
M.I.T., was used for this purpose. For the high 
power factor samples a General Radio impedance 
bridge, operating at 1000 c.p.s., was used; for 
the low power factor samples a high precision, 
variable-frequency bridge was used. Frequencies 
of from 1000 to 3000 c.p.s. were used. The fre- 
quency at which the conductivity of solutions is 
measured is irrelevant, at least in the acoustical 
range. A variation of conductivity with frequency 
in this range could be expected only in hetero- 
geneous systems. 


3. EVALUATION OF CONDUCTIVITY DATA 


First, the way in which conductivity data were 
interpreted will be shown. Electrolytes in alco- 
holic and similar solvents in a range of D (di- 
electric constant) of around 20-30 are usually 
weak electrolytes that are only partially dissoci- 
ated. In applying Debye-Hiickel’s theory, it is 
possible to determine both the dissociation con- 
stant K and the limiting conductance Ao from a 
plot of conductivity o vs. molar concentration c. 
How this is best done was shown by Fuoss and 
Kraus.’ 

This particular method can no longer be 
applied in the low dielectric constant range 
(around 3) of the present research, since the 
degree of dissociation and equivalent conductance 
are too small to allow a reasonable extrapolation 
for zero concentration. Thus A» cannot be deter- 
mined in these cases. On the other hand, one has 
the advantage that the ionic concentration is so 
small as to make the activity coefficient practi- 
cally unity. In other words, one need not con- 
sider the Debye-Hiickel effect. 

A relatively simple plot of measured data then 
will reveal whether the electrolyte behaves as a 
normal weak electrolyte, or deviates from this 
behavior by formation of complex ions, such as 
were found by Fuoss and Kraus in solutions of 
tetraisoamylammonium nitrate in dioxane-water 
mixtures. This simple plot is suggested by the 
approximate equation, valid for very weak elec- 


trolytes: 
A’c=Ac?K. (1 ) 


1933) M. Fuoss and C. Kraus, J. Am. Chem. Soc. 55, 476 
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TETRABUTYLAMMONIUM CHLORIDE 


! 


| 


DIELECTRIC CONSTANT = 3.2 
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0.2 0.4 0.6 0.8 


CONCENTRATION MOLS PER LITER, 107° 


Fic. 2. Conductivity-concentration plot for tetrabutyl- 
ammonium chloride solutions in 3 percent water in di- 
oxane. 


If, then, 1/A? is plotted against c, a straight 
line going through the origin should be obtained 
in case of a normally behaving weak electrolyte, 
the slope of which line is equal to 1/Ac?K. 

A word should be added here as to the con- 
ductivity of the solvent in question without any 
added electrolyte. The correct procedure prob- 
ably is to deduct this value from the conductivi- 
ties measured and in most (although not all) 
cases this correction was applied before the 
curves were plotted. However, in case of higher 
conductivities this correction is too small to 
matter, and it was found by frequent checks that 
even in the low conductivity range in question 
the conductivity of the solvent as a rule was 
below 20 percent of that of the weakest solution 
investigated. 

A typical plot is shown in Fig. 2, referring to 
tetrabutylammonium chloride in a solvent of 3 
percent water in dioxane (dielectric constant 
= 3.2). Table I shows the measured data of spe- 
cific conductivity o and the values of A=10*0/c. 
The plot shows a very good approximation to a 
straight line at low concentrations, proving that 
the electrolyte is a normally behaving weak 
electrolyte in the solvent with D=3.2. The slope 
of the line is 4.610’. The point representing 
the highest concentration lies definitely below the 
straight line, as was observed quite generally. 
This is due, as mentioned, to ionic association. 

While a separation of the term Ao?K into its 
factors is not possible from these data, an approxi- 
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LITHIUM CHLORIDE 


DIELECTRIC CONSTANT= 2.8 


0.1 0.2 0.3 
CONCENTRATION MOLS PER LITER ,10°> 


Fic. 3. Conductivity-concentration plot for lithium chloride 
solutions in 2 percent water in dioxane. 


mate estimate of the value of K can be made by 
taking Ap~60, which is the correct order of mag- 
nitude as known from data where a direct deter- 
mination is feasible. With this assumption K ~6 
X10-" is obtained. While this value is only esti- 
mated, it is probably quite close, say within 25 
percent of the correct value. 

A second example is shown in Fig. 3. It refers 
to lithium chloride solutions in a solvent of 2 
percent water in dioxane (D=2.8). The measured 
data are shown in Table II. 

The plotted data lie around a straight line 
although not so closely as in Fig. 2. 

The slope of the line in Fig. 3 has the value 
8.4X10"° and, as before, K~3.3X 10-5. 

In certain instances, particularly in some of 
the silver picrate solutions (to be discussed later) 
plots showed an upward trend as the concentra- 
tion increased. In such cases the drawing of a 
straight line seemed less justified than in the 
majority of plots. 


4. EVALUATION OF E.M.F. DATA 


Next we come to the e.m.f. measurements. 
Here the following general behavior was ob- 
served. Above D =3 the measurements were rela- 
tively easy to make and equilibrium was reached 
fairly quickly. Below D=3 the e.m.f. data, read 
on the potentiometer, often changed rapidly in 
the beginning, then slowly, and finally reached 
equilibrium. This process took 30 minutes in 
some cases. It could be established fairly cer- 
tainly that this was not due to the characteristics 
of the apparatus, but rather to the properties of 
the interphases. When the final value was ob- 


tained, reversal of the direction of the cell caused 
little or no variation of the measured data. 

It seemed as if D=2.6 was about the limit 
down to which reliable e.m.f. data were obtain- 
able. Below this value the figures became erratic 
and not reproducible. The reason for this is not 
definitely established, but it seemed that the 
sluggishness of the electrode reaction was the 
cause. Intentional polarization of the electrodes 
during measurements seemed to promote estab- 
lishment of an equilibrium. It is quite possible 
that further research, aimed at using certain 
physical or chemical agencies to overcome the 
slowness of the electrode reaction, will facilitate. 
work in this range. 

While the accuracy of the reading was 1 mv, 
the reproducibility did not reach this limit. The 
order of magnitude of the potential differences 
was around 30 mv and the reproducibility of this 
value was about +10 percent. By preparing 
entirely fresh samples and repeating the measure- 
ments, occasionally even greater deviations were 
observed. This type of behavior is, however, 
quite general with dielectric liquids, so it is not 
surprising here. The physical status depends on 
too many single factors, not all of which are, in 
most cases, under control. 

In plotting e.m.f. data it should be recalled 
that the voltage E of a concentration cell com- 
posed of monovalent ions is 


E=2t(RT/F) log (f1y1/fev2), (2) 


where y=ionic concentration, f=activity coeffi- 
cient, indices 1 and 2 referring to the two half- 
cells, R=gas constant, 7=temperature in K°, 
F=Faraday constant, and t=transference num- 
ber of the non-reversible ion. 

We now make use of the data obtained from 
conductivity measurements. Since 


y=(Ke)}, 
we have with K ~10-” and c~10-, y~ 10-8, and 
the factors f become very nearly unity, if the 


TABLE I. Conductivity data for N(C4Hg)«Cl in 
3 percent water in dioxane. 


A 
24 X10-3% 


15.5 
8.0 
5.0 
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equations of Debye-Hiickel are applied. More- 
over, the ratio of ionic concentrations can be 
replaced by the ratio of the square roots of the 
molar concentrations. Hence we have, instead of 
Eq. (2): 

E=59 logio (¢1/C2), (3) 


giving E in mv, the factor RT/F at room tem- 
perature having the value of about 59 mv. 

If, then, E is plotted vs. logio ¢1/c2, a straight 
line going through the origin should be obtained, 
having a slope of 59¢. From this, the transference 
number ¢ is known. 

As an example, measurements paralleling those 
shown in Table I are given in Table III. The 
corresponding plot is given in Fig. 4. In spite of 
the noticeable deviation of the measured points 
from the straight line shown, the e.m.f. data 
prove satisfactorily the validity of Eq. (2) in this 
range. This means that the fundamental laws of 
electrochemistry seem to hold even for solvents 
of a nature approaching that of insulating liquids. 
The slope of the line is 23, and thus the trans- 
ference number of the organic cation in its 
chlorine salt is 0.39. 

In supplementing Fig. 3, Table IV gives the 
e.m.f. data with the LiCl solutions, and Fig. 5 
the corresponding plot. In spite of the low di- 
electric constant (2.8) the straight line require- 
ment is well fulfilled, the slope being 22 and the 
transference number for the lithium ion in LiCl 
in this solvent 0.37. 


5. NUMERICAL RESULTS 


Proceeding in the manner indicated in the last 
two sections, a systematic study was carried out 


TaBLe II. Conductivity data for LiCl in 
2 percent water in dioxane. 


c A 
3.2 10-5 23 X 0.73 X 10-3 
1 x10 32 0.32 
10-4 64 0.20 


TABLE III. E.m.f. data for N(C4Hg9)4Cl in 
3 percent water in dioxane. 


ce E, mv 
3.2 10-3 10-4 34 
3.2 3.2 23 
3.2 10 13 
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Fic. 4. E.m.f.-concentration plot for tetrabutylammonium 
chloride solutions in 3 percent water in dioxane. 


on a number of electrolyte and solvent com- 
binations. 

Two types of solvents were chosen; either 
xylene with a small percentage of ethyl alcohol, 
or dioxane with a small percentage of water, or 
ethy! alcohol. 

The electrolytes were also of two types. In one 
series the anion was chloride and correspondingly 
the electrodes were silver-silver chloride. The cor- 
responding cations were: tetraethylammonium, 
tetrabutylammonium, 3-9 bisdimethylaminophe- 
nazothionium (methylene blue), lithium, and 
hydrogen. In the other series the cation was 
silver, and the electrodes silver, too; the anions 
were picrate and nitrate. 

The deduced values for the dissociation con- 
stant K and the cation transference number ¢, are 
summarized in Table V. It should be restated 
that the former is only an estimate, although 
probably a close one. The latter is uncertain 
within the experimental errors, and because of 
the questionable applicability of Eq. (3). 

Figures 6 and 7 show the results for silver 
picrate as electrolyte. The abscissae are dielectric 
constants, and the ordinates the logarithms of dis- 
sociation constant and cation transference num- 
ber, respectively, for the three different solvent 
systems. 


6. DISCUSSION OF RESULTS 


In discussing the results, the chief conclusion 
to be stressed is that e.m.f. measurements in the 
range of dielectric constants as low as 2.6 and 
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Fic. 5. E.m.f.-concentration plot for lithium chloride 
solutions in 2 percent water in dioxane. 


specific conductivities as low as 10-! mho/cm 
are possible and that the data are sufficiently 
reproducible to allow a quantitative evaluation 
of the results. Whether the limit in dielectric 
constant has been reached or not cannot be 
decided from the material on hand. A continua- 
tion of the research along these lines is planned. 

A question of applicability of the method to 
physico-chemical determinations of both scien- 
tific and technical nature, as inferred in the 
introduction, cannot yet be answered definitely. 
It is apparent that more work of an exploring 
and fundamental nature will have to be carried 
out before the method can be tried, for instance, 
for the determination of ionic concentrations in 
insulating liquids. 

Regarding the specific results obtained in this 
work on the dissociation constant, it can be seen 
that the dielectric constant of the solvent has a 
primary influence upon its magnitude, a fact 
well known from data on solutions of higher 
conductivity and accounted for by all theories 
based on the electrostatic attraction of ions. 
Apart from this, it can be seen that the alcohol- 
xylene system has a higher dissociative power 
with all electrolytes than has the water-dioxane 
system. This indicates a specific influence on the 
part of the solvent. 

In this connection it should be pointed out 
that mixtures of solvents probably cannot be 
treated as if they were homogeneous, since a 
separation near the ions will take place, the mix- 
ture being enriched with regard to the component 
with the higher dielectric constant. Homogeneous 
liquids, however, do not appear to be suitable as 
solvents for electrolytes. Thus chloroform with 


D=4.7 was tried but its solubility for electrolytes 
was found to be insufficient. It appears that small 
amounts of admixtures of higher polarity are 
quite essential for ionization and solubility (the 
two being to a certain extent parallel) in insu- 
lating liquids. Sticher and Piper® arrived at a 
similar conclusion in their study of insulating 
liquids bombarded by electrons. 

The nature of the electrolyte will also affect the 
dissociation constant as can be seen by comparing 
the different electrolytes in the system com- 
prising 3 percent water in dioxane. Organic 
electrolytes have a higher dissociation constant 
than inorganic ones, and the more complex the 
molecular structure the greater is the difference. 
This may be seen also from the high K values of 
methylene blue. This circumstance might yield a 


TABLE IV. E.m.f. data for LiCl in 2 percent 
water in dioxane. 


E, mv 


34.5 
20.5 
11.0 


TABLE V. Deduced values (approximate) of dissociation 
constant K and cation transference number ¢, in solutions 
of low dielectric constant. 


Dielectric 
Electrolyte Solvent constant K 


N(C2Hs)«Cl 10% alcohol in 3.2 16 X10-" 
xylene 

N(Ce2Hs)sCl 5% ethyl alcohol in 25 X1075 

N(CoHs)4Cl 2.5 X10-2 


xylene 
3% water in dioxane 
N(C2Hs5)4Cl 60 X10715 
N(C4Ho9)4Cl 6 X10-2 


2.7 
3.2 
2% water in dioxane 2.8 
3% waterindioxane 3.2 
N(CsH¢)4Cl 2% waterindioxane 2.8 120X10-% 0.40 
CisHisNsSC1-3H20 10% alcohol in 3.2 0.4 X107% 
xylene 
CisHisN3SC1-3H20 5% = alcohol in 2.7 0.61072 
xylene 
LiCl 10% ethyl alcohol in 3.2 9X<10-2 
2.7 
3.2 
2.8 
3.2 
3.2 
3.1 


xylene 
LiCl 5% ethyl alcohol in 3.5 X107- 
xylene 
LiCl 3% water in dioxane 0.2 X10- 
LiCl 2% water in dioxane 3.5 X10-% 
HCl 10% ethyl alcohol in 5.5 X1072 
xylene 
HCl 0.7 X10712 
AgO-CcéH2(NO2)3 10% alcohol in 60 X10745 
ioxane 
AgO -CeéH2(NO2)s 8% alcohol in 2.85 4.51075 
10% in 3.2 5 X1072 
xylene 
AgO -Ce6éH2(NO2)s 5% ethyl alcohol in 2.7. 8.51075 


xylene 

AgO -CeH2(NOz)3 4% alcohol in 2. 1 X10745 
xylene 

AgO -CeéH2(NOz2)s 3% waterindioxane 3.2 0.7 X1072 

AgO -CsH2(NO2)s 2% water indioxane 2. 20 X10-% 0.1 

AgO -CeéH2(NO2)s 1% waterindioxane 2. 45 X1078 

NO. 5% waterindioxane 3. 23 X10-2 
4% water in dioxane 
3% water in dioxane 


3% water in dioxane 


— and J. D. Piper, Ind. Eng. Chem. 33, 1567 
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Fic. 6. Dissociation constant vs. dielectric constant for 
silver picrate solutions in various solvents, as shown. 


clue that would help in understanding the rela- 
tively high conductivities obtained with bom- 
barded insulating liquids. Balsbaugh,® too, ex- 
plains the high power factor of his hydrocarbon 
oils, oxidized in the presence of Cu with small 
amounts of O, by electrolytic dissociation. 

From the straight line relation of the con- 
ductivity plot it follows that the ions are single 
ones in the lower range of concentration (below 
10-* molar), but above this range multiple 
ions in the sense of Fuoss and Kraus might 
appear. 

Concerning the figures on transference numbers 
(ratio of cation mobility to the sum of cation and 
anion mobility), it is noted that, with few excep- 
tions, they are less than 0.5. While this value 
seems to be plausible for electrolytes having 
large organic cations, it does not seem quite 
plausible for silver picrate, where one would 
expect ¢, figures larger than 0.5. 

Several possible explanations may be advanced. 
One is that the use of Eq. (2) in the same form 
as known from aqueous solutions is not correct. 
It is possible that a correction in this equation 
will have to be introduced that will correspond- 
ingly change the transference figures. 

Another possibility, and one that deserves con- 


C. Balsbaugh, A. G. Assaf, and J. L. Oncley, Ind. 
sha ‘Chem. 34, 93 ¢ (1942). 
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Fic. 7, Cation transference number vs. dielectric constant 
for silver picrate solutions in various solvents, as shown. 


2.6 


sideration, is that while Eq. (2) may be correct, 
the use of Eq. (3) is not justified. One way of 
avoiding this difficulty would be to compute the 
ratio behind the logarithm directly from the 
measured conductivity data. Because of some 
objections, of both experimental and theoretical 
nature, that could be raised against such a pro- 
cedure, the method of using Eq. (3) in computing 
the data in Table V was given preference, but in 
future work of this type this point will have to be 
checked. While the results derived according to 
both methods generally do not differ much, 
there is a divergence in some cases. Thus ¢, for 
silver picrate in water-dioxane mixtures would 
become as high as 0.40 and that for silver nitrate 
in the same mixtures as high as 0.50, if this 
second method of computation were used. 

A third possibility is that the deduced figures 
are correct and explainable by a certain type 
of ionic association. One can assume that the 
undissociated molecules of silver picrate form 
aggregates of two (or more) single molecules, and 
that the silver ions become attached to these 
aggregates, thus attaining a relatively large size 
and a small ¢.. The author has computed the 
equilibrium obtained under these assumptions, 
and it can be shown that plots in the manner of 
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TABLE VI. Ionic dimensions in A units, calculated for 
silver picrate solutions in alcohol-xylene. 


D a Ye Ya 
5 3.2 6.1 3.2 2.9 
aa 5.9 4.2 1.7 

2.6 5.7 4.6 Ei 


Figs. 2 and 3 should have an upward concavity, 
: as has indeed been observed experimentally on 
several solutions of silver picrate. The explana- 
tion is also in keeping with the observed de- 
"i crease of ¢, with decreasing dielectric constant 
! (Fig. 7), because formation of aggregates is 


favored in solvents of low D values. These com- 
putations are not, however, shown here, as they 
need further verification before publication would 


d be justifiable. 

One could attempt to use values of K and ¢, 
F to compute ionic radii. The theory of Bjerrum,’ 
developed later by Kirkwood," establishes a rela- 


W. M. GARRISON AND M. BURTON 


of ¢. in conjunction with Stokes’ equation can be 
used to compute ionic radii. A few figures are 
shown in Table VI, using Bjerrum’s equations. 
They refer to silver picrate in alcohol-xylene and 
include the closest distance of approach a and 
the ionic radii r, and r, in A units. 

An interesting conclusion is that the dimen- 
sions are several angstroms. Thus the ions appear 
to be of molecular size. While the trend for r, is 
in keeping with the formation of aggregates, rq 
shows an unaccountable decrease with a decrease 
of dielectric constant. The figures substantiate 
only partly the explanation suggested above, 
based on large positive ions, and indicate that 
the transference numbers may have to be revised. 

In any case the data obtained seem to be of in- 
terest and warrant a continuation of the research. 

The author wishes to thank several of his 
colleagues in the Research Department of The 
Detroit Edison Company both for discussion and 


tion between K and the distance of nearest ionic 
approach. When the latter is known, the figures 


10 J, A, Kirkwood, J. Chem, Phys. 2, 767 (1934). 


assistance, in particular. Howard A. Boltz for 
having built the necessary electrical measuring 
apparatus. 
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The photolysis of propionaldehyde at \A2537, ~2900, and ~3200A has been studied by 
i the Paneth mirror method. The alkyl radicals produced are exclusively ethyl; no atomic 
i hydrogen can be detected by the guard mirror method in amounts exceeding 2 percent of the 
H total number of mirror-active particles present. The yield of C:H; radicals has been deter- 
mined and compared with the yields of CO, Hz, and C2H¢. The results are most readily in- 
terpreted in terms of competing primary decomposition processes producing free radicals in 
j one reaction and ultimate molecules in another. Whereas the former reaction increases in im- 
portance at shorter wave-lengths, the latter becomes less significant. On the other hand, an 
explanation for the results can be given on the basis of a mechanism involving only the pro- 
duction of free radicals in the primary process at all wave-lengths. The role played by excited 
propionaldehyde molecules (of life which may be >10-* sec.) and by freshly formed (energy- 
rich) C;H; and HCO radicals on their first collision is also considered. The mechanism of energy 
transfer within the propionaldehyde molecule is discussed and it is shown that the products of 
the primary decomposition are probably determined by geometric considerations and by the 
relative heights of the dissociation levels involved in the normal state. 


it has been generally recognized since the work 
of Leighton and Blacet! that such other sub- 


T is not possible to write a single over-all 
equation for the photolysis of propionalde- 
hyde. Although the principal equation is 


CO, 


1 P. A. Leighton and F. E. Blacet, J. Am. Chem. Soc. 54, 
3165 (1932). 
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stances as hydrogen, butane, and polymer are 
produced in the process. It was such products 
that created the suspicion that free radicals were 
involved and that prompted the repeated efforts? 
to detect them. 

Recent work by May, Taylor, and Burton’ 
employing the Paneth mirror method has suc- 
cessfully shown that free radicals are produced 
in the photolysis of well-purified propionalde- 
hyde. Although the same authors were able to 
demonstrate by the guard mirror method‘ that 
the production of free hydrogen atoms at 2537A 
was certainly less than 2 percent of the total free- 
radical production (and perhaps non-existent), 
they were not able to decide conclusively the 
nature of the radicals; i.e., whether they were 
methyl, ethyl, or a mixture of the two. The 
research reported herein was intended initially 
to clear up this portion of the problem. In the 
course of the work data were obtained which are 
here used to indicate the probable nature of the 
primary process. It will be shown also that the 
investigation throws some light on the mecha- 
nism of energy transfer within the molecule, as 
well as on possible secondary processes. 


EXPERIMENTAL 


The work involved in this investigation may 
be summarized as follows: (1) check of the half- 
life of radicals produced by the photolysis of 
propionaldehyde, (2) analytical determination of 
the radicals involved, (3) analytical determina- 
tion of the free radical to carbon monoxide ratio, 
and (4) analytical determination of the more 
volatile gases produced in the photolysis, both in 
the presence and the absence of lead mirrors. 
Except for part (1) the effects of variation of 
wave-length and temperature were also examined 
in each case. In addition, and incidental to the 
main purposes of this work, an effort was made 
to detect production of free hydrogen atoms at 
the longest wave-lengths employed, namely, at 
about 3200A. 


Light Sources 


No effort was made to maintain real constancy 
of the light sources other than to fix their posi- 


2T. G. Pearson and J. H. Purcell, J. Chem. Soc. 1151 
(1935); W. West, J. Am. Chem. Soc. 57, 1931 (1935). 
a one” Taylor, and Burton, J. Am. Chem. Soc. 63, 249 
1). 
4M. Burton, J. Am. Chem. Soc. 58, 1645 (1936). 
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tions accurately by stops. Unless otherwise indi- 
cated, the light intensities at the irradiated zone 
were probably constant within about 10 percent 
for any given wave-length. 


2537A 


The source was a low pressure spiral lamp of 
the type previously described.* The Hanovia 
Chemical and Manufacturing Company, manu- 
facturers of the lamp, state that 98 percent of 
the ultraviolet light is concentrated at 2537A, 
leaving a possible maximum of 2 percent at 
1850A. The energy output was controlled by a 
constant current transformer, the primary of 
which was held constant by means of an Adjust- 
a-volt transformer. 


2900A 


The light source was a General Electric H-6 
lamp in a quartz envelope. Spectrographs showed 
that when a chlorine-mercuric chloride filter’ 
was employed (as in this work) the range trans- 
mitted was 2800-2967A with a maximum at 
~2900A corresponding to about 80 percent of 
the energy. 


3200A 


The light source was the H-6 lamp described 
above with a Pyrex filter in place of the chlorine- 


mercuric chloride one. Photographs of the 
transmitted mercury lines showed complete 
absorption by the Pyrex filter below 3000A with 
a broad maximum of transmitted light at ~ 3106— 
3291A. The longer wave-lengths transmitted are 
known to be without substantial effect in the 
photolysis of the aldehydes. 


Temperature Control 


For the 30° runs an electric fan manually con- 
trolled was used to keep the temperature con- 
stant to +1°. 

For the 100° runs at 2537A, the spiral lamp 
was placed inside an asbestos shield just large 
enough to enclose the lamp and the irradiated 
zone. Heat from a micro gas burner was used to 
supplement the heat of the lamp and to maintain 
constancy of temperature. In the work at 2900 
and 3200A, steam escaping through two nozzles 


5 Cf. E. J. Bowen, J. Chem. Soc. 76 (1935). 
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onto the irradiated zone, delimited by asbestos 
shields, maintained the temperature at 100+1°. 

Except for the irradiated zone the reaction 
tube in each case was at room temperature, 
usually 18-25°. 


Chemicals 


Propionaldehyde was purified, stored, and pre- 
pared for use according to the method described 
by May, Taylor, and Burton.* The corrected 
boiling range of the fraction used was 48.1-— 
48.3°C which may be compared with the range 
47.8-48.0°C reported by them. 

Mercuric bromide used in separating the 
diethyl mercury from the excess propionaldehyde 
was purified by reprecipitating Eimer and 
Amend C.P. mercuric bromide from absolute 
alcohol. After drying, the mercuric bromide was 
kept at 100°C under vacuum for several hours. 


1. HALF-LIFE EXPERIMENTS 


In this part we worked with apparatus like 
that already described* and, following that 
method, obtained a value for the half-life (i.e., 
time required for half the radicals to disappear) 
of 4.0X10-* sec. which, considering the uncer- 
tainty of the visual method for following (lead) 
mirror removal, is in good agreement with the 
value of 3.9X10-* sec. previously given. The 
actual data obtained are not here reviewed. 
Since the earlier result is based on more data, we 
believe it to be more reliable. 


2. DETERMINATION OF THE NATURE OF THE 
RADICAL PRODUCED 


In this part we followed the technique em- 
ployed by F. O. Rice and his co-workers.* The 
high vacuum techniques already established were 
maintained throughout this work. Free radicals 
produced by streaming the propionaldehyde at 
an input pressure of 2.5 mm through the irradi- 
ated zone, were allowed to react with a mercury 
mirror. The resultant product was collected, 
together with unaffected propionaldehyde, in a 
reservoir cooled at liquid-air temperatures. The 
amount of mercury product collected in any one 
prolonged run was insufficient for analysis. Con- 


6 F, O. Rice and K. K. Rice, The Aliphatic Free Radicals 
(Johns Hopkins Press, Baltimore, 1935). 
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sequently, after closing and opening the appro- 
priate stopcocks, the propionaldehyde was sepa- 
rated from the mercury alkyl compound by 
raising the temperature of the reservoir to that 
of a chlorobenzene mush (~ —45°) and distilling 
the propionaldehyde back to the storage reservoir 
maintained at liquid-nitrogen temperature. Sev- 
eral runs using ~250 g of propionaldehyde 
were made in this way at any one condition of 
light and temperature (in the irradiated zone). 
The trap was then removed for recovery and 
analysis of the mercury alkyl. 


Purification and Identification of Alkyl 
Mercury Compounds 


The trap containing the salts was removed and 
attached to a cold finger as described by Rice 
and Rice.* The trap was evacuated and heated to 
100°C in a water bath. Cold water was passed 
through the cold finger and the alkyl mercuric 
bromide sublimed, condensing on the finger. The 
system was always closed during an actual 
sublimation, in this way insuring equilibrium and 
preventing fractionation. After the sublimation, 
the trap was removed and cleaned of excess mer- 
curic bromide. The apparatus was reassembled 
and the alkyl mercuric bromide was resublimed 
back onto the bottom of the trap, by passing 
steam through the finger while the trap was kept 
at 0°C. This procedure was repeated ten times. 

Identification of the alkyl mercuric bromide 
was made by analyses for carbon, hydrogen, and 
bromine using the micro technique of Niederl and 
Niederl? with some modifications. The micro 
combustion tube used in the carbon-hydrogen 
analyses was of the usual design but contained 
~3.5 cm of gold sponge*® near the exit end to 
remove any mercury vapor. Details of the anal- 
ysis of the alkyl mercury compounds are given in 
Table I. The ratio H : C : Br is included because 
the possibility of a slight contamination by mer- 
cury must be considered. The ratio calculated in 
this manner is based on the carbon analyses, 
which are the most accurate. Evidence for ethyl 
radicals as the only radicals formed in any 
measurable amount seems conclusive. 


™Niederl and Niederl, 
Analyses (John Wiley and Sons, New York, 1938). 
8 M. Furter, Mikrochemie 10, 66 (1931). 
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TABLE I. Details of analysis of alkyl mercury 
bromide samples. 
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length Percent Ratio 
Sample H :€: Or 
2 2537 1.56 7.66 23.7 0.204 1 3.10 
4a 2537 1.75 7.99 26.4 219 1 3.29 
5a 3200 1.77 7.76 25.0 228 1 3.19 
C.:H;HgBr 1.61 7.76 25.8 .208 1 3.32 
(theor.) 
CH;HgBr 1.02 4.06 27.1 sao 61 «(66ST 
(theor.) 


3. FREE RADICAL: CARBON MONOXIDE RATIO 


The method of counting free radicals, de- 
veloped by Feldman, Ricci, and Burton,? was 
followed in this work. In brief the method is to 
deposit an excess of lead as a sensitive mirror, to 
pass the gas containing free radicals (e.g., par- 
tially photolyzed propionaldehyde) over the 
mirror and to collect the products, and to 
measure analytically the amount of lead trans- 
ported. On the basis of the assumption that the 
lead was transported exclusively as lead tetra- 
ethyl, the amount of free ethyl reacting with the 
lead during the period of the photolysis could 
then be determined. It will be shown in the dis- 
cussion that a satisfactory standard value to 
which the. production of free radicals can be 
related is the yield of carbon monoxide in the 
experiment; the method for determining the 
latter is indicated in the next part. In Table III 
the yield of ethyl is expressed as a C2Hs/CO 
ratio. The experiments listed are numbered in 
chronological order. All experiments actually 
completed are included in the table. 


4. MORE VOLATILE PRODUCTS 


When the gaseous products were to be 
analyzed, the high vacuum required during a run 
(of about five minutes) was maintained with a 
mercury piston. At the end of the run, when the 
flow of propionaldehyde vapor had been stopped 
and while the receiver was still maintained at 
~77°K, pumping was continued until all the 
fixed gases produced had been transferred to a 
gas analysis system of the Saunders-Taylor 
type.’° In order to insure complete removal of 


— Ricci, and Burton, J. Chem. Phys. 10, 618 


(1 
© K. W. Saunders and H. A. Taylor, J. Chem. Phys. 9, 
616 (1941), 
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such gases from the propionaldehyde, pumping 
was interrupted several times while the solid in 
the receiver was repeatedly melted and refrozen. 
This process was continued until no further gas 
could be drawn off in this way. 

When the analysis of the first fraction, con- 
sisting of CO and He, had been completed, a 
second fraction was drawn off in a similar way 
at the temperature of sec. butyl chloride mush 
(~138°K)." No effort was made to draw off less- 
volatile products, such as butane, because of 
the difficulty of separation from the propion- 
aldehyde.! 

Although, in general, the Saunders-Taylor 
technique was followed, slight modifications 
proved desirable. In the original technique the 
procedure is simply to burn the gas on a hot 
platinum wire, drawing the products through 
gold leaf (de-ozonizer), and then to absorb the 
water and CO, formed in magnesium perchlorate 
(“anhydrone’’) and freshly dried KOH, respec- 
tively. A difficulty with the technique involves 
the latter substance, which must not be too dry, 
for the rate of absorption of CO, will then become 
too slow. We found it necessary to use anhydrone 
for drying after the absorption of CO, in KOH 
even when supposedly dry KOH was used. The 
explanation is simply that this apparatus is 
capable of measuring microscopic amounts of 
water vapor and, when they are not carefully 
removed, they have a measurable effect con- 
siderably disturbing to the accuracy of the 
analysis. Having discovered this, we thereafter 
used only semi-dry KOH and dried the residual 
gas. 

In the case of ethane, Saunders and Taylor 
found it necessary to apply a correction factor, 
since they consistently obtained analyses indica- 
ting C2H¢.s. We have found that in this latter 
case the difficulty is at the beginning of the 
analysis. Traces of moisture present in the 
oxygen must be removed before it is measured 
and mixed with the sample prior to combustion. 
Under such circumstances “‘pure”’ ethane gives a 
value of m in C,Hens2 ranging from 1.91 to 2.10 
(depending on the method of calculation), and 
no arbitrary correction need be applied. 

In some analyses of the first fraction, a 


" Cf. D.S. Herr and W. A. Noyes, Jr., J. Am. Chem. Soc. 
54, 3165 (1932). 
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TABLE II. Details of ethane analyses. 


nin ConHa +2 


Obs. vol. Calc. vol.* CO2z COz 
Expt 1 Ve Vi V2 
128 4.33 4.16 1.96 2.04 
129 6.42 6.26 2.05 2.10 
140 4.01 4.15 1.98 1.91 
141 4.17 4 2.00 1.97 

2 Xcontraction—CO2 


3 


separate determination of H2 and CO was made 
(e.g., in experiments 117, 112, and 118). How- 
ever, after it was determined that the fraction 
contained only these two gases and that the 
determination of Hz and CO by calculation from 
the contraction over anhydrone was reliable, the 
absorption of CO2 was omitted. Similarly, once 
it had been determined by analyses that the 
138°K fraction contained only CsHe, but few 
analyses on this fraction were conducted; usu- 
ally, the volume was simply measured and taken 
to be ethane. The data on the basis of which it 
was concluded that the 138°K fraction was solely 
ethane are summarized in Table II. 

The results obtained in the analyses in the 
various runs at wave-lengths of ~3200A, 
~2900A, and 2537A at temperatures of 25° and 
100° are summarized for experiments with lead 
mirrors in Table III and for experiments without 
mirrors in Table IV. Experiments are numbered 
in chronlogical sequence. 


Tests for Atomic Hydrogen at 3200A 


The guard mirror technique used by May, 
Taylor, and Burton* was emulated by us in an 
effort to detect atomic hydrogen in experiments 
at 3200A. In each experiment a heavy lead mirror 
4 cm wide was deposited at the edge of the 
irradiated zone. A faint, narrow standard anti- 
mony mirror was then deposited about 5 cm from 
the irradiated zone. The time of removal of the 
antimony mirror in the course of a usual experi- 
ment (photolysis of propionaldehyde streaming 
at an input pressure of 2.5 mm) was then noted. 
In three separate experiments it was found that 
a barely detectable fading of the antimony 
mirror, corresponding to no more than 10 percent 
removal, could be observed in 30 minutes. It had 
been previously established that at approximately 
the same free radical concentration (cf. Table 
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III) at 2537A similar antimony mirrors without 
the lead guard mirror were removed in duplicate 
experiments in less than two minutes. In other 
words, the free radicals present at 3200A were 
~ 150 times as effective in removing antimony as 
was the fraction which leaked by the guard; the 
concentration of atomic hydrogen seems to be 
even less than the maximum 2 percent reported 
by May, Taylor, and Burton.’ The reaction 


C,Hs;CHO-C;H;CO+H 


may therefore be neglected as a possible source 
of molecular hydrogen in these experiments. 


DISCUSSION 


Absorption of light by propionaldehyde in- 
volves only excitation of the carbonyl group.” 


(1) 


TABLE III. Summary of data with mirrors.¢ 


Time, H: co Ph CHs CoH, 


Expt. A,A min. cm cm cm 
119 3200 30° 3 0.42 7.95 35 0.31 4.71 
120 3 0.34 6.38 29 33 4.11 
128 3 0.53 8.90 39 31 6.07 
148° 1.5 043 6.92 36 3606 «(3.45 
149» 15 @42 731 33 5.04 
129 100° 3 0.62 9.25 37.5 288 6.42 
130 2900 30° 1.5 0.36 4.04 266 .454 2.25 
131 1.5 043 4.43 30.2 .477 2.69 
132 3 0.94 8.98 57.2 438 5.20 
150° 1.5 0.36 4.80 340 49 2.44 
151° 1.5 0.31 3.57 26.2 2.03 
133 100° 3 0.93 9.24 56.6 42 5.99 
134 3 0.90 9.10 53.5 40 5.40 
101 2537 30° 3 0.74 8.64 48 38 

102 3 0.71 8.27 45 37 

103 3 0.75 8.70 46 36 

104 3 0.77 8.83 

105 3 0.68 8.42 

117 0.57 6.15 3.97 
127 3 0.62 6.60 4.62 
142 3 0.71 7.94 43 313 542 
143 3 0.60 7.43 40 371 4.97 
106 100° 3 1.00 11.10 44 .274 

107 3 0.94. 10.46 

114 15 0.43 5.47 22.5 .284 

139 1.5 0.36 5.10 23.5 .318 3.41 
140 1.5 047 5.89 29.5 .345 4.01 


® Lead mirrors in each run were placed at the edge of the irradiated 
zone. The gas pressures are recorded in cm of mercury and were meas- 
ured at 23°C for a volume of 0.488 cc in the Saunders-Taylor apparatus. 
> 1 cm illuminated zone. 


12 Cf. G. K. Rollefson and M. Burton, Photochemistry and 
the Mechanism of Chemical Reactions (Prentice-Hall, Inc., 
New York, 1939), p. 223. 
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Thereafter, the energy may be made ineffective 
by fluorescence 


C,.H;CHO*—C,H;CHO+ hy, (2a) 


or by molecular deactivation 


(2b) 


or it may be made effective as the result of an 
energy shift within the molecule leading either to 
rearrangement and an ultimate molecule decom- 
position or to rupture into free radicals." 

The possible primary decomposition processes 
involved in this case are all predissociation 
mechanisms of different types. The results that 
may be expected depend in part upon the sta- 
bility (i.e., the ‘‘natural” life) of the excited 
state. In general, provided sufficient energy is 
available, a rupture into radicals seems favored 
by its simplicity. Whether the molecule decom- 
poses by rupture or by rearrangement will 
depend on (a) the amount of energy, (b) the 
time available, i.e., the life of the excited 
molecule, and (c) the extent of rearrangement 
required for the ultimate-molecule mechanism. 
Such considerations have led to conclusions that 
decomposition of complex molecules into radicals 
is favored by short wave-lengths and that 
ultimate-molecule decomposition, if it involves 
a profound rearrangement of the component 
parts, occurs with higher probability in a dis- 
torted field; i.e., as an induced predissociation."* 


The Free Radical Process 


Analyses of the mercury alkyl product show 
that the only alkyl radicals produced are C2Hs. 
Furthermore, the guard-mirror experiments of 
May, Taylor, and Burton* at 2537A and our own 
experiments at 3200A indicate the absence of H 
atoms in any significant quantity. Thus, the only 
free radical decomposition process we have to 
consider is 


C,H;CHO*—C;H,’+ HCO’. (3) 


The prime marks (’) after the radicals indicated 
that they are formed with excess energy depend- 


(1938); Burton and G, K. Rollefson, J. Chem. Phys. 6, 416 
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TABLE IV. Summary of data without mirrors. 


Time, co He CoHe 
Expt. A, A T min. cm cm cm 
124 3200 30° 8.27 0.50 7.01 
125 7.99 52 7.11 
144° 4.42 29 3.61 


1.5 
136 1.5 4.50 61 3.72 
146° 9 6.41 64 5.57 
147° 9 6.49 59 4.95 
137 100° 1.5 4.96 66 4.00 
138 1.5 4.75 69 3.90 
118 2537 30° 1.5 5.41 45 3.65 
121 1.5 6.01 46 3.94 
122 100° 3 6.79 54 4.70 
123 3 6.92 58 4.91 
141 3 6.70 52 4.17 


© Incident light at half intensity. 


ing on the wave-length." This energy is available, 
e.g., as activation energy, up to the first collision. 
Thereafter, the radicals reach thermal equi- 
librium very quickly. 


Radicals are produced at all wave-lengths— 
although not always to the same extent. The 
maximum yield of free radicals, in terms of the 
C.H;/CO ratio, occurs at intermediate wave- 
lengths (cf. Fig. 1), but this does not necessarily 
mean that the ratio of free radical to ultimate 
molecule processes reaches a maximum at such 
intermediate wave-lengths, nor even, indeed, 
that the wave-length 2900A has any special sig- 
nificance, for more thorough coverage of wave- 
lengths might relocate the maximum. 


Proposed Mechanism 


The data are not such as to warrant a com- 
pletely detailed consideration of the possible 
mechanism at all wave-lengths. However, certain 
aspects are sufficiently striking to merit special 
attention. The results at 2537A are susceptible 
to the simplest interpretation. Here a possible 
satisfactory mechanism can be contrived from 
reaction 3 followed or accompanied by the reac- 


14 The amount of excess energy is estimated on the basis 
of the value assumed for the strength of the C.H;CHO 
bond. If we assume ~75 kcal. [cf. M. Burton, J. Chem. 
Phys. 7, 1072 (1939) ], the excess energies are: at 2537A, 
37 kcal.; at 2900A, 22 kcal.; at 3200A, 15 kcal. 
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3200 2900 2537 


Wave-length, A 


_ Fic. 1. C:Hs/CO ratio as a function of wave-length. 
© Without mirror, 30°C, full intensity ; © Without mirror, 
30°C, half intensity; @ With mirror, 30°C, full intensity; 
e With mirror, 30°C, narrow illuminated zone; ® With 
mirror, 30°C, half intensity; o Without, mirror, 100°C, full 
intensity; m With mirror, 100°C, full intensity. 


tions shown. 


(4) 


2 Cth (5) 

2 CO. (6a) 
or ? 

2 HCO» H;CO+CO. (6b) 


Reactions 6a and 6b are intended to represent 
over-all reactions rather than precise mecha- 
nisms. It is possible that formaldehyde is formed 
by reaction 6b and that a portion of it then 
decomposes to give the over-all effect 6a. Our 
data cannot distinguish such reaction steps. 
There are also other possible reactions which 
require only casual mention, e.g., 


or ? 
(7b) 


If reactions 7a and 7b occur at all, their rates 
will be affected by the absence or presence of a 
mirror since the latter would remove the CH; 
radicals required. In turn, the participation of 
HCO radicals in reactions 6a and 6b would be 
affected. At 2537A, as can be seen from Figs. 2 
and 3, the C2Hs/He ratios remain constant,!® 


18 Tt should be noted that data involving Hz are more 
uncertain than others because of the greater relative ex- 
perimental error involved in the determination of the small 
volumes of hydrogen. 
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within experimental error, in the absence or the 
presence of a mirror. It is for this reason that we 
eliminate reactions 7a and 7b from further 
detailed consideration. 

Reaction 4 is an induced predissociation reac- 
tion of a type suggested by Taylor and Burton"* 
in the case of the nitric oxide-catalyzed pyrolysis 
of acetaldehyde; it is characterized by the fact 
that an odd-electron molecule furnishes the 
required intense distorted field. Evidence that 
reactions of this kind are not uncommon appears 
recently from several sources. For example, O. K. 


3200 2537 


Wave-length, A 


Fic. 2. C2Hs/CO ratio as a function of wave-length. 
For legend see Fig. 1. 


12.0 


10.0 


2900 2537 


Wave-length, A 


Fic. 3. C2He/H: ratio as a function of wave-length. 
For legend see Fig. 1. 


3200 


16H. A. Taylor and M. Burton, J. Chem. Phys. 7, 414 
(1939). 
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PHOTOLYSIS OF PROPIONALDEHYDE 


Rice and Hayden! have stated that the rela- 
tively stable HCO radical may be an even more 
effective chain carrier than CH; in the high 
temperature photolysis of acetaldehyde. Also, 
although Burton, Ricci, and Davis'* found that 
free radicals are produced in the pyrolysis of 
acetaldehyde even at the lowest temperatures 
studied, Morris'® obtained the apparently con- 
tradictory result that mixed deutero-methanes 
were not found as products in the pyrolysis of 
mixtures of normal- and deutero-acetaldehyde. 
The contradiction follows from the assumption 
of a Rice-Herzfeld®® chain-continuing step, e.g., 


CH3+ CD;CDO—CH;D +CD;CO. 


However, if the chain-continuing step is an 
induced predissociation, as in reaction 4, there is 
no contradiction in the two sets of results. 

Induced predissociation of this type requires 
an activation energy which is ordinarily available 
only at high temperature. However, the radicals, 
when produced, possess considerable excess 
energy (i.e., they are quite ‘‘hot’’) and retain 
this energy until their first collision. Thus, if the 
first collision is with a molecule like C.H;CHO 
in this case, there may be enough energy for the 
induced predissociation. Recently O. K. Rice and 
Hayden” have estimated ~11 kcal. for the 
energy of the chain-continuing step in acetal- 
dehyde, a figure which indicates that the excess 
energy available in R’ in this case is of the order 
of magnitude required to produce the effects 
noted, particularly at shorter wave-lengths. For 
this reason, the apparent ‘‘chain length”’ of this 
step (which can in no event exceed unity accord- 
ing to the explanation given) must decrease as 
the wave-length is increased. We shall see that 
it is this theoretical requirement which suggests 
the necessity for an ultimate molecule step, at 
least at longer wave-length, 


C,H;CHO*—C;He+ CO. (1 1) 


Reaction 5 is offered without further discussion 
to account for the butane reported by Leighton 


(1943), K. Rice and W. L. Hayden, J. Chem. Phys. 10, 445 
se Burton, Ricci, and Davis, J. Am. Chem. Soc. 62, 265 
1 J. Morris, J. Am. Chem. Soc. 63, 2535 (1941). 
20F.O. Rice and K. F. Herzfeld, J. Am. Chem. Soc. 56, 
284 (1934). 
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TABLE V. Calculated and observed C2H;/CO ratios 
at 2537A. 


C2H3/CO 


Discrepancy 
% calc. 


—0.20 
— .26 


139 — 33 


140 


Average 


and Blacet.! We have not been concerned with 
this feature of the yield in our work. 

Reactions 6a and 6b are adopted as the only 
reactions involving two formyl radicals in this 
case. Frankenburger and his co-workers,”! work- 
ing with the mercury-photosensitized reaction 
between hydrogen and carbon monoxide (at 
pressures considerably higher than employed by 
us), found formaldehyde and glyoxal, formed 
presumably by reactions involving HCO, in the 
ratio of 4 : 1. On the other hand, Zemany” has 
adduced evidence that in the photolysis of 
glyoxal at pressures of a few millimeters the reac- 
tion between HCO radicals yields principally 
CO, He, and HCHO and a negligible quantity of 
(HCO).. Since our conditions are more nearly 
similar to those of Zemany and, in any event, 
the neglect of the possible formation of glyoxal 
does not introduce any serious error into our 
considerations, we make the simplifying assump- 
tion that the only significant reactions involving 
two free HCO radicals are 6a and 6b. The internal 
consistency of this assumption in our work can 
be seen by an arithmetical examination of the 
data at 2537A. 

If we limit ourselves only to the reactions so 
far considered, we see that possible reactions at 
2537A producing our measured products are: 
3 producing R’; 4 and 11 producing C2H¢ and 
CO; 6a producing Hz and CO; and 6b producing 
CO. Let us represent the total cm of C2H¢ by x, 
CO by y, and He by z.% The number of cm 
produced in any particular reaction or group of 


2 Frankenburger, Klinkhardt, Steigenwald, and Zimmer- 
man, Zeits. f. Elektrochemie 36, 757 (1930). 

2 P. Zemany and M. Burton, unpublished work. 

23 Moles of a gas are proportional to the pressure regis- 
tered by a fixed volume in the Saunders-Taylor apparatus. 
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reactions will be indicated by the appropriate 
subscript. Thus, 

N41 = 
and 

Veo = —X—2z2. 2 


The number of cm decomposed to give free radi- 
cals is 2(y —x —z) and consequently the calculated 
C.H;/CO ratio (fear) should be 2(y—x—z)/y. 
The results at 2537A are compared with the 
observed ratios = C2H;/CO in Table ITI) in 
Table V. In the last column the discrepancy is 
noted as a fraction of the calculated value. This 
fraction should be equal to the fraction of radicals 
which have disappeared before reaching the 
mirror. This latter may be calculated from our 
knowledge of the half-life of the C2.H; radicals 
(~4.0X10-* sec.), the length of the irradiated 
zone (10 cm), the distance of the mirror from the 
edge of the irradiated zone (0 cm), and the 
streaming velocity (1120 cm/sec.). The result is 
31 percent, a value which is in better agreement 
with the average discrepancy shown in Table V 
than can reasonably be expected from the 
accuracy of the analyses (+5 percent). 


Other Reactions 


It must be emphasized that the data so far 
considered furnish no criterion for a choice 
between reactions 4 and 11. Either one or both 
may be the source of C2H¢. Further conclusions 
necessitate examination of the data at longer 
wave-lengths, the striking features of which are 
best illustrated in Figs. 1 to 3. For example, the 
C:H;/CO ratio (Fig. 1) reaches a maximum at 
some intermediate wave-length. If reaction 4 is 
the principal source of C:He, the result is easily 
accounted for; we conclude merely that at longer 
wave-length reaction 4 is less frequent since less 
energy is available in R’. The existence of a 
maximum at intermediate wave-lengths can 
then be interpreted on the basis of a precon- 
ceived notion that the more complicated process 
of rearrangement involved in reaction 11 takes 
place with greater probability at longer wave- 
length.'* This conclusion is supported also by the 
CO/H: and C2H¢/He ratios which increase 
markedly at shorter wave-lengths, indicating the 
incidence of a mechanism which does not 
yield Hoe. 
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This conclusion of itself is not sufficient to 
account for the other results. The C2:H»/CO 
ratios at 2900A and 3200A are clearly sensitive 
to the presence or the absence of a mirror (cf. 
Fig. 2). Such a result indicates that C.H; 
radicals at these wave-lengths (but not at 2537A) 
react with some substance, most probably 
C.H;CHO, to yield C2H¢ without yielding CO. 
Obviously, the reaction cannot involve the 
(“‘hot’’) C.H;’ radicals since they do not survive 
a first collision ; i.e., 


C.H;’+M-C:H;+M or products. 
Neither can it involve a reaction such as 


residue, (12) 


where the assumption is that the activation 
energy is very low, for on such an assumption 
C.H./CO should increase also at 2537A. 

Possible acceptable explanations of the results 
may involve reactions such as 


residue, 
HCO+C.H;CHO*—CO + residue, 


(8) 
(9) 


Residues—Polymers, (10) 


and, in addition to those, 


(13) 
The latter reaction for example, might account 
for the increases in CO/He and C2H¢/H: already 
noted at 3200A. However, a detailed considera- 
tion of any of these reactions is certainly not 
justified on the basis of the data at hand. The 
statement should, however, be made that if 
excited molecules such as C:H;CHO* are in- 
volved in reactions such as 8 and 9 (where 8 
accounts for the C.H¢/CO effects already noted), 
then the half-life implicitly required for the ex- 
cited state is >10-* sec., as judged from the time 
required for the radicals to pass from the irradi- 
ated zone into the remoter portions of the reaction 
tube.” 


*4 Note the approximate constancy of CO/H: ratios indi- 
cated at 2537 and 2900A. 

26 Cf. R. Spence and W. Wild, J. Chem. Soc. 590 (1941), 
who have recently arrived at a similar conclusion regarding 
the excited state of acetone. 
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Conclusion on Mechanism of Photolysis 


The major conclusion from the data here 
presented is that a free radical decomposition 
occurs at all wave-lengths. An ultimate molecule 
decomposition may also occur but, if it does, its 
probability is greatest at the longer wave- 
lengths. On the other hand, it is possible to 
devise an exclusively free radical process (involv- 
ing a reaction such as 13) which may account for 
the results at longer wave-lengths. The data are 
not sufficiently extensive to warrant a critical 
evaluation of the possible mechanisms. 


Mechanism of Energy Transfer 


The simplest possible mechanism after reac- 
tion 1 would seem to be a spontaneous predis- 
sociation transition starting from the excited 
state represented by C.H;CHO* and ending with 
any repulsive or weakly attractive state whose 
dissociation level is less than 90 to 112 kcal. 
above the ground state. According to the bond 
strengths we have assumed,?® H atoms could 
then appear at 2900 and 2537A and free CH; 
radicals might be expected at all wave-lengths. 
The failure to observe the latter at all indicates 
that geometric factors are important. 

According to Mulliken,?’ the excited electron 
is a 2p non-bonding one on the oxygen atom. 
The experimental results here reported would 
indicate that although this excited state can 
interact with a state involving excitation of an 
adjacent bond it cannot do so with a state involv- 
ing excitation at a more remote bond. The failure 
to observe atomic hydrogen in any substantial 
amount must then be attributed to another 
factor. 

The conditions bearing on the decomposition 
into CsHs and CO have already been discussed 
for the parallel case of acetaldehyde.” In regard 
to the alternative decompositions yielding H or 
C>H; in the primary step, it is not unreasonable 


6 M. Burton, J. Chem. Phys. 6, 818 (1938). See, however, 
Andersen, Kistiakowsky, and Van Artsdalen, ibid. 10, 305 
(1942); also D. P. Stevenson, ibid. 10, 291 (1942). 

27 R. S. Mulliken, J. Chem. Phys. 3, 564 (1935). 

(1939), W. Davis and M. Burton, J. Chem. Phys. 7, 1075 
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to assume that the ~20 kcal. difference in the 
dissociation levels may be reflected by a large 
difference in the heights of the intersections of 
the repulsive (or perhaps weakly attractive) 
states representing dissociation into H and C2Hs, 
respectively, with an excited state of propion- 
aldehyde (a state, incidentally, which may be 
different from that to which the original excita- 
tion occurred). Such a difference might markedly 
favor the predissociation process yielding C2H; 
radicals?® but it is impossible to arrive at any 
better interpretation of the results without a 
more precise knowledge of the potential energy 
relationships. For example, conditions have been 
indicated where (unlike this case) decomposition 
to products on the higher dissociation level may 
be favored.?° 

That ~20 kcal. may be more than sufficient 
to inhibit a decomposition even when there is 
sufficient energy available is well illustrated by 
the case of crotonaldehyde. In that case, the 
conjugated double bond may strengthen the 
R—CHO bond by perhaps 10 kcal.*® No photo- 
decomposition can be detected at room tem- 
perature at wave-lengths from 3130 to 2380A.* 
Neither the R—CHO bond nor the remote 
CH;—C bond is broken even though more than 
enough energy is available for either process. 
The fact that the CH;—C bond in this case is 
not broken indicates that in the case of a “‘long- 
chain” aldehyde we are not dealing with possible 
competition between two or more possible 
mechanisms (i.e., breaks at various C —C bonds) 
but with only one possible rupture process," the 
break of the C—C link adjacent to the excited 
atom. 
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(1939), K. Rollefson and M. Burton, J. Chem. Phys. 6, 674 

30 Cf. the case of glyoxal, reference 26. 
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The experimental heat capacity data for crystalline Cle, CO2, SO2, SCO, N2O, C2Ne, CeHe, 
C:H,, and CH;Br have been compared with the results of calculations made employing a semi- 
theoretical method suggested by Lord, Ahlberg, and Andrews and amplified by Lord. Good 
agreement between the calculated and experimental results is obtained. 


HE purpose of this paper is to examine the 

experimental heat capacity data for a 
number of simple molecular solids in terms of the 
theoretical views advanced by Andrews'* and 
Lord, Ahlberg, and Andrews!’ and more recently 
by Lord.? According to this point of view the 
heat capacity of a molecular solid at constant 
pressure C, is representable by an equation of 
the form 


Cr+ La 1) PT, (1) 


with C,=Cy1)+C.u, a and b being constants 
assumed independent of 7. Here C,,1) represents 
the contribution of the translational and rota- 
tional modes of vibration of the molecules, 
acting as units, to the heat capacity of the solid 
at constant volume. This is given by Cy 1 
=nD(6/T), where n is the number of rotational 
and translation degrees of freedom possessed by 
the molecule in the gas state, and D(6@/T) is the 
usual Debye heat capacity function. The contri- 
bution of the internal vibration of the molecules 
to the heat capacity (at constant volume) is 
given by C,g)= >; E(0:;/T), where 7 represents 
the number of such internal degrees of freedom 
per molecule and E(6/T) the Einstein heat 
capacity function, 6; being calculated from the 
spectroscopically determined internal frequencies 
of the molecule in the gas state. 

In the absence of theoretical means of deter- 
mining 0, a, and b at the present time we have 

* Present address: U. S. Rubber Company, Detroit, 
Michigan. 

¢ Contribution from the Cryogenic Laboratory, Depart- 
ment of Chemistry, The Johns Hopkins University. 

1¢ D. H. Andrews, Proc. Roy. Acad. Amsterdam 29, 744 
(1926) ; * Lord, Ahlberg, and Andrews, J. Chem. Phys. 5, 
649 (1937). 


*R. C. Lord, J. Chem. Phys. 9, 693 (1941). Paper I of 
this series. 


contented ourselves with an examination of the 
validity of Eq. (1) by making an analysis of the 
experimental heat capacity data for a number of 
simple molecular solids in terms of the above- 
mentioned point of view. The molecular solids 
chosen for analysis were Cle, CO2, SOx, SCO, 
N2O, C.No, CeHe, and CH3;Br. These 
molecules have been chosen because the requisite 
vibrational analyses have been carried out for 
the gas state and also precise studies of the heat 
capacities of the solid state have been made. 


METHOD OF COMPARISON 
A. Chlorine 


As an illustration of the technique employed 
in obtaining values for 6, a, and b in Eq. (1) from 
the experimental heat capacity data and known 
internal frequencies of the molecule we shall 
consider solid chlorine in some detail. 

Clearly there are three desiderata: (1) a 
temperature range in which C,—C, and C,, are 
both very small or zero; (2) a temperature range 
in which C,—C, and C,,1) are large while C, 1 is 
small; (3) a temperature range in which C,—C, 
and C1) are large while C,,z) has very nearly 
reached a value approaching its upper limit. 
The first temperature range permits the deter- 
mination of 9, the second a, and the third b. In 
any actual case these conditions will not be 
entirely fulfilled and it will usually be necessary 
to find an approximate value of a in order to 
obtain a value for 6.8 

The heat capacity of solid chlorine has been 


3Such an approximate value may be found by using 
the Nernst-Lindemann melting point equation C,—C, 
ae T (m.p.). Zeits. f. Elektrochemie 17, 817 
1911). 
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measured by Giauque and Powell’ with good 
accuracy from 15°K to almost the melting 
point (172°K), and the frequency of vibration 
of the molecule in the gas state is also known.°® 

Since C,«7 for Cle below 60°K is 0.001 cal./ 
mole/deg. or less we may use this range to 
evaluate 6. The Lindemann melting point 
formula® gives an approximate value for a=1.1 
X10 cal.-?. Employing this value of a to 
obtain (C,—C,), we have calculated 6p (using 
n=5) using the experimental values of C, 
(interpolated from a smooth curve) in the range 
15-50°K. The resulting 6 values show a tempera- 
ture trend increasing from approximately 141 at 
15° to 146 at 50°. In terms of the analysis made 
in a previous paper* this observed upward trend 
of @ with increasing temperature is to be inter- 
preted as meaning that the ratio of the Debye 
6’s for rotation and translation is not unity but 
nearer to 1.5, and that the apparent 6 will reach 
a limiting value of about 150. This value (@= 150) 
was then chosen for calculating C,:z) throughout 
the entire temperature range (15-170°K) in 
Eq. (1). 

The more exact calculation of a may then be 
carried out by choosing a temperature such that 
the term aC, 1) is large compared to the term 
This usually requires that C,<7) be small 
compared to Cz) since in practice we find a 
and b to be of the same order of magnitude. 
Choosing T=100°K with C, (observed) = 10.10 


TABLE I. 


Coty Cen Cp—Cr Cp (cale.) Cp (obs.) 


SL 


SSO 


The numerical results listed in Tables I and III are expressed in 
calories. One calorie (15°) =4.183 international joules. The gas constant 
R was taken to be 1.986 cal. deg.~! mole. 


*W. F. Giauque and T. M. Powell, J. Am. Chem. Soc. 
61, 1970 (1939). 
5A. Elliott, Proc. Roy. Soc. London A127, 638 (1930). 
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cal. deg.-! mole, C,=C.¢1) + Cou) = 8.92 +0.04 
we obtain a=1.2X10~ cal.~?. 

Having obtained a value for a we may evalu- 
ate the constant 6. Here we desire C,,;) to be as 
large as possible. Selecting T=160°, where 
Cour) =0.315 cal. deg.-! mole~!, we solve Eq. (1) 
directly to obtain b=5.9X10~ cal.-?. 

These values of 6, a, and 6 have been used to 
calculate C, throughout the temperature range 
15-170°K. Table I summarizes the calculation. 
The next to last column shows A= C, (observed) 
—C, (calc.) for various temperatures. The agree- 
ment, on the whole, is quite good. The chief 
discrepancy lies in the temperature range 
15-50°K. This is due, of course, to the choice of 
6=150 as a value to be used at all temperatures. 
Clearly one may improve the general fit of the 
experimental data by using a somewhat lower 
value for 6, say 6=148. Since has reached 
85 percent of its maximum value at 80°K such 
a change in @ will leave a and 6 virtually un- 
changed while improving the fit at low tempera- 
tures. This is shown in the last column of Table I 
where A*=C, (obs.) —C, (calc.) and the param- 
eters 0=148°, a=1.2X10~ cal.-?, and b=6 
X10-? cal.-? have been used. 


TABLE IT. 


Source of data 


Heat Frequency 


Molecule capacity assignment 


& 


1 Giauque and Egan, J. Chem. Phys. 5, 45 (1937). 

2 Giauque and Stephenson, J. Am. Chem. Soc. 60, 1389 (1938). 

3 Kemp and Giauque, J. Am. Chem. Soc. 59, 79 (1937). 

4 Blue and Giauque, J. Am. Chem. Soc. 57, 991 (1935). 

5 Ahlberg, Blanchard, and Lundberg, J. Chem. Phys. 5, 539 (1937). 

6 Huffman, Parks, and Daniels, J. Am. Chem. Soc. 52, 1547 (1930). 

7 Nernst, Ann. d. Physik 36, 395 (1911). 

8 Ruehrwein and Giauque, J. Am. Chem. Soc. 61, 2940 — 

®’ Brown and Manov, J. Am. Chem. Soc. 59, 500 500 (1 937 

10 Egan and Kemp, J. Am. Chem. Soc. 59, 1264 (1937). 

11 Egan and Kemp, J. Am. Chem. Soc. 60, 2097 (1938). 

12 Adel and Dennison, Phys. Rev. 44, 99 fg Hl ibid. 43, 716 (1933). 

13 Badger and Bonner, Phys. Rev. 43, 305 (1933 

14 Bartunek and Barker, Phys. Rev. 48, 516 (193 5). 

15 Plyler and Barker, Phys. Rev. 38, 1827 (1931); 41, 369 (1932). 

16 Langseth and Lord, Kgl. Danske Vid. Sels. Math-F ys. Medd. 16, 
75, No. 6 (1938), Table 18. The missing frequencies in*this Table, Nos. 
14 and 15, make only a trifling contribution to C,. They have been 
estimated to be 1600 and 1160 cm™, respectively. 

17 Woo and Badger, Phys. Rev. 39, 932 (1932). See also Burcik and 
Yost, J. Chem. Phys. 7, 1114 (1939); Stitt, ibid. 1115 (1939). 

18 Dennison and Wright, Phys. Rev. 38, 2077 (1931). 

19 Bonner, J. Am. Chem. Soc. 58, 34 (1936). 

20 The frequencies used by Egan and Kemp (reference 11) to com- 
pare the thermal and spectroscopic entropies were used. See also 
Stevenson and Beach, J. Chem. Phys. 6, 25, 108 (1938). 


10 
the 
the 
of 
ve- 
ids 
‘O, 
ese 
site 
for 
eat 
ved 
om a b 
wn 
CO2 171 0.0110 0.030 (1) (12) 

all SO: 172 0.0098 (13) 

sco 131 0.0122 0.020 (3) (14) 

N:O 160 0.0118 0.034 (4) (15) 

are K CoH. 170 0.0141 (10) (19) 
1ge 15 0.00 40.14 +0.11 150 0.0095 0.01 (11) (20) 

30 06 3.72- 25 17 
-C, 40 ‘17 5.54 19 M1 

50 ‘31 6.89 10 02 
rly 60 ‘477.89 05 
70 002 ‘64 8.63 05 ‘00 
nit. 80 007 82 9.21 02 —0.02 

00 019 ‘o1 ‘02 
100 037 19 10.13 1 —0.03 06 

110 (064 39 10.47 ‘05 08 
In 120 10 54 1083 10.87 ‘04 ‘08 

130 14 87 11.31 11.29 02 05 
be 140 20 1S 11.73 ‘00 ‘03 

150 25 48 12.18 12.20 +40.02 00 
ary 160 32 81 12.63 12.68 05 +0.02 

170 38 19 13.12 13.17 05 ‘02 
to | 
een 
sing 
-C, 
817 


W. F. BRUCKSCH, JR. 


AND W. T. ZIEGLER 


TABLE III. 


CO2 


Cp(obs) ©p(calc) Cp(obs) 


Cp(obs) 


0.54 
1.23 
2.14 
3.08 
4.69 
6.10 


Cpr(cale) 4 


CeHe 


Cp(obs) Cp(calc) 4 Cp(calc) 


Cp(calc) 4 


0.75 +0.11 
1.62 ‘ 
2.65 


0.91 +0.01 


0.63 +0.05 
.10 


14.00.12 


For the molecular solids studied in this paper 
we have always adjusted the @ to give a good 
fit at low temperatures in the manner described 
above. 

B. Other Molecules 


The other molecular solids whose heat capa- 
cities have been analyzed in the above manner 


are listed in Table II together with the source of 
the heat capacity data and frequency assign- 
ments employed in making the comparisons. In 
addition, Table II presents the values used for 
(the number of translational and rotational 
degrees of freedom of the molecule) 6, a, and 
for each solid. These have been used to make the 
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|| SO2 sco N:O 
°K Cp(calc) 4 Cp(obs) Cp(calc) 4 p\calc) 
15 0.52 +0.02 0.83 0.61 +0.22 1.37 1.10 +0.27 0.67 0.63 +0.04 
20 1.16 07 1.68 1.37 31 2.61 ~—-2..21 40 1.50 1.38 12 
25 2.01 13 2.72 ~=2.39 33 3.78 3.43 35 2.46 2.34 12 
30 2.95 13 3.79 3.49 30 4.84 4.56 28 3.45 3.33 12 
40 4.70  —0.01 5.78 5.59 19 6.62 6.35 27 5.13 5.13 .00 
50 6.11 01 7.36 7.28 08 7.70 7.61 09 6.52 6.52 .00 
60 7.18 7.17 +0.01 8.62 8.55 07 848 849 —0.01 7.56 7.57. —0.01 
70 7.97 7.98 —0.01 9.57 9.52 05 9.08 9.14 06 8.32 8.35 03 
80 8.58 8.61 03 | 10.32 10.29 03 9.62 9.66 04 8.95 8.97 02 
90 911 9.12 01 | 10.93 10.91 02 | 10.09 10.12 03 9.44 9.47 03 
100 9.53 9.54 O01 | 11.49 11.45 04 | 10.53 10.54 01 9.90 9.94 04 
110 9.92 9.93 01 | 11.97 11.95 02 | 10.95 10.96 01 | 10.32 10.37 05 
. 120 | 10.30 10.28 +0.02 | 12.40 1243 -0.03 | 11.40 11.39 +0.01 | 10.77 10.79 02 
130 | 10.67 10.62 05 | 12.83 12.88 05 | 11.85 11.85 00 | 11.25 11.24 +40.01 
. 140 | 11.04 10.99 05 | 13.31 13.36 05 11.72 11.70 02 
150 | 11.39 11.36 03 | 13.82 13.84 02 12.19 12.20 —0.01 
160 | 11.77 11.75 02 | 14.33 14.33 00 12.71. 12.73 02 
170 | 12.17 12.17 00 | 14.85 14.85 .00 13.30 13.28  +40.02 
180 | 12.61 12.60 1 | 15.42 15.38 +0.04 13.98 13.87 Al 
190 | 13.07 13.05 02 | 16.02 15.93 .09 
C2N2 CS2 CHsBr 
°K {Co(obs) Cpicatc) 4 ([Cp(obs) 3 
15 | 0.86 1.64 1.24 +040) 0.92 0.68 1.54 0.90 +0.64 
20 | 1.76 2.82 244 .38| 1.84 1.94 —0.10] 1.52 2.56 1.94 .62 | 
25 | 2.82 3.92 3.71 .21| 3.00 3.20 .20| 2.54 246 .08| 3.60 3.18 42 
30 | 3.78 3.70 .08| 4.97 484 13] 4.24 446 .22| 3.56 3.62 —0.06| 4.54 4.44 .10 
40 | 5.57 5.52 .05| 665 659 .06| 647 6.58 .11| 5.67 5.81 .14| 6.37 6.55 —0.18 
50 | 7.02 694 .08| 7.83 7.79 .04| 814 815 7.39 7.63 7.93 818  .25 
60 | 820 8.11 .09| 862 9.32 9.28 +0.04/ 8.92 9.05 9.16 9.33 .17 
70 | 9.13 9.08 .05| 9.28 9.26 .02]10.16 10.20 —0.04| 10.18 10.18  .00}10.02 10.14  .12 
80 | 9.97 9.94 .03| 9.89 9.83  .06/10.85 10.90 .05| 11.55 11.11 +0.44| 10.79 10.79 .00 
90 | 10.71 10.73 —0.02| 10.49 10.35 .14/11.44 11.53 13.29 11.84 1.45|11.42 11.31 +0.11 
100 | 11.41 11.45 .04/11.07 10.86 .21/11.98 12.16  .18/17.05 12.52 4.53}11.92 11.73 .19 
110 |12.08 12.12 .04] 11.61 11.37. .24/12.61 12.77 .16 12.34 12.09 .27 
120 | 12.73 12.77 .04/12.09 11.89 .20/13.32 13.36 .04 12.67 12.43 .24 
130 | 13.35 13.39 .04/12.52 12.42 .10| 14.02 14.03 .01 12.96 12.76 .20 
140 | 13.91 13.99  .08|12.92 12.94 —0.02| 14.76 14.73 +0.03 13.21 13.05 .16 
150 |14.47 14.58 .11/13.28 13.48 .20/15.53 15.48  .05 13.44 13.37 07 
160 |15.02 15.14 .12}13.63 14.01 .38| 16.31 16.28  .03 13.66 13.69 —0.03 
170 |15.58 15.69 17.10 17.14 —0.04 1338 
180 | 16.20 16.23 .03 18.01 18.05 .04 
190 | 16.82 16.77 +0.05 18.99 19.05 .06 
200 |17.45 17.29 20.07 20.10 .03 
210 | 18.11 17.79  .32 21.21 21.25 .04 
220 |18.82 18.29  .53 22.30 22.38 .08 
230 | 19.54 18.78  .76 23.52 23.56 .04 
240 | 20.27 19.27 1.00 24.80 24.80 .00 
250 26.13 26.10 +0.03 
260 27.75 27.44 31 
270 29.55 28.83 .72 


HEAT CAPACITIES OF MOLECULAR LATTICES 


calculations of C, listed in Table III. In all 
tables the experimental C, values listed have 
been obtained from a smooth curve through the 
data given in the reference cited. 


DISCUSSION OF RESULTS 


An examination of Tables I and III shows that 
Eq. (1) does indeed represent the experimental 
heat capacity data with good accuracy, especially 
for the molecules Clz, COs, N2O, and CeHe, 
where the disagreement is of the same order of 
magnitude as the experimental errors in C,. In 
all the molecular solids studied the agreement is 
poorest for very low temperatures (below 30°K). 
As has already been pointed out this disagree- 
ment is to be expected, since, in general, the 
arbitrary choice of a single Debye function to 
represent both the translational and torsional 
motions of the molecules is certainly inadmissible 
in this temperature range. That the agreement 
can be quite good to temperatures as low as 4°K 
has been shown by Lord, Ahlberg, and Andrews! 

for benzene. 

' It should be pointed out that the good 
agreement between the calculated and experi- 
mental values of C, for these molecular solids is 
very dependent upon the correct choice of n 
(Table II). Thus, good agreement between the 
calculated and experimental values of C, for 
crystalline SCO is obtained only if the molecule 
is assumed to be linear (n=5). This choice is in 
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agreement with the conclusion of Eucken and 
Schaefer® that an angular model’ for the SCO 
molecule is incompatible with ultrasonic meas- 
urements of the heat capacity of the gas. 

The discrepancies which arise for C2N2 and 
C2H, in the neighborhood of the melting point 
are probably to be attributed to the fact that as 
melting becomes imminent the motion of the 
molecules comprising the solid is no longer 
harmonic, their amplitudes becoming very large. 
An analysis of the heat capacity data for 
crystalline SiH,y, SFs, and CCl, exhibited 
similar disagreements in the neighborhood of the 
melting or transition temperatures. 

From Table II it may be seen that the values 
of the constants a and b are of the same order of 
magnitude. Due to the fact that the solid usually 
melts before the contribution of the internal 
molecular vibrations to C, becomes large the 
constant 0 is less accurately determined than a. 
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as 2 Eucken and K. Schaefer, Chem. Abst. 36, 2767 

7J. Wagner, Zeits. f. physik. Chemie B48, 309 (1941). 

8 A. Eucken, Zeits. f. Elektrochemie 45, 126 (1939), has 
carried out a somewhat different analysis of the experi- 
mental C, data and has concluded that the rapid increase 
of C, t the higher temperatures is due to the character of 
the torsional motion of the molecules in the lattice. 
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Erratum; Contact Angles at Liquid-Liquid-Air Interfaces 
[J. Chem. Phys. 10, 623 (1942) ] 


Fox 
Department of Chemistry, Columbia University, New York, New York 
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Some Factors Affecting the Reduction of Silver Oxalate* 


T. H. JAMES 
Kodak Research Laboratories, Eastman Kodak Company, Rochester, New York 


(Received September 8, 1942) 


(1) Silver oxalate precipitates age rapidly, and the process is not greatly impeded by ad- 
sorption of wool violet or 3,3’-diethyl-9-methylthiacarbocyanine. (2) Gelatin markedly de- 
creases the rate of reaction between silver oxalate and hydroxylamine in slightly acid solution. 
A break is obtained in the rate-gelatin curve corresponding to a definite adsorption layer. The 
thickness of this layer is the same for silver oxalate, silver chloride, and silver thiocyanate. 
(3) The dye, 3,3’-diethyl-9-methylthiacarbocyanine, has a much smaller effect than gelatin 
upon the reaction rate. (4) The effect of bromide ion upon the gaseous products of the reaction 
between silver oxalate and hydroxylamine in alkaline solution is considerably smaller than in 
the corresponding reactions of silver chloride and thiocyanate. 


HE thermal decomposition of silver oxalate 

has been the subject of several investiga- 

tions.’ This reaction is markedly autocatalytic 

in nature, the silver acting as the catalyst. 

Nuclei of silver sulfide also effectively catalyze 
the decomposition.® 

It is possible to prepare photographic emul- 
sions of silver oxalate which can be “developed” 
by heat alone. However, it is also possible to 
develop silver oxalate emulsions with normal de- 
veloping agents in solutions of rather weak ac- 
tivity. Arens and Eggert’ found that silver 
oxalate-gelatin plates resemble silver chloride 
emulsions in sensitivity and spectral range, but 
knowledge of the reaction between silver oxalate 
and developing agents is scant. 

Some points of interest have been found in the 
reaction of silver oxalate with hydroxylamine. 
The oxalate represents a different type of salt 
from the chloride® and thiocyanate® previously 
studied, and differences in surface conditions and 


* Communication No. 882 from the Kodak Research 
Laboratories. 
_ 1J. Y. Macdonald and C. N. Hinshelwood, J. Chem. Soc. 
127, 2764 (1925). 

2 A.F. Benton and G. L. Cunningham, J. Am. Chem. Soc. 
57, 2227 (1935). 

3 J. Y. Macdonald, Nature 137, 152 (1936) ; J. Chem. Soc. 
832 (1936). 

4 J. Y. Macdonald and R. Sandison, Trans. Faraday Soc. 
34, 589 (1938). 

5S. E. Sheppard and W. Vanselow, J. Am. Chem. Soc. 
52, 3468 (1930). 

6S. E. Sheppard and W. Vanselow, U. S. Patents No. 
= (October 9, 1934); No. 2,129,242 (September 6, 
(1931) Arens and J. Eggert, Phot. Korr. 67, No. 8, 17 
8 T. H. James, J. Am. Chem. Soc. 62, 536, 1649 (1940). 
®°T. H. James, J. Chem. Phys. 10, 646 (1942). 


adsorption are to be expected. Oxalate has no 
tendency to form complexes with silver ion’® so 
that no solvent effects are present due to excess 
oxalate ion. 

The precipitates employed in the present work 
were prepared by the addition of 5.00 ml of 
0.10M potassium oxalate to 15 ml 0.0333M 
silver nitrate. The procedure was similar to that 
followed in the preparation of the chloride and 
thiocyanate precipitates.*® Reaction with hy- 
droxylamine was carried out at 20°C in a solution 
buffered with sodium acetate-acetic acid mix- 
tures. 

The silver oxalate precipitates age rapidly and, 
unlike silver bromide or chloride, the aging 
process is not seriously impeded by adsorption of 


TABLE I. Variation of rate with age of precipitate. 


10H. Schafer and R. Abegg, Zeits. f. anorg. allgem. 
Chemie 45, 293 (1905). 
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A 
pH R i/t 
0.1 0.40 0.26 j ) 
| 41 | J 
1.0 11 .10 
1.0 21 
1.1 19 14 
15 Az 
4.0 .23 
17 .23 19 
19 .26 .20 
20 .28 .22 
42.5 .28 .23 
1.0 .044 .049 
.045 
; 1.2 .042 .044 
47 .050 .050 
65 .047 .053 
70 .050 .055 


REDUCTION 


OF SILVER OXALATE 


Fic. 1. Reduction of AgeC.0, by NH2OH. Effect of gelatin and 
dye at pH 5.27. -o—o-, R; -x—x-, 1/t. 


wool violet or the thiacarbocyanine dyes. Small 
amounts of 3,3’-diethyl-9-methylthiacarbocya- 
nine (ca. 0.05 mg) were completely adsorbed 
from solution by freshly prepared oxalate precipi- 


tates, but part of this dye was desorbed within a 
few minutes. The dye remaining adsorbed to the 
oxalate after ten minutes could be extracted with 
water. 

Larger quantities of wool violet were adsorbed 
by the fresh precipitate, but again desorption 
occurred within a few minutes. The specific 
surface of a precipitate which had aged for 1 
hor was near the minimum value. Three precipi- 

(fates aged for 1 to 2 hours adsorbed 0.045 +0.003 
ging dye; a precipitate aged 66 hours adsorbed 
004 mg of dye. 
he reduction of silver oxalate by hydroxyl- 
amine is markedly accelerated by the silver 
formed in the reaction. In neutral or acid solu- 
tion, nitrogen is almost exclusively the oxidation 
product of the hydroxylamine, indicating true 


TABLE II. Specific surface and adsorption of gelatin. 


Relative Relative 
specific average 
surface gelatin 


W. V. Iv* Gelatin 
Salt mg mg mg 


AgeC:0, 0.043 0.15-0.20 
A 135 0.105 - 6 
45 2.0 -3.0 


* 3,3’-diethyl-9-methylthiacarbocyanine chloride. 


silver catalysis." The variation of reaction rate 
with the age of the precipitate is given in Table I, 
in which R is the slope of the reaction curve when 
reaction is 10 to 20 percent complete (approxi- 
mately the maximum rate), and ¢ is the time 
required to complete 3 percent of the reaction. 

A definite drop in reaction rate is observed 
between 0.1 hour and 1.0 hour. Further aging of 
the precipitate, however, produces a small but 
definite increase in reaction rate.* It is improbable 
that thermal decomposition of the oxalate during 
the aging process could account for the observed 
increase in reduction rate. The increase is prob- 
ably due to a desorption of excess oxalate ion. 
This implies a further perfecting of the surface 
which cannot be reliably established by measure- 
ments of dye adsorption. 

It is of interest to note that the reaction rates 
of the 0.1-hour-old precipitates are almost as 
great as those of 0.5-mole silver nitrate under the 

1 T, H. James, J. Am. Chem. Soc. 64, 731 (1942). 

* Individual rate determinations with the previously- 
studied silver chloride seldom varied by more than 2-3 
percent from the mean. Determinations with silver oxalate 
precipitates which had aged at least 17 hours did not vary 
more than 5 percent from the mean, but the deviation in 
the case of precipitates aged only 1-2 hours was sometimes 
considerably greater. It is to be noted, however, that in no 
case were the rates for the 1-2-hour-old precipitates as 

reat as the smallest measured rates for precipitates aged 17 
ours or more. Relative initial pH values for either set of 


aging experiments recorded in Table I did not vary by 
more than 0.02 unit. 
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same conditions (i.e., R=1.10, 1/t=1.0), even 
though the silver-ion concentration supplied by 
the nitrate solution is 8.3X10-*, or 308 times 
that which the oxalate could supply at equi- 
librium. 

Adsorption of gelatin to the silver oxalate pro- 
duces a marked effect upon the reaction rate. 
This is shown by the data plotted in Fig. 1. 
The gelatin curve has the same shape as that 
obtained previously for silver chloride,* and 
shows a definite break. Some decrease in rate is 
obtained beyond the break, but the decrease is 
small and the break is sufficiently well defined to 
warrant the assumption that it corresponds to 
’ the completion of a definite gelatin layer. The 
thickness of this layer is the same for the oxalate, 
chloride, and thiocyanate, as shown by the data 
collected in Table II. 

The dye, 3,3’-diethyl-9-methylthiacarbocya- 
nine, has little effect upon the reaction of the 
oxalate, as shown in Fig. 1. This is in sharp 
contrast to the effect of the dye upon silver 
chloride,* which is considerably greater than that 
of gelatin. Furthermore, this dye produces no 
optical sensitization of the oxalate.* The much 
weaker adsorption of the dye to the oxalate than 
to the chloride serves to emphasize the role which 
the anions play in securing adsorption of cationic 
dyes of this type. The mechanism of adsorption 
of gelatin is evidently quite different, and does 
not depend upon primary adsorption of anions. 
The anionic dye, erythrosin, produces a greater 
effect upon the rate of reaction of silver oxalate 
than the thiacarbocyanine, but not as great an 
effect as gelatin. 

At pH = 10.2, silver oxalate, like silver chloride, 
yields considerable nitrous oxide in its reaction 
with hydroxylamine. Unlike the chloride, how- 


TABLE III. Effect of bromide ions on the oxidation product 


pH = 10.2. 
Br-(m moles) 0.00 0.01 0.02 0.04 0.10 
% N20 60 54.5 47.55 40.5 34.0 


* This — with the observations of S. E. Sheppard and 
W. Vanselow that carbocyanine dyes apparently do not 
sensitize silver oxalate emulsions. 


T. H. JAMES 


4 
60 
50 
7 
40 
30 
20 = 
oo! 002 006 O10 
m moles 


Fic. 2. Effect of Br~ on reaction of 0.5 moles Ag2C20,, 
Agel, and AgCNS. oo, AgeC:0,; A A, AgCNS; xx, 
gCl. 


ever, small amounts of bromide ion do nat reduce 
the nitrous oxide yield to a very low figure, al- 
though some decrease occurs. Data for the 
oxalate are given in Table III, and Fig. 2 com- 
pares the effects of bromide upon the oxalate, 
chloride, and thiocyanate. 

The reaction of hydroxylamine at pH = 10.2- 
10.8 with silver chloride or thiocyanate® to which 
small amounts of bromide ion have been added 
yields chiefly the same product as that obtained 
with pure silver bromide, namely, nitrogen. 
A layer of silver bromide is formed about the 
chloride or thiocyanate particles? which prac- 
tically converts the reaction of these compounds 
to the silver bromide type. In the case of silver 
oxalate, however, the oxide yield is not much 
lower than would be expected for a simple mix- 
ture of bromide and oxalate in the proportions 
indicated. A silver bromide layer evidently is not 
formed about the silver oxalate particles, In 
agreement with this conclusion, Sheppard and 
Vanselow observed that treatment of silver 
oxalate emulsions with very dilute bromide solu- 
tions did not promote dye adsorption or optical 
sensitization. 

2 T. H. James, J. Phys. Chem. (in press). 


13S. E. Sheppard and W. Vanselow, private communi- 
cation. 
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An Empirical Equation for Thermodynamic Properties of Light 
Hydrocarbons and Their Mixtures 


II. Mixtures of Methane, Ethane, Propane, and n-Butane* 


MANSON BENEDICT, GEORGE B. WEBB, AND Louis C. RUBIN 
Petroleum Research Laboratory, The M. W. Kellogg Company, Jersey City, New Jersey 


(Received September 22, 1942) 


An empirical equation is given for the isothermal variation with density of the work content 
of hydrocarbon mixtures in the gaseous or liquid state. From this fundamental equation are 
derived (a) an equation of state, (b) an equation for the fugacity, (c) an equation for the 
isothermal variation of the enthalpy, and (d) an equation for the isothermal variation of the 
entropy. These equations summarize P-V-T-x properties of the gaseous or liquid phase, liquid- 
vapor equilibria, critical properties, latent heats of evaporation, and heats of mixing. The 
equation is a generalization to mixtures of the equation for pure hydrocarbons given in the 
first paper of this series. Methods are proposed for determining the parameters in the equation 
for mixtures from the parameters in the equation for pure hydrocarbons. A comparison is 
made between observed properties of mixtures and those predicted by the equations. These 
comparisons show that the equations provide a satisfactory practical summary of the volumetric 
properties and liquid-vapor equilibria of mixtures of methane, ethane, propane, and n-butane. 


A. INTRODUCTION 


HE object of this paper is to describe an 

empirical equation which summarizes the 
following groups of thermodynamic properties of 
either gaseous or liquid mixtures of light hydro- 
carbons: (a) Pressure-volume-temperature-com- 
position data, (b) Liquid-vapor equilibria, (c) 
Latent heats of evaporation and heats of mixing. 
The equation employed for this purpose is a 
generalization to mixtures of the equation for 
thermodynamic properties of pure hydrocarbons 
given in the first paper! of this series. 

The use of a single equation to summarize 
these thermodynamic properties of mixtures has 
a number of practical advantages. It permits the 
condensation of experimental results obtained at 
a variety of temperatures, pressures, and com- 
positions into a single compact expression. The 
formulation by a single equation of several types 
of experimental measurements, such as P-V-T 
data, heats of mixing and liquid-vapor equilibria, 
is necessarily thermodynamically consistent. 
Calculation of fugacities, partial volumes or 
partial enthalpies, by differentiation with respect 

* Some of the material of this paper was presented at the 
Symposium on “Reactions and Equilibria in Chemical 
Systems under High Pressure” at the 100th meeting of the 
American Chemical Society at Detroit, Michigan, Sep- 


tember 9-13, 1940. 
RO eam Webb, and Rubin, J. Chem. Phys. 8, 334 


to composition, can be carried out more readily 
by an equation than by graphical means. Use of 
the same equation for gaseous and liquid phases 
permits treatment of the critical region, in 
which satisfactory representation by conven- 
tional means is difficult. Finally, the most useful 
feature of an equation for mixtures is that it 
permits reliable prediction of thermodynamic 
properties of mixtures of any composition what- 
ever, a feature of great practical value because 
of the impossibility of determining experimen- 
tally the properties of every mixture of interest. 
Historical 

Because of these useful aspects of a thermo- 
dynamic equation for mixtures, a number of 
rules or equations for predicting such thermo- 
dynamic properties have been proposed in the 
past. Although these former proposals are more 
limited in scope than the equation to be de- 
scribed in this paper, several of them have been 
useful in developing the present equation. These 
will be described briefly. For a more general dis- 
cussion of methods for predicting thermodynamic 
properties of mixtures, reference may be made to 
review articles. 

7 Wien and Harms, Handbuch der Experimental Physik 
(Akademische Verlagsanstalt, Leipzig, 1929), Vol. 8. 


Beattie and Stockmayer, Reports on Progress in Physics 7, 
195 (1940). 
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The essential principle of van der Waals theory 
of the continuity of the gaseous and liquid states 
is the application to both states of the same 
equation for the variation of work content with 
density and composition. The same principle 
has been followed in developing the equation of 
this article. From such an equation for the work 
content the temperatures, pressures, and com- 
positions of coexistent liquid and gaseous 
mixtures may be evaluated. Van der Waals? 
equation of state led to the correct qualitative 


aspects of phase equilibrium between liquid and’ 


vapor, but did not provide an accurate quan- 
titative description, because of its failure to agree 
with P-V-T properties of actual mixtures, 

A number of equations of state for gas mixtures 
more accurate than van der Waals have been 


developed, among which those of Keyes* and_ 


Beattie and Bridgeman‘ are perhaps the best 
known. These equations are suitable for gases up 
to their critical density, but they cannot be used 
for P-V-T properties of gases at higher density or 
for liquids, and cannot therefore be used to 
predict liquid-vapor equilibria. Beattie and 
Ikehara® have given a detailed account of equa- 
tions of state for gas mixtures and have described 
various methods for expressing the dependence 
of parameters in these equations on composition. 

General thermodynamic principles underlying 
the formulation of properties of mixtures have 
been discussed by Gillespie® and by Beattie? who 
showed that all thermodynamic properties of 
mixtures might be computed from certain 
limiting properties of the pure components at 
low pressures and experimental P-V-T data or an 
equation of state for the mixture. 

A more fundamental approach to the predic- 
tion of the’propertiesof mixtures is afforded by the 
statistico-mechanical treatment of J. E. Mayer.® 
In principle, this approach makes possible the 
calculation of an equation of state for gas mix- 
tures from the interaction energies of aggregates 
of molecules making up the mixture. In practice, 

2 J. van der Waals, Die Continuitét des Gasformigen und 
Fliissigen Zustandes (Barth, Leipzig, 1899). 

3 F. G, Keyes, Proc. Nat. Acad. Sci., Wash. 3, 323 (1917). 

‘J. A. Beattie and O, C. Bridgeman, Proc. Am. Acad. 
Sci. 63, 229 (1928). 

5 J. A. Beattie and S. Ikehara, Proc. Am. Acad. Sci. 
64, 127. (1930). 

6 L. J. Gillespie, Phys. Rev. 34, 352 (1929). 


7]. A. Beattie, Phys. Rev. 31, 680 (1928). 
8 J. E. Mayer, J. Phys. Chem. 43, 71 (1939). 
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too little is known of the molecular properties of 
complex substances, such as hydrocarbons, for 
this approach to be of quantitative value. Its 
principal applicability to the present problem is 
in the limitations it imposes on the functional 
form of equations of state for gas mixtures. 

Each of the references cited above has been 
helpful in extending to mixtures the group of 
equations for pure hydrocarbons, given in the 
first paper of this series.' Our treatment resembles 
that of van der Waals,” in that a single equation 
is proposed for both the liquid and gaseous 
phases. It differs from it in that more exact cor- 
respondence with experimental data is obtained 
through use of a more complicated expression, 
which may be regarded as an extension of the 
equations of Beattie and Bridgeman! and Beattie 
and Ikehara.® 


B. EQUATION FOR MIXTURES 
Definition of Terms 


The fundamental equation for thermodynamic 
properties of mixtures makes use of the same 
function as has been used! for pure hydrocar- 
bons, namely, the residual work content A, For 
mixtures, the residual work content is defined as 
the difference between the actual work content A 
of one mole of mixture, and the work content of 
one mole of the mixture in the hypothetical ideal 
gas state at the same temperature, density, and 
composition. The work content is also known as 
the Helmholtz free energy or Gibbs ~-function ; 
it is related to the internal energy E, entropy S, 
and absolute temperature 7, by the equation: 


A=E-TS. (B1) 

A precise definition of the residual work content 
is afforded by the equation: 

A=A-—RT In d—lim (A—RT Ind), (B2) 


where R is the gas constant and d is the molal 
density. Symbols introduced in this section are 
listed in Table I. 

The statistico-mechanical treatment of Mayer 
may be used to show that the limit of Eq. (B2) 
exists and has the following value: 


lim (A—RT In d)=>; x,LRT In x; 
+lim (4;—RT Ind)},_(B3) 
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car- 
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TABLE I. Nomenclature. 


Work content of mixture, per mole 
Residual work content, per mole 
Work content of ith component, per mole 
Molal density 

Internal energy of mixture, per mole 
Free energy of mixture, per mole 
Fugacity of ith component 
Enthalpy of mixture, per mole 
Number of moles of ith component 
Total number of moles (2; 7;) 
Pressure 

Universal gas constant 

Entropy, per mole 

Absolute temperature 

Volume 

Mole fraction of ith component 
Chemical potential of ith component 


~ 


A oe of pure ith component in. hypothetical 
: ideal gaseous state, at unit pressure. 


where A; is the work content of one mole of the 
pure ith component and x; is its mole fraction. 
By combining Eqs. (B2) and (B3), an equation 
for the work content of a mixture is obtained: 


xi_RT In dx; ‘ 

+lim (A;—RTInd)]}. (B4) 
The value of Eq. (B4) is that it segregates the 
work content into the part within the brackets, 
which may be evaluated from properties of the 
pure components at low pressure, and the re- 


sidual work content A, which may be evaluated 
from P-V-T data for the mixture according to 


the equation: 
4 P—RTd 
A= ——wa. (B5) 
@ 


Conversely, the pressure of 2 mixture is related 
to the residual work content by the equation: 


P=RTd+d*(dA/dd)r, (B6) 


The fugacity of the ith component in a mixture 
fi is obtained from the residual work content as 
follows: 


RT In In dRTx;, (B7) 


where n; is the number of moles of the ith com- 
ponent and N=)}°;n;. Since the conditions for 
liquid-vapor equilibria may be expressed solely 
in terms of the pressure, temperature, and fu- 
gacity, it follows that the pressures, temperatures, 
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and compositions of coexistent liquid and vapor 
phases may be evaluated from the residual work 
content alone. 

Other themodynamic properties of mixtures 
may be evaluated from corresponding properties 
of the pure components in the ideal gas state at 
unit pressure and the residual work content 
according to the following equations: 


Energy : 


0A/T 
) + (B8) 
OT 


Enthalpy : 
IT 
--r(—) +P/d—RT 
oT d,z 


+2: xH,°. (B9) 
Entropy : 


S=—(0A/0T)a,z 
xiLR In dRTx;—S;"]. (B10) 
Free Energy : 
F=A+P/d—RT 
+>: x,_RT In dRTx;+F;°]. (B11) 
Chemical Potential : 
(@NA/dn,)y, 
+RT In dRTx;+F;. (B12) 
The properties of the pure components, E,°, F;°, 
H;®, and S;°, may be calculated either statisti- 
cally, from molecular data, or with the third law, 
from calorimetric measurements down to low 
temperatures. This paper is not concerned 
further with these properties of pure com- 
ponents; its primary object is to develop an 
equation expressing the dependence of the 


residual work content of a mixture on density, 
temperature, and composition. 


Development of Equation 


In the first paper of this series,! it was shown 
that the residual work content of a pure hydro- 
carbon varied with temperature and density 
according to the equation: 


A =(BoRT—Ao—Cy/T?)d 
+(bRT —a)d?/2+aad5/5 
cd*7>1—exp (—7d?) exp (— 7d?) 
yd? 


| (B13) 
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TABLE II. Empirical equations for thermodynamic prop- 
erties of mixtures of hydrocarbons. 


Pressure 
P=RTd+(BoRT—Ao— Co/T?)d?+ (6RT —a)d* 


(1+ yd") exp (— rd). (a) 
Fugacity 
RT |n f:=RT In dRTx; 
+[(Bo+Boi)*RT —2(A 0A oi)! —2(CoCoi)#/T? ]d 
+3/2[RT(8b;)! — (a’a;)! Ja? 


3d?(c*c;)'[ 1—exp (— yd?) _ exp (—yd*) 


T? ya? 


—exp (—yd 
Enthalpy 
H=2; x:H®+(BoRT —2A9—4Co/T?)d 


ya? 


exp (—-yd?) 
2 


Entropy 
R |n dRTx;) — (BoR+2Co/T?)d 


S= xi( 
2cd*{ 1—exp (—yd?) _ exp (— vd?) 
—bRd?/2+ Ta [ vd 5 |. (d) 


+d? exp (- (c) 


* This equation assumes i? peg of Bo (Eq. B25). For 
Lorentz combination of Bo (fq. B 26) the coefficient marked with an 
asterisk should be Bao +(Boi)§]3/4. 


It has been found experimentally that if an 
equation of state of a particular form applies to 
each of the components of a mixture, an equation 
of the same form usually applies to every mixture 
of the components, provided they do not react 
chemically. As a hypothesis, it is therefore as- 
sumed that Eq. (B13) expresses the variation 
with density and temperature of the residual 
work content of mixtures, as well as of pure 
hydrocarbons. 

Development of an equation for the residual 
work content of mixtures is then reduced to 
formulating the dependence on composition of 
the eight parameters of Eq. (B13): Bo, Ao, Co, 
b, a, c, a, and y. In this formulation, use is made 
of the following result of the statistico-mechan- 
ical treatment of Mayer.* The residual work 
content of a mixture may be developed in the 
power series : 


WEBB, 


AND RUBIN 


where Q’ and Q” are functions of the tem- 
perature. From this result, it follows that six of 
the parameters of Eq. (B13) must depend on 
composition as follows: 


Bo=L xix ;Bois, 

ang xX ;A 

xix ;Coij, 
b= > ijk, 


ijk 
XX 


ijk 


(B15) 
(B16) 
(B17) 
(B18) 
(B19) 


and 
(B20) 


C= 
ijk 

In these equations, constants with repeated 
suffixes, such as Boi, refer to pure components 
and their values are known. Those with dis- 
similar suffixes, such as Boy2, are a property of the 
mixture, arising from interaction between unlike 
molecules; these may be termed interaction 
constants. 

In this way, development of an equation for 
the residual work content of mixtures may be 
largely reduced to assignment of numerical 
values to the various interaction constants. 
When sufficient experimental data are available, 
they may be evaluated from the data; ordinarily, 
it will be preferable to calculate them from the 
properties of the pure components. The following 
types of equations have been assumed for this 


purpose : 
for b, c, anda: bij = 

for Agand Cy: = 
and for Bo, either : Bo:;= (Bo: +Bo;)/2, 


(B21) 
(B22) 
(B23) 


or: 
(B24) 


Bois=((Boi)*+ (Bos)*]8/8. 
Here, repeated suffixes, such as b;;;, have been 
contracted into a single suffix, such as b;. Choice 
of these equations was made for the following 
reasons. Equations (B22), (B23), and (B24) have 
been used for corresponding parameters in the 
Beattie-Bridgeman equation of state. Equation 
(B23) is termed “‘linear combination of Bo’’ and 
Eq. (B24) is termed ‘Lorentz combination of 
B,.”’ Equation (B21) was chosen so that a simple 
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equation for the fugacity in a mixture would 
result. 

These hypotheses about the interaction con- 
stants, and two additional hypotheses about y 
and a, result in the following equations for the 
dependence on composition of the parameters in 
Eq. (B13): 


Bo=> x:Bo;. (Linear) 

Bo=X (Boi)! + (Boj) }*/8. 

(Lorentz) (B26) 
(B27) 
(B28) 
(B29) 
(B30) 
(B31) 
(B32) 
(B33) 


(B25) 


b=(Li 
xi(ai)* 


From the foregoing hypotheses, consistent 
equations for all thermodynamic properties of 
gaseous or liquid mixtures of hydrocarbons may 
be derived. Equations for some of the more 
common functions are given in Table II. In these 
equations, parameters without the suffix (,), 
such as Bo, refer to the mixture and are evaluated 
from Eqs. (B25) through (B33). Those with the 
suffix 7 refer to pure components. Fractional 
powers of the constants required in the above 
equations for methane, ethane, propane, and 
n-butane are tabulated in Table III. These con- 
stants have been taken from Table II of the 
first paper of this series.'* 


TABLE III. Fractional powers of constants for 
pure hydrocarbons. 


Ethane 
0.0627724 


n-Butane 
0.124361 


Methane Propane 


0.0973130 
0. 


0.0624358 


30.047 44.062 
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TABLE IV. Summarized comparison of equation with 
P-V-T data for methane-n-butane mixtures. 


25.06 50.44 74.90 25.06 5044 74.90 
(b) = (b) (b)  (b) (c) 


0.5-7 :1.25-8 1.25-10 1.25-10 —1.25-8 _1.25-10 1.25-10 


A deviatio 


1.32 
1.52 
1.64 


| 
| 


S25 


eessssssss 
| 


* Beattie, Simard, and Su, J. Am. Chem. Sor. 61, 24, 26 (1939). 

b See reference 9. 

¢ Kvalnes and Gaddy, J. Am. Chem. Soc. 53, 394 (1931); Michels and Nederbragt, 
Physica 3, 569 (1936). 


C. COMPARISON OF EQUATION WITH P-V-T DATA 
FOR MIXTURES OF HYDROCARBONS 


Introduction 


To substantiate the assumptions employed in 
deriving the equations for mixtures, these equa- 
tions will be compared with representative ex- 
amples of the more precise measurements of 
thermodynamic properties of mixtures of me- 
thane, ethane, propane, and n-butane. This 
comparison will be made for both the Lorentz 
and the linear methods of calculating Bo so that 
the more satisfactory of the two methods may be 
found. P-V-T data will be taken up first; then 
liquid-vapor equilibria will be considered. 


Methane-n-Butane 


Probably the most precise and extensive 
measurements of P-V-T properties of mixtures 
of light hydrocarbons are those of Beattie, 
Stockmayer, and Ingersoll? on the following 
mixtures of methane and n-butane: 
Mole percent methane 25.06 50.44 
Mole percent n-butane 74.94 49.56 
These measurements extended from 75 to 300°C 
and from 30 to 300 atmos. A comparison of these 
data with pressures calculated by both the 
Lorentz and linear forms of the equation of state 
is summarized by Table IV. This table shows 


74.90 
25.10 


To =273.13°K. 
R =0.08207 (liter)-(atmos.)/(°K)-(g mole). 


® Beattie, Stockmayer, and Ingersoll, J. Chem. Phys. 9, 
871 (1941). 


|| 
of 
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Mole 
0 
data (a) 
ty 
15) range, 
i 
16) 
17) 75 049 0.17 
100 0.22 «0.13 
125 0.22 
18) 150 032 
200 047 
19) 225 055 — 
250 0.53 
275 0.58 
300 0 0.52 —_ 
20) All 0 0.42 0.38 
nts 
dis- 
the 
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for 
be 
‘ical 
nts. 
ble, 
rily, 
the 
ving 
this 
B21) 
B22) 
B23) 
B24) 
been 
Bo 0.0426000 
hoice (By)! 0.349250 0.397426 | | 0.499147 
wing 1.36198 2.03852 2.62150 3.17564 
0.150233 0.423783 0.712921 0.996410 
have 0.150075 0.223217 0.282310 0.341990 
(a); 0.366924 0.701466 0.982254 1.23471 
the 0.136530 0.319997 0.505277 0.681415 
0.0774600 0.108628 0.148324 0.184391 
tion (a) 0.0499143 a 0.0846781 0.103269 
atl Molecular 
and Weight 16.031 58.078 
on of 
imple 
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TABLE V. Deviations of pressures calculated by the equation of state from observed pressures for gas containing 
50.435 mole percent methane, 49.565 mole percent n-butane. 


Density, 2.001 2.502 3.002 
moles/liter 
Temp., °C 


1.501 


3.502 
Pressure, normal atmospheres 


4.003 4.503 5.004 6.004 7.005 


100 


125 


that the Lorentz form of the equation of state 
for mixtures fits the experimental data almost as 
well as the equation of state fits the data for pure 
methane and pure n-butane. The average devi- 
ation of 0.42 percent for the mixtures, over the 
temperature range 75-300°C, may be compared 
with 0.34 percent for pure n-butane, from 150 to 
300°C, and 0.39 percent for pure methane, from 
75 to 200°C. This agreement is very satisfactory. 
The linear form of the equation of state, with 
an average deviation of 1.56 percent, although 
less accurate than the Lorentz form, is in good 
enough agreement with the data for most prac- 
tical purposes. 

Space does not permit a complete comparison 
of all experimental points with the data. A 
partial comparison is afforded by Table V, for 
the equimolal mixture. Here the published molal 
composition and densities, which were based on 
1941 molecular weights, have been converted to 
1937 molecular weights, the latter having been 
used in obtaining the constants of Table III. 


Except for a uniform tendency for pressures cal- 
culated by the Lorentz form to be slightly lower 
than observed pressures at this composition, 
there are no marked systematic deviations from 
the Lorentz form of the equation of state. At the 
other two compositions, pressures calculated by 
the Lorentz form are slightly high. Pressures 
calculated by the linear form are high at all 
points at every composition. 


Methane-Ethane 


Michels and Nederbragt!® have made careful 
measurements of P-V-T properties of methane- 
ethane mixtures, over a rather narrow tem- 
perature range. In Table VI their observations 
are compared with those predicted by the Lorentz 
and linear forms of the equation of state. This 
comparison shows that for the gaseous phase, 
both the Lorentz and linear forms of the equation 
of state agree with the data within experimental 


10 A. Michels and G. W. Nederbragt, Physica 6, 656 
(1939). 


8.006 9.006 10.007 
Av. % 
Dev. 
= Obs. 31.40 36.24 
% Dev., Linear —0.25 —0.39 0.32 
Lorentz +0.29 +0.30 0.29 
|__| Obs. 34.61 40.25 50.54 59.83 68.37 76.42 84.30 92.18 100.29 118.58 140.30 170.37 214.98 284.73 
% Dev., Linear —0.26 -—0.35 —049 —0.65 -—0.78 -0.98 -—120 —1.51 —189 -—245 -3.50 -—3.67 -—3.85 -—3.20 1.73 
Lorentz +0.26 +0.30 +044 +057 +0.76 +1.02 41.06 +1.02 +058 +0.72 +1.06 +0.91 0.76 : 
150 Obs. 37.73 44.14 56.11 67.21 77.67 87.80 97.84 107.98 118.58 142.04 171.22 210.53 266.72 352.37 
% Dev., Linear —0.37 -041 -046 -—058 -—0.72 -0.88 -—112 -146 -—1.80 -—2.62 -—3.12 -3.08 -—2.87 -—2.53 1.57 
Lorentz +0.13 +0.23 +043 +058 +0.72 +085 +091 +087 +082 +053 +044 +0.70 +091 +099 0.65 
175 Obs. 40.84 47.99 61.57 74.42 86.76 98.90 111.08 123.55 136.62 165.90 202.23 250.71 319.09 
% Dev., Linear —0.37 -042 -049 —0.60 -—0.76 -—0.95 -—122 -155 -190 -—2.57 -—2.95 -—2.83 -—2.45 1.47 
Lorentz +0.15 +0.21 +037 +050 +0.60 +068 +0.68 +062 +0.29 +0.24 +4053 +0.89 0.48 
200—=(t««. Obs. 43.92 51.82 66.99 81.49 95.68 109.77 124.14 138.86 154.51 189.50 232.85 291.10 
% Dev., Linear —0.36 -—045 -—0.69 —113 -—138 -—2.04 -—2.72 -—3.10 —2.68 1.46 
, Lorentz +0.14 +0.23 +037 +038 +043 +042 +041 +0.27 +0.20 -—0.08 -0.18 +0.37 0.29 
. 225 Obs. 46.99 55.61 72.33 88.61 104.70 120.79 137.31 154.50 172.71 213.81 264.80 
% Dev., Linear —0.30 —0.34 -—044 —0.62 -0.78 —1.66 -197 -—253 —2.77 1.26 
Lorentz +0.19 +0.25 +036 +041 +048 +041 +031 40.25 +4014 -—0.07 —0.06 0.27 
250 Obs. 50.02 59.36 77.62 95.57 113.48 131.60 150.29 169.85 190.66 237.81 296.17 
% Dev., Linear -—0.30 -—0.32 -—0.66 -—0.86 -—113 -—144 -—174 -2.04 -2.54 -—2.76 1.29 
Lorentz +0.18 +0.27 +033 +035 +030 +020 -—0.03 -0.22 —0.22 0.23 
. 275 Obs. 53.02 63.08 82.87 102.50 122.23 142.36 163.24 185.18 208.55 261.77 
% Dev., Linear —0.30 -032 -045 —-0.68 -0.91 -—118 -—149 -—180 -—2.11 —2.59 1.18 
Lorentz +0.17 +0.25 +0.33 +030 +0.28 +0.20 +0.09 -—0.04 -—0.19 —0.37 0.22 
300 Obs. 56.02 66.79 88.12 109.45 130.94 176.12 226.40 
% Dev., Linear —0.29 -—030 -—042 -—0.64 —0.95 —1.55 —2.19 0.91 
Lorentz +0.18 +0.27 +0.34 +032 +0.21 —0.02 —0.33 0.24 
Av. % Dev., Linear 0.31 0.36 0.54 0.82 0.96 1.04 1.35 1.64 1.99 2.57 3.03 3.07 2.89 2.87 
; Lorentz 0.19 0.26 0.37 0.43 0.48 0.54 0.46 0.46 0.41 0.38 0.29 0.58 0.95 0.95 
( Over-all Av. % Dev., Linear 1.39 
Lorentz 0.42 
iJ 
1 
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TABLE VI. Deviations of pressures calculated by equation of state from observed pressures for 
mixtures of methane and ethane. 


Phase _ Gas 

Density, moles/liter 3 

Mole’ Temp., 
%CHe °C 


Pressure, normal atmospheres 


Liquid or retrograde gas 
4 5 9 10 11 13 
Av. % 


deviation 


19.96 0 Obs. 
% Dev., Linear 
Lorentz 


25.35 Obs. 
% Dev., Linear 
Lorentz 


Obs. 
% Dev., Linear 
Lorentz 


Obs. 
% Dev., Linear 
Lorentz 


Obs. 
% Dev., Linear 
Lorentz 


Obs. 
% Dev., Linear 
rentz 


Obs, 
% Dev., Linear 
Lorentz 


Obs. 
% Dev., Linear 
Lorentz 


Obs, 
% Dev., Linear 
Lorentz 


Obs. 
% Dev., Linear 
Lorentz 


Obs. 
% Dev., Linear 
Lorentz 


49.84 


Obs. 
% Dev., Linear 
Lorentz 


Linear 
Lorentz 


Average % deviation, 


Linear 
Lorentz 


Over-all av. % deviation, 


57.76 
—3.39 
—2.03 


+0.17 


75.83 
—2.10 
—0.98 


65.89 
—0.49 


69.69 . 
—1.37 
—0.37 


error. At the few points given for the liquid 
phase, the Lorentz form is in better agreement 
with the data. 


Methane-Propane 


Figure 1 shows a graphical comparison of pres- 
sures observed by Sage, Lacey, and Schaafsma™ 
for an equimolal mixture of methane and propane 
with those calculated by the linear form of the 
equation of state. These experimental points are 
more irregular than those considered above, so 
that the comparison with the equation is less 
conclusive. Nevertheless, practical agreement 


a be ae Lacey, and Schaafsma, Ind. Eng. Chem. 26, 214 


with the data is found in both the gaseous and 
liquid regions. 


Conclusion 


The comparisons given above suggest that the 
Lorentz form of the equation of state for mix- 
tures may be used to calculate pressures of 
gaseous or liquid mixtures of light hydrocarbons 
with an average error under 0.5 percent. The 
linear form of the equation of state is in less 
satisfactory agreement with the data, since it 
leads to pressures about 2 percent too high. In 
practical applications, where this error may not 
be serious, the linear form is recommended 
because it is easier to use. 
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ev. 
32 18.59 30.36 
29 40.27 +0.13 1.25 
40.27. +0.20 0.83 
21.04 36.03 46.22 52.90 57.31 
+4033 4035 +019 —0.05 0.22 
40.24 4039 4048 +034 0.33 
49.84 23.38 41.28 $4.93 
40.17 +4031 +0.51 0.33 
40.21 +4036 +0.60 0.39 
47 
40.06 0 19.38 33.42 43.15 
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25.35 21.77 38.78 $2.00 
+0.00 +010 +0.21 0.10 
40.05 4018  +0.35 0.19 
1.26 
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1.29 
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Fic. 1. P-V-T data for equimolal mixture of methane 
and propane. 


D. COMPARISON OF EQUATION WITH LIQUID- 
VAPOR EQUILIBRIA FOR MIXTURES 
OF HYDROCARBONS 


Methane-Propane 


Since the equation of state gives a satisfactory 
representation of P-V-T properties, it is to be 
anticipated that the basic equation for mixtures 
will be equally successful in predicting liquid- 
vapor equilibria. 

In Fig. 2, the observed equilibrium points for 
the methane-propane system at 55°C " are com- 
pared with the border curve predicted by the 
equation, with linear combination of Bo. Except 
for the dew point curve at pressures above 60 
atmospheres, observed and calculated data coin- 
cide within 1 mole percent methane; the agree- 
ment above 60 atmospheres is still quite satis- 
factory. The predicted border curve has the 
correct qualitative characteristics in the region 
of critical phenomena and retrograde condensa- 
tion above 60 atmos. Furthermore, the agree- 
ment between the observed and_ calculated 


3 


8 


PRESSURE, NORMAL ATMOSPHERES 


/ 


— CALCULATED 


O SAGE,LACEY ann SCHAAFSMA 
@ CALCULATED CRITICAL POINT 
® OBSERVEDCRITICAL POINT 


MOLE METHANE 


Fic. 2. Liquid-vapor equilibria for the system methane- 
propane at 55°C. 


critical points is noteworthy : 


Ob- Calcu- 
served lated 


82.00 81.31 


Critical pressure—atmospheres 
43.2 41.4 


Critical composition—mole percent CH, 

The method used in calculating the coordinates 
of a pair of points on the border curve represent- 
ing a liquid and a vapor phase in equilibrium 
makes use of the following principle: Assume a 
liquid composition x;. Find by trial a liquid 
density d”, a gaseous density d°, and a gaseous 
composition y; such that the three conditions for 
liquid-vapor equilibrium, 


P(d”, x1)=P(d%, 91), 
RT In fi(d", x1) =RT In fi(d%, 
RT In f3(d", x1) In y1); 


(D1) 
(D2) 
(D3) 
are satisfied. Then P is the pressure at which 
liquid of composition x; is in equilibrium with 
vapor of composition 


The critical point is calculated by finding by 
trial composition x; and pressure P at the tem- 
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TABLE VII. Distribution coefficients for the methane-n-butane system. 


Observed equilibrium conditions 


Mole fraction 
methane 


Liquid 


Methane 


Calculated 
by equation 
Temp. Pressure 
°F siIb./sq. in. 


Vapor Linear Lorentz 


Linear 


Distribution coefficients 
n-Butane 


Calculated 


% Deviation by equation % Deviation 


Lorentz Linear Lorent 


70 400 
800 
1200 
1600 


TABLE VIII. Distribution coefficients for the ethane-n-butane system. 


Observed equilibrium conditions 


Mole fraction 
ethane 


Liquid 


Ethane 


Calculated 
by equation 
Pressure 


Temp. 
°F Ib./sq.in. 


Linear Lorentz 


Linear 


Distribution coefficients 
n-Butane 


Calculated 


% Deviation by equation % Deviation 


Lorentz Linear Lorentz Linear Lorentz 


0.1748 200 . 3.290 
2.385 


1.602 


8. 
8 

1, 
6. 


5 
8 
2 
58.7 
28.8 
81.8 
26.0 
88.4 
33.4 
74.5 
8 
2 
5 


358 5885 832 38 


Av. % Dev. 


—1.84  —1.84 
—1.59 


—0.34 


0.5122 
0.7042 
0.8743 


perature T (in this case 55°C) such that the 


conditions, 
In fi) 


=0, (D4) 


In fi) 


0x? a 
are satisfied. 

The type of comparison between equation and 
data used for the methane-propane system re- 
quires laborious calculation. In the other mix- 
tures of light hydrocarbons to be examined, a 
simpler, but less exact, method has been used. 
For these other mixtures, distribution coefficients 
(K,) are calculated by the equation: 
G 


fi 
/ di 


Ox, 


(DS) 


(D6) 


where f;”/x; and f:%/y; are evaluated at the 
observed equilibrium compositions. In an exact 
treatment, f;“/x; and f;%/y; should be evaluated 
at the calculated equilibrium compositions, at 
which 


(D7) 


For ideal solutions and in cases where observed 
distribution coefficients equal those calculated by 
Eq. (B24), there will be no difference between 
results of approximate Eq. (D6) and the exact 
Eq. (D7). For the mixtures now to be considered, 
distribution coefficients evaluated from (D6) 
will not differ from the true calculated values by 
more than 2 or 3 percent. Use of Eq. (D6) has 
the practical advantage that it is unnecessary to 
solve by successive approximation for the cal- 
culated compositions. These distribution coef- 
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0.1379 0.8831 6.402 6.203 5.523 +3.11 +13.7 0.1359 0.1291 0.1307 +5.00 +3.83 
0.2804 0.9174 3.269 3.230 2.956 41.19 +4958 0.1149 0.1117 0.1146 42.78 +0.26 
0.4161 0.9178 2.216 2.205 2.074 +6.41 0.1403 0.1392 -0.1441 40.78 
0.5450 0.8942 1.641 1.647 1.589 —0.37  +3.17 0.2318 0.2279 0.2373 +168 
160 400 —-0.0885 0.6233 7.018 6.994 6.252 +0.34 +10.9 0.4142 -2.61 
800 0.2119 0.7451 3.516 3.551 3.269 -100 +703 0.3236 0.3267 0.3326 —0.96 
1200 0.3328 «0.7593 2.309 -114 +442 0.3605 0.3588 0.3681 40.47 
1600 0.4696 0.7285 1.551 1.560 1.523 -0.58 41.81 05118 0.5107 0.5226 40.21 —0.16 
250 800 0.1519 0.4147 2.731 2.706 2.570 40.92 +5.90 0.6906 0.7045 0.7096 —2.01 —2.75 
1200 0.3020 0.4101 1.358 1.330 1.316 43.09 0.8453 0.8548 0.8597 -1.12  —1.70 
P Av. % Dev. 1.10 6.60 1.76 2.20 
(Vapor Ob, 
A 3.142 +1.70 +6.14 0.5029 0.5122 PF 
2.299 +2.61 46.12 0.6932 0.7046 
0.281 1.608 mm 1.574 +0.37. +2.11 0.8713 0.8748 
0.4509 0.896 1.987 1.986 1.933 40.05 +2.72 0.1894 0.1891 0.1873 40.16 +1.11 
0.8075 1 1.791 1.774 1.744 40.95 +262 0.3506 0.3579 0.3557 —208  —-1.46 
0.703 1 1.559 1.562 1.531 -0.19 41.80 0.5409 0.5438 0.5436 -0.54 
0.6577 0.962 1.436 1.450 1.432 40.89 +212 0.1110 0.1087 0.1056 42.07 +4.86 
hane- 0.9195 1.398 1.395 1.381 40.21 +1.22 0.2352 0.2387 0.2332 -1.49 +40.85 
0.8635 1 1.313 1.313 1.305 +0.00 +0.61 0.3988 0.3964 0.3910 40.60 +1.96 , 
0.7695 1 1.170 1.167 1.168 +0.26 +0.17 0.6734 0.6802 0.6724 -1.01 40.15 
0.8220 0.9665 1.176 1.171 1.168 +0.68 0.1882 0.1894 0.1808 +43.93 
0.9395 1.143 1.140 1.137 +0.26 40.52 0.3399 0.3430 0.3306 -0.91 42.74 
0.8765 1.066 1.064 1.066 +0.19 +40.00 0.6938 0.7038 0.6916 -144 40.32 
aicu- 
lated Fe 0.62 2.06 1.13 1.67 
41.4 
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Fic. 3. Mollier chart for a natural gas of the following composition: methane—89.47 mole percent ; ethane—7.91 mole 
percent ; propane—1.92 mole percent ; n-butane—0.70 mole percent. 


ficients are also called K values, and are often 
used in practical calculations of liquid-vapor 
equilibria. 


Methane-n-Butane 


An investigation of liquid-vapor equilibria in 
the system methane-n-butane has been described 
by Sage, Hicks, and Lacey." A comparison of 
their observed liquid-vapor equilibria with the 
predictions of the equation is presented in Table 
VII. The observed distribution coefficients 
quoted in this table are taken from Table I of the 
above authors. The points chosen for comparison 
were selected to include the extreme tempera- 
tures studied experimentally and to cover as 
wide a range of pressure as would permit accurate 
determination of distribution coefficients from 
the observed dew- and bubble-point data. Table 
VII shows that distribution coefficients calcu- 


12 Sage, Hicks, and Lacey, Ind. Eng. Chem. 32, 1085 
(1940). 


lated from the linear form of the equation are in 
very good agreement with the experimental data, 
since the average deviation, 1.10 percent for 
methane and 1.76 percent for n-butane, are of 
the order of the experimental error. Distribution 
coefficients calculated from the Lorentz form of 
the equation are in less satisfactory agreement 
with the observations, since the average devia- 
tions for this form are 6.60 percent for methane 
and 2.20 percent for n-butane. 


Ethane-n-Butane 


Liquid-vapor equilibria in the system ethane- 
n-butane have been reported by Kay." To 
facilitate comparison with the equation, these 
data have been plotted on large-scale tem- 
perature-composition diagrams at pressures of 
200, 400, 600, and 800 Ib./sq. in. and have been 
interpolated to the liquid compositions studied 


13 W. B. Kay, Ind. Eng. Chem. 32, 353 (1940). 
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conditions and the ‘‘observed”’ distribution coef- 
ficients calculated from them are listed in Table 
VIII. In this table a comparison is made between 
the ‘‘observed’”’ distribution coefficients and 
those predicted by the linear and Lorentz forms 
of the equation; the results are similar to those 
found for the methane-n-butane system. The 
predictions of the linear form of the equation are 
in excellent agreement with the data, since the 
average deviations, 0.62 percent for ethane and 
1.13 percent for m-butane, are of the order of the 
experimental error. The Lorentz form of the equa- 
tion is less satisfactory. 


Conclusion 


The examples cited above show that the linear 
form of the equation is more satisfactory than 
the Lorentz form for predicting liquid-vapor 
equilibria. Since the linear form is simpler, and 
since it is fairly satisfactory for calculating 
P-V-T properties, the linear form is preferred in 
most applications of the equation for calculating 
thermodynamic properties of mixtures of hydro- 
carbons. 


E. USE OF EQUATION 


Three examples will be given of the use of the 


‘linear form of the equation to predict thermo- 


dynamic properties of mixtures of hydrocarbons 
for which no experimental data are available. 


Mollier Chart 


In compressing or expanding natural gas, cal- 
culations of pressure and temperature changes 
and input or withdrawal of work and heat are 
effected most conveniently with the aid of a 
Mollier chart, exhibiting the variables, pressure, 
temperature, enthalpy, and entropy. Figure 3 
illustrates such a chart for a natural gas of the 
composition indicated. This chart has been cal- 
culated from the following equation : 


Pressure: Eq. (a) of Table II; 
Enthalpy: Eq. (c) of Table II; 
Entropy: Eq. (d) of Table II. 


The functions H® and S° occurring in these equa- 
tions were evaluated from the limiting specific 
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experimentally. These interpolated equilibrium — 


heat at zero pressure C,° as follows: 


T 
C,"dT, (E1) 
T 
S*= | C,%T/T, (E2) 
r 
where ° 
To=273.13°K. 


It has been assumed that C,° for light hydro- 
carbons is given by the equation : 


C,°=A +BT, 


where A and B are: 


(E3) 


A(L-atmos./g mole-°K) B(L-atmos./g mole-°K?) 


Methane 0.2231 0.000468 1 
Ethane 0.1764 0.001170 
Propane 0.3718 0.001239 
n-Butane 0.4255 0.001766 


Enthalpy-Composition Diagram 


Design of distilling columns operating at high 
pressures is complicated by the fact that dis- 


19.000 


ENTHALPY, BTU/LB.MOLE 


—— TIE LINE 
D-DISTILLATE 
F-FEED 
B-BOTTOMS 


MOLE% ETHANE. 
Fic. 4. Enthalpy-composition diagram for the system 
ethane-propane at 600 Ib./sq. in. 
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Fic. 5. Liquid-vapor equilibria for the system methane- 
ethane-propane at 800 Ib./sq. in. and 131°F. 


tribution coefficients and latent heats vary with 
composition. Despite these complications, col- 
umns separating two-component mixtures can be 
accurately designed with the aid of an enthalpy- 
composition diagram, by the method of Ponchon 
and Savarit.“ Such a diagram for the system 
ethane-propane at 600 pounds per square inch 
absolute is illustrated in Fig. 4. Compositions of 
coexistent phases were calculated in the same 
manner as for Fig. 2, and enthalpies in the same 
manner as for the Mollier chart. The example 
worked out in Fig. 4 is for a column of the fol- 
lowing characteristics : 


Feed Overhead Bottoms 


Temperature, °F 141 82 196 
Mole % Ethane 50 95 5 


Stripping 


11 
+ 13,200 


Enriching 


Number theoretical trays 
Heat input, B.t.u./lb. mole prod. 


7 
— 12,800 


Three-Component System 


Calculation of liquid-vapor equilibrium in 
multicomponent systems is difficult because of 
144 C, S. Robinson and E. R. Gilliland, Elements of Frac- 


tional Distillation (McGraw-Hill Book Company, New 
York, 1939), p. 107. 


the large number of simultaneous equations to 
be solved by trial. To illustrate this type of cal- 
culation, equilibrium compositions of liquid and 
vapor phases of the system methane-ethane- 
propane at 800 Ib./sq. in. and 131°F are shown 
in Fig. 5. 

Conditions for equilibrium in this problem are 
expressed by the following equations: 


P(d®, ¥2) = 800 Ib./sq.in. 

x1, = 800 Ib./sq.in. 
RT In fi(d%, = RT In fi(d”, x1, x2). 
RT In fo(d%, ¥2) = RT In fo(d¥, x1, x2). 
RT In f3(d%, yi, v2) = RT In f(d”, x1, x2). 


(E4) 
(ES) 
(E6) 
(E7) 
(E8) 


The subscripts 1, 2, 3 refer to methane, ethane, 
and propane, respectively. To find the composi- 
tions of a pair of coexistent phases, one of the 
mole fractions, x; for example, is fixed. It is then 
necessary to find by trial the composition of the 
gas phase 91, yo, the density of the gas phase d°, 
the density of the liquid phase d”, and the mole 
fraction of ethane in the liquid x2, such that the 
five equations of condition are satisfied. 

The results of such a calculation are illustrated 
in Fig. 5. The full line of this figure represents the 
border curve computed for this system; the 
dotted lines represent tie lines connecting liquid 
and vapor phases. The occurrence of critical 
phenomena in this system at the specified tem- 
perature and pressure is demonstrated by the 
approach of the tie lines to tangency with the 
border curve as the mole fraction of ethane in 
the system is increased. 
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Letters to the Editor 


HIS section will accept reports of new work, provided 

these are terse and contain few figures, and especially 
few halftone cuts. The Editorial Board will not hold itself 
responsible for opinions expressed by the correspondents. 
Contributions to this section should not exceed 600 words in 
length and must reach the office of the Managing Editor not 
later than the 15th of the month preceding that of the issue in 
which the letter is to appear. No proof will be sent to the 
authors. The usual publication charge ($3.00 per page) will 
not be made and no reprints will be furnished free. 


Velocity of Sound as a Bond Property 


R. T. LAGEMANN AND J. E. Corry 
Emory University, Emory University, Georgia 
November 16, 1942 


ECENTLY in this Journal Rao! has pointed out that 

the function V=v,!/*M/d, where v, is the velocity of 

sound, is a constant characteristic of a pure liquid, and is 

additive and constitutive in that it may be computed by 
summing increments assigned to atoms and bonds. 

An alternative procedure for computing this constant is 
to assume that there is an additivity only of what we 
choose to call ‘‘bond velocities.” Thus for the paraffins we 
have (n—1)(C—C)+(2n+2)(C—H)=V, where the con- 


TABLE I. Bond velocities. 


Bond Bond increment Bond Bend increment 
C-H 95.2 Cc-Cl 230 

4.25 C=C 129 

c-O 34.5 186 
O-H 99.0 


TABLE II. Experimental and calculated molecular sound velocities. 


Molecular sound velocity 


Calculated Calculated 

from bond from Rao’s 

Liquid Experimental increments increments 
Pentane 1160 1159 1160 
Hexane 1356 1354 1355 
Heptane 1545 1549 1550 
Octane 1746 1743 1745 
Methyl acetate 851 830 843 
Ethyl acetate 1037 1025 1038 
Propyl acetate 1211 1220 1233 
Butyl acetate 1408 1414 1428 
Amy] acetate 1598 1609 1623 
Methyl alcohol 421 419 454 
Ethyl alcohol 624 614 649 
Propy] alcohol 806 809 844 
Butyl alcohol 1004 1003 1039 
Amy] alcohol 1198 1198 1234 
Benzene 979 971 945 
Toluene 1170 1166 1140 
Xylene 1362 1360 1335 
Cymene 1731 1750 1725 
Acetone 781 766 769 
Diethylketone 1161 1155 1159 
Methylhexylketone 1720 1739 1744 
Methylene chloride 649 650 649 
Ethylene chloride 846 845 844 
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stant V, which varies from liquid to liquid, may be regarded 
as analogous to the molecular refractivity and termed the 
“molecular sound velocity.”” On this assumption equations 
have been set up for the members of several homologous 
series using the experimental values of V listed by Rao, 
and solved for the bond increments. The method of least 
squares was employed on the paraffins to find the C—C 
and C—H increments which were then used to secure other 
bond velocities from other series. This has resulted in the 
tentative average values given in Table I. 

In Table II are listed the experimental values of V, those 
computed using the atomic and bond values given by Rao, 
and those computed with the bond values of Table I. 
The present assumption yields as good agreement with 
experiment as does Rao’s assumption, and is, moreover, 
more simple. It neglects, of course, single-double bond 
resonance and any variations in bond values in different 
molecules or different series. 

Certain relations of the bond velocities to other bond 
properties have been found. These will be discussed at 
another time. 


1M. Rama Rao, J. Chem. Phys. 9, 682 (1941). 


Quantization of Molecules 


KASIMIR FAJANS 


Department of Chemistry, University of Michigan, 
Ann Arbor, Michigan 


November 12, 1942 


HE difference in the internuclear distance* of He 

(0.74A) and Liz (2.67A) can be understood by assum- 
ing a different quantization of the electron pair binding the 
2H? and 2Li*, respectively. The comparison with the ideal 
unpolar Bohr model led! to the assignment »=1 for He 
and m=2 for Liz for the principal quantum number of 
these electrons. 

The idea of W. Kossel and G. N. Lewis of complete 
noble gas shells and the concept of interaction between 
atomic cores? and electrons prove to be useful for the quan- 
tization of other covalent molecules. 

The Hz molecule is a diatomic analog of the He atom, 
inasmuch as both particles have »=1 and have weak 
external fields. The difference between them is due to the 
splitting of the He** into 2H*. 

Elements of atomic numbers 2 to 6 do not form gaseous 
diatomic molecules at N.T.P. The analogy between phys- 
ical properties of Nz (and CO) and argon has been em- 
phasized by Langmuir.* He proposed a model consisting 
of two N** cores, and two electrons near them, surrounded 
by a (cubical) noble gas configuration of 8 electrons. As 
has been shown‘ on the basis of crystal structure, molar 
volume, and refraction, the isoelectronic particles C2"; 
NC-; Ne; CO; NO*; (O.++?) have outer shells which 
deviate only slightly from the spherical symmetry of the 
noble gases. 

Since Nz has 10 electrons beyond its two K groups, one 
can compare Ne: with Ne by splitting the charge of Ne™+ 
into 2N5*. In accord with the general principle® that the 
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splitting of a nucleus within a given electronic shell leads 
to a loosening of the latter, one finds: the molar refractions 
of He and Hp are 0.5 cc and 2.0 cc; those of Ne and Ne 
are 1.0 cc and 4.0 cc, respectively. 

This loosening of the electronic system is accompanied 
by a strengthening of the external field of the molecule; 
with respect to the boiling point, Hz is nearer to Ne than 
to He, and Ne is nearer to Kr than to Ne. 

The significance of the ten electron shell is also shown 
in a maximum of binding strength (maximum in force 
constant, minimum in distance), e.g., in the series C2, Ne, 
F2; BeO, BO, CO, NO, O2; C2, CoH», 

The above facts and the quantum formula of the neon 
atom 15?2s?2p%, lead one to sharpen Langmuir’s picture of 
Nz by allocating the principal quantum number toa 
group of 2 electrons and m =2 to a group of eight electrons. 
It has been emphasized! that the principal quantum num- 
ber of the electrons nearest to the bonded cores can be 
expected to depend on the spacial extension of the elec- 
tronic system of the latter. The fact that N** is much 
smaller than Li*+ and the internuclear distance in Ne 
(r=1.09A) is also smaller than in Liz (r=2.67A, n=2) 
gives additional support to the quantum formula of N2 : 
KK; 1°28. A subdivision of the eight electron shell has to 
be based on spectroscopic considerations.?° 

While each of the two pairs of K electrons of Ne is 
quantized with respect to one of the two nuclei N**, the 
third pair with n=1 and the 8 electrons with »=2 are 
quantized with respect to the field of both N**. Therefore 
the above does not contradict the Pauli principle which in 
its original form applies to monatomic particles. 

* The data on diatomic molecules used in this and the following two 
letters are from G. Herzberg (see reference 2b). Some of the other data 
are from L. Pauling, The Nature Of The Chemical Bond (Cornell Uni- 
versity Press, Ithaca, New York, 1940), second edition. 

1T. Berlin and K. "Fajans, Je Chem. Phys. 10, 691 (1942). 

2a. The methods of applying quantum theory to molecules de- 
veloped by Heitler, London, Pisa and Slater as well as that of 
Hund, Lennard-Jones, and Mulliken compare the molecules with 
atoms. b. See G. Herzberg, aideraler Spectra and Molecular Structure 
(Prentice-Hall, Inc., New York, 1939). 

31, Langmuir, J. ‘Am. Chem. Soc. 41, 901 (1919). See also W. Kossel, 
Ann. d. Physik 49, 360 (1916) 

4K. Fajans and T. Berlin, Bitfalo Meeting of the American Chemical 
Society, September, 1942. The detailed papers will be published soon. 

5 E.g., the molar refraction increases in the series Ne, HF, H20, 
H3N, HC from 1.0 to 6.5 cc. (See N. Bauer and K. Fajans, J. Am. 


Chem. Soc. December, 1942.) The refraction of CeHs is 10.3 that of 
is 12.9 ce. 


Difficulties in the Valence Bond Theory* 


KasimiR FajANs 


Department of Chemistry, University of Michigan, 
Ann Arbor, Michigan 


November 12, 1942 


HE subdivision of the 10 valence electrons of Ne into 

the groups of 2 and 8 electrons (see III) is not in 
accord with its usual electronic formula : N : : : N : indi- 
cating a triple bond. 

It can be shown that this contradiction adds to the 
difficulties of Kekulé’s valence bond theory and its usual 
electronic interpretation. 

In N2* the force constant is smaller, the internuclear 
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distance (1.117A) is larger than in Ne. Therefore, the de- 
tached tenth electron strengthened the binding of the two 
N*5* cores. One has to conclude that the other 9 electrons 
are also bonding since the ionization process eliminates the 
electron which is most loosely bound and there is (see III, 
footnote 4) a close interrelation between the strength with 
which the electrons are bound and that with which they 
bind the cores. Thus, all ten and not merely six electrons 
take part in the binding of the cores.'? 

The usual valence bond formula of Nz fits into the series 
F-—F, O0=0, N=N which one is inclined to extrapolate to 
C=C. The increase in dissociation energy into atoms from 
F, (65 kcal.) to Ne (170 kcal.) could be considered as a 
further support of these formulae, since the energy of the 
carbon-carbon bond increases in the series single, double, 
triple bond.* From this point of view one would expect that 
the bond C=C is still stronger than N=N or —C=C-. 
Contrary to that, the dissociation energy of C2 into atoms 
is 83 kcal., ie., of the order of magnitude of the single 
bond only. Moreover, C2 has a very strong external field 
and thus has no noble gas character. One has also to recall 


that the electronic formula :O::O: has already been 
disproved by the paramagnetism of O2 and that other 


formulae have been proposed for it, e.g., :0:0: by G. N. 


Lewis.‘ 

Thus among the valence bond formulae of the mentioned 
four diatomic molecules, °F :F: is the only one which is 
not in disagreement with the behavior of these elementary 
substances. 

The usual electronic formula of Nz and many others 
based on the idea of completion of a noble gas shell by 
sharing electrons also encounter the following difficulty. 
Only a limited number of atoms or radicals can assume a 
noble gas configuration by an exothermic process: H, F, 
Cl, Br, I, OH, CN, and perhaps some of the polyatomic 
radicals. 

The oxygen ion O~, however, is unstable in the free 
gaseous state and its electronic shell has to be stabilized, 
‘e.g., by H+ in OH™ or by other cations (V). This is also 
true for N*-, C+, or B5-. Thus, in a symmetrical molecule 
Ne one could assume a temporary polarity (II) connected 
with a completion of one N* stabilized by N** from the 
other atom. The above consideration contradicts, however, 
the assumption that in O2, Nz, C2 both atoms at the same 
time complete their octets. 

The contention (see Pauling, III*) that many molecules, 
even the simplest, are the result of the resonance between 
a number of forms with different combinations of single, 
double, and other bond types can be considered as a strong 
criticism of the valence bond theory. For it means the 
sacrifice of the initial aim of the theory to represent one 
experimentally homogeneous molecular species by one for- 
mula. 

The quantum formula of Ne, KK; 1228, (see III) is free 
of these difficulties. It considers all ten electrons as com- 
mon to both cores and resembles in this respect the method 
of molecular orbitals used (III, 2a) in the discussion of 
spectral data. There are, however, many types of sub- 
stances to which neither a valence bond formula nor an 
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electronic shell common to the whole molecule can be 
assigned. They are discussed in the following letter. 


* This is communication IV on “Electronic Structure of Molecules.” 
The other communications in the J. Chem. Phys. will be quoted by 
the following Roman numbers: I. K. Fajans and N. Bauer, 10, 410 
(1942); II. T. Berlin and K. Fajans, 10, 691, 1942; III. K. Fajans, pre- 
ceding letter; V. K. Fajans, following letter. ; 

1 This bonding has to be considered as due to the electron density 
present between the cores. As has been shown (III, 4) it is possible on 
this basis to interrelate properties of molecules which do not differ 
with respect to quantization. 

2It should be noted that considering the molecule as the result of 
the interaction between atoms and not between cores and electrons, it 
appears to be possible (III, 2b) to correlate the spectroscopic properties 
of Ne with the view that eight of the initial electrons of the atoms 
acted as bonding, two as antibonding, and that in the balance the bond 
between the atoms is due to six electrons. 

3K. Fajans, Ber. d. D. Chem. Ges. 53, 643 (1920); 55, 2826 (1922). 

4G. N. Lewis, Chem. Rev. 1, 231 (1924). 
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N any attempt to reach a uniform point of view con- 

cerning the electronic structure of molecules, one has 

to consider the electric interaction and the quantization. 

Using a few characteristic examples we shall combine 

these two factors into a principle which seems to have 
broad application to chemical binding. 

The oxide MgO can be considered as ionic. However, 
O= in the free state is unstable with respect to O~ and e~ (I) 
and, therefore, the principal quantum number of its eighth 
electron is m= ©. When combined with Mgt? the value of 
n is reduced, in all probability, to n=2. This process re- 
quires energy and, thus, the neon state of O= has to be 
stabilized by the binding process. 

This point of view can be applied to an electron pair. 
Left to itself the pair would dissociate. However, it can be 
stabilized by positive charges, in H~, He, Lit, as well as 
in He, Lig and represents a spacial distribution of negative 
charge. Depending on the field strength of the positive 
charges, the pair can be expected to assume a larger or 
smaller size with unchanged principal quantum number n, 
but can also differ in quantization as has been concluded 
(II) for He and Lis. 

One characteristic example of this analogy between O- 
and an electron pair will be given here. The P atoms in 
P,Ow form a tetrahedron, 6 of the O atoms are situated 
near the middle of its 6 edges, and each of the remaining 4 
is attached directly to one P atom. The P-O distances! 
for the 4 oxygens are shorter than those for the 6 oxygens, 
which can be easily understood if one derives the molecule 
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from P5+ and O~. The one-sided polarization leads to 
shorter distances than a many-sided polarization.2 The 
structure of P,O¢ is similar to that of P,Ow suggesting 
that the 40~ have been replaced by 4 unshared electron 
pairs. Finally, in Ps, which also has a tetrahedral structure, 
the 60™ of P,O¢ can be considered to be replaced by elec- 
tron pairs, each pair shared by two P**. 

Coordination compounds are characterized by the bind- 
ing of stable molecules or of radicals stabilized by the 
binding process.| The remarkable group of isoelectronic 
substances Ni(CO),, Co(CO)3(NO), and Fe(CO)s(NO)s is 
diamagnetic. The C—O distance (1.15+0.03A)* in these 
compounds differs slightly from that in free CO (1.128A), 
the N-O distance (1.11+0.04A) is near that in NOt+ 
(1.073A). We suggest that the compounds contain de- 
formed CO and NO, both of which have great stability 
as particles with ten electrons (III). Accordingly, the 
metals are present in these compounds in the isoelectronic 
form Ni, Co~, and Fe™, having the electron configuration 
K?L°(M-+ N)'8, In the normal stable paramagnetic state of 
the Ni atom, the electrons beyond the LZ shell are sub- 
divided into 3s?3p%3d* and 4s*. The transfer of the two N 
electrons into the M shell would complete it and make it 
diamagnetic. We therefore conclude that this state of 
higher energy is stabilized by the energy of binding be- 
tween particles with 10 and 18 electrons. 

The boron hydrides can be fitted into the above scheme, 
e.g., BsHio can be formulated as BH3-BeH,4-BHs. The 
radical BH, is a member of the isoelectronic series No, 
HCN, HCCH, H2BBH: in which the successive splitting 
of the nucleus (III) leads to the loosening of the system. 
Similarly, BsH» is BsH¢- BHs, the former of the two radicals 
having an 18 electron shell.* ; 

The above offers examples of stabilization of radicals 
with various numbers of electrons, 2(2e-), 6(BHs), 8(O™), 
10(B2H,), 18(Ni) within ionic, covalent, and coordination 
compounds. In papers to follow this binding principle will 
be extended to many other cases. 

I appreciate the valuable discussions with Professor 
Norman Bauer, University of New Hampshire, and Mr. 
Theodore Berlin. 


* This is communication V on ‘Electronic Structure of Molecules.’’ 
See communication IV, preceding letter. 
_ IL. R. Maxwell, J. Opt. Soc. Am. 30, 374 (1940). As will be shown 
in a detailed paper, the gradation of the distances in Ps, PsOs, PsOu0 
can be understood from the above point of view. 

2K, Fajans, J. Chem. Phys. 9, 378 (1941). 

+ The binding is due to dipoles and polarizability. 
( 9a) O. Brockway and J. S. Anderson, Trans. Faraday Soc. 33, 1233 

_ ‘The proposed structures are in agreement with the electron diffrac- 

woo of S. H. Bauer and L. Pauling, J. Am. Chem. Soc. 58, 2403 
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